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FOREWORD 


The origins of this series were a number of discussions in the Education Com- 
mittee and in the Council of the International Union of Pure and Applied 
Biophysics (TUPAB). The subject of the discussions was the writing of a text- 
book in biophysics; the driving force behind the talks was Professor Aharon 
Katchalsky, first while he was president of the Union, and later as the honor- 
ary Vice-president. 

As discussions progressed, the concept of a unified text was gradually 
teplaced by that of a series of short inexpensive volumes, each devoted to a 
single topic. It was felt that this format would be more flexible and more 
suitable in light of the rapid advances in many areas of biophysics at present. 
Instructors can use the volumes in various combinations according to the 
needs of their courses; new volumes can be issued as new fields become 
important and as current texts become obsolete. 

The International Union of Pure and Applied Biophysics was motivated to 
participate in the publication of such a series for two reasons. First, the 
Union is in a position to give advice on the need for texts in various areas. 
Second, and even more important, it can help in the search for authors who 
lave both the specific scientific background and the breadth of vision needed 
to organize the knowledge in their fields in a useful and lasting way. 

The texts are designed for students in the last years of the standard univer- 
sity curriculum and for Ph.D. and M.D. candidates taking advanced courses. 
They should also provide a suitable introduction for someone about to begin 
tesearch in a particular field of biophysics. The Union is pleased to collabo- 
fate with the Cambridge University Press in making these texts available to 
students and scientists throughout the world. 


Franklin Hutchinson, Yale University 
Watson Fuller, University of Keele 
Lorin J, Mullins, University of Maryland 
Editors 


PREFACE 


initial reaction upon being invited by the Editorial Committee of the 

iernational Union of Pure and Applied Biophysics to write a monograph on 

ciples of membrane transport was a mixture of flattery and terror. I was 
rally flattered by the invitation but terrified by the seeming enormity of 
task and the effect that acquiescence would have on my already over- 

idened schedule. When this initial reaction gave way to calm objectivity, I 

d myself: With so many books appearing annually on the subject of 

brane transport, is there a justifiable need for yet another? 

_ My answer to this question, which was confirmed y of my col- 
s, was that there still is a need for the type of ‘ener would have 
rately longed for in 1959 when, after completing sidency in In- 

mal Medicine, I joined the Biophysical Laboratories of the Harvard Medical 
as an initiate in the field of membrane transport; my contacts with 
rgraduates, graduate students, postdoctoral fellows, and junior scientists 
ting the past two decades have reinforced this feeling. 
In short, this monograph is intended for the “initiate.” Its goal is to intro- 
iiuee some of the basic principles that govern solute and water transport 
weross membranes and to emphasize the foundations and reasoning that 
wnderlie these principles in a relatively brief and readable form. It is not in- 
tended to supplant the more advanced and comprehensive treatments that 
be found in a number of books, chapters, and review articles, but rather 
to provide the readers with an appreciation of basic concepts and approaches 
that will ease their way into the more advanced literature. In order to achieve 
this goal, I have attempted to develop this monograph from “first principles” 
guming only that the reader has had some previous exposure to elementary 
differential and integral calculus, chemistry, and physics. All final expressions 
wre derived, inasmuch as possible, step by step, and particular care has been 
taken to see to it that a consistent set of conventions is carefully defined and 
adhered to throughout. There are a number of basic points that are repeated 
throughout this monograph; this may irritate the more advanced reader but, 
ut the same time, could be valuable to the beginner. In any event, this intro- 
duetory text is designed to teach and 1 am convinced that repetition is an 
fesential ingredient of pedagogy. 

xi 


xii Preface 


If I have succeeded in achieving my intended goal, a large measure of the 
success is owed to my teachers, colleagues, and students. 

I am particularly indebted to: 

Professor Arthur K. Solomon, who opened the door to my scientific career 
and who permitted (urged) me to teach the first offering of a course entitled 
“Bioenergetics and Membrane Transport” (Biology 119) at Harvard College 
in 1967. This was a memorable experience which forced me to come to grips 
with theoretical principles from the vantage point of a teacher rather thana 
user. Much of the organization and many of the developments of this text are 
based on my lecture notes for that course. 

The late Professor Peter F. Curran, a long-time combination teacher- 
collaborator-friend, who was exceptionally capable of employing theory to 
provide a logical, coherent, and concise “picture” of experimental data. 

And, the late Professor Aharon Katchalsky, a profound thinker, a gifted 
educator, and a constant source of inspiration; one cannot lavish too much 
praise upon his personal and intellectual attributes. 

The tragic and untimely deaths of Peter Curran and Aharon Katchalsky, 
within a brief span of time, came as terrible blows. I am privileged to have 
shared their friendship and collegiality, and owe them both a great deal. 

In addition, I am grateful to Drs. Alvin Essig, Raymond Frizzell, and 
Stephen Thompson for their careful reading of this manuscript and their 
corrections, cogent suggestions, and critical comments. 

I wish to thank my secretary, Susan Frizzel, who nurtured this effort from 
first draft to final product and am eternally indebted to her dog, Fagin, for 
not having completely devoured Chapter 4 during an unguarded moment. (I 
reject my son’s suggestion that he found it “too dry”!) 

Finally, I owe a special brand of gratitude to my wife, Harriet, and my 
children, Jeffery and Kenneth, who have always provided me with the 
warmth, encouragement, and peace of mind to enable me to “do my thing.” 
My debt to them and my parents can never be repaid. 


S.GS. 


1 Some basic principles of 
thermodynamics: The relations 
between flows and forces 


The formal description of transport across membranes rests firmly on the 
unassailable citadel of thermodynamics, that branch of science concerned 
with the changes that accompany the transition of a system (a defined por- 
tion of the universe selected for study) from one set of conditions (state) to 
another.’ In this chapter we will develop the basic principles that deal with 
the transitions of a system from one state (initial) to another (final), which 
result from the transfer of matter within the system.” 


11. Properties of state 


Any macroscopic system at equilibrium can be described com- 
pletely by an equation of state, which relates a set of parameters or properties 
that are functions only of the current state of the system and are independent 
of its past history. It follows that any change in a property of state is uniquely 
determined by the initial and final conditions of the system and is indepen- 
dent of the path (or mechanism) of transition between these twe states. 

A property of state can be formally defined as follows. If Z is a function of 
other properties of state, for example, Z = f(U, V, W,...) and if 


dZ=L-dU+M-dV+N-dW+--- 


where L = (8Z/dU)y, Wawa? M= (82/8 Vu,w,...5 N= (0Z/OW)u, Vosc0% and so 
on, then @Z is an exact (or total) differential of f(U, V, W,...) and Z isa 
property of state.> The reason, of course, is that if these relations hold, then 
1+ dU, M- dV, N- dW, and so on, represent the individual contributions of 
ihe changes in U, V, W, . . . to the total change in Z and, since addition is a 
eommutative operation, the total change in Z does not depend upon the 
order or sequence of the changes in U, V, W, ... . For example, referring to 
Figure |.1, if a system undergoes a change from state A to state C asa result 
of changes in properties U and V, the change in Z will be the same if the 
thange in U preceded the change in V (pathway 1) or if both U and V change 
umultaneously (pathway 2), that is, 
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C (Ug Vor Wy --) 


<-— 


B(Us,Vp,Wa---) 
A(Up,VasWa---) 


l= 
Figure 1.1. U,V, W,... are properties of state. If a system under- 
goes a transition from an initial state A toa final state C as a result 
of changes in U and V, the properties of the final state are indepen- 
dent of the path. 


AZ=Zc- Zz =f(Uc, Veo, W,,.. )-fUs, V4, Wa,- A ) 
= [Zp - Z4)+ Zc - Zz)] 

Finally, the fact that dZ is an exact differential means that there will be no 
change in Z if the system is displaced from state A and then restored to state 
A regardless of the paths (mechanisms and durations) of displacement and 
restoration (e.g., referring to Figure 1.1, A to B and back, A to C and back, 
A to B to C and back, A to C to B and back); that is, the closed integral, 
§ dZ = 0 (Green’s theorem). 

The state properties of a homogeneous system can be divided into two 
categories, namely, intensive properties and exfensive properties. Intensive 
properties are those characteristic of the total system, as well as any arbi- 
trarily selected part of the system, and they are nonadditiye; for example, 
pressure (P), temperature (7), concentration (c), and electrical potential (); 
these properties reflect the “intensity” of the system. The extensive proper- 
ties of a homogeneous system are measures of the size or content of the 
system and are additive; for example, volume (V), number of moles (7). 
Obviously, the total volume of a system is the sum of the volumes of all sub- 


Table 1.1. Conjugate properties of state 


Intensive property Extensive property 


Pressure (P) Volume (V) 
Electrical potential (w) Charge (e) 
Temperature (T) Entropy (S) 
Chemical potential (1) Moles (7) 


jons (parts) of the system whereas the temperature of these subdivisions 
the same. Further, it is obvious that the ratio of any two extensive prop- 
is an intensive property (e.g., n/V = c), and it can be shown that any 
tensive property can be expressed as a function of all other intensive prop- 
ies of asystem (.g., the perfect gas law PV=nRT can be written P= cRT). 
Vivery intensive property is paired with a conjugate extensive property and 
difference in an intensive property is a driving force for the flow or displace- 
nt of its conjugate extensive property. The conjugate properties of state 
h which we will be concemed are listed in Table 1.1. Clearly, a difference 
pressure (AP) is a driving force for the displacement of volume (V); a dif- 
fence in electrical potential (Ay) is a driving force for the flow of charge 
, a8 we will see, a difference in temperature (AT) is a driving force for the 
placement of entropy (S); and for the moment let us define a driving force 
the flow of uncharged matter () as a difference in chemical potential 
). An explicit expression for the chemical potential will be derived below. 
ther, for reasons that will become evident, we refer to a difference in an 
nsive property as the conjugate driving force for the flow of its paired 
ugate) extensive property. Finally, it should be noted that the product 
an intensive property and its conjugate extensive property has units of 
tgy or work. Thus, PdV is the work that must be performed to displace 
amount of volume (dV) against a pressure (P); W de is the work that must 
performed to bring an amount of charge (de) from a region of zero elec- 
potential (ground state) to a region where the electrical potential is p; 
d as we will see, u,dn; is the work that must be performed to transfer a 
mber of moles of an uncharged substance (dn;) from a region where the 
ical potential of 7 is zero to one where it is y;. 


The first and second laws of thermodynamics 


The entire edifice of thermodynamics is founded on two laws or 


watulates that simply summarize a vast amount of human experience. The 


t law, often referred to as the law of conservation of energy, states that 
total internal energy of a system (Z) can only change as a result of the 


‘1 or loss of heat (Q) and/or performance of work (W) on or by the system, 
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—w 


dw= plv -pde “E pide 
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that is, 


dE =dQ- dw (1.1) 


where dQ is positive when heat is gained by the system and dW is positive 
when work is performed by the system on its surroundings. # is an extensive 
property of state determined by the total potential and kinetic energy of the 
components within the system, and dé is an exact differential of an equation 
of state. However, Q and W are not properties of state inasmuch as the 
amount of heat gained and work performed by a system when it changes 
from one state to another depends upon the nature (pathway, mechanism) 
of the transition. 

For a system that is capable of exchanging matter with its surroundings, 
equation (1.1) can be written more explicitly as follows: 


dE =dQ- PdV+wW det >. ny; dn, (1.2) 
i 

where P dV is positive when the system performs work by expanding (dV) 

against an external pressure P; y de is positive when an amount of charge de 

is transferred into a system whose electrical potential is y; and yz; dn; is posi- 

tive when an amount (moles) of matter dn; is transferred into a system where 

the chemical potential of i is yz; (see note 9). 

The first law is a “bookkeeper’s delight’’; it states that a change in the in- 
ternal energy of a system must be the sum of the exchange of heat with its 
surroundings and the work performed on or by its surroundings. Since the 
heat gained or lost by the system must be equal to the heat Jost or gained by 
the surroundings, and since the energy lost by a system as a result of the 
performance of work on its surroundings must be equal to the energy gained 
by the surroundings, the internal energy of the universe is conserved; E is 
simply transferred from one system to another but the total remains con- 
stant. However, the first law provides no way of distinguishing between 
spontaneous or natural processes and processes that can never take place 
spontaneously; according to this law, all processes are equally possible provid- 
ing they obey equation (1.1). 

The second law of thermodynamics is the culmination of many efforts to 


describe the fact that all natural or spontaneous processes take place only 


in the direction toward equilibrium and come to a halt when this time- 
independent state is reached! Work must be performed on a system to dis- 


place it from equilibrium; such a displacement cannot occur spontaneously 
and, for this reason, spontaneous processes are referred to as “irreversible.” 
This monumental achievement was accomplished by introducing the con- 
cepts of reversible processes, namely, processes that take place infinitely 
slowly so that the transition from one state to another may be treated as if 
the system passes through an infinite series of equilibrium states. For such an 
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d process, it can be shown that, although dQ is not an exact differen- 
the ratio (dQ/T) is an exact differential. Thus, a new property of state 
, the entropy (S), which for a reversible process is defined by the 
n 


dS = dQ/T (1.3) 


§ is a property of state, a change in S is dependent only on the initial 
final states of the system and is independent of the intervening path. 
for a given set of initial and final conditions, dS is the same if the tran- 
between these states is the result of a spontaneous (irreversible) process 
wn idealized reversible process. However, for a given transition, the heat 
1 the surroundings is greater (or the amount of heat taken up by the 
from its surroundings is less) when the process is irreversible than 
it is reversible.* Thus, for an irreversible process 


ds > dQ/T (1.4) 


fed another way, if the same change in state is carried out reversibly and 
rsibly 


dS = (dQ/ rey > (dQ/ iy) ee 


1 is, the heat that is actually absorbed in any real process is less than what 
have been absorbed had the process been reversible. The difference 
N (dQ/T),cy and (dQ/T)jnrey was referred to by Clausius as “the un- 
nsated heat.” 

ations (1.3) and (1.4) for reversible and irreversible processes, respec- 
, comprise a formal statement of the second law of thermodynamics. 
The direction of a spontaneous or irreversible process must satisfy the in- 
ty given in equation (1.4). Since S is an extensive property of state, dS 
he divided into two parts, that is, 


dS =d,S+d,S 


dS represents the change in entropy of a system caused by exchange 
eat with the surroundings and d,S represents the internal production or 
tion of entropy (the “uncompensated heat”). For an irreversible process 
N+ a\S>dQ/T, where d,S =dQ/T, so that djS > 0; for a reversible pro- 
d,S + dS = dQ/T, so that d,S = 0. Thus, an irreversible process can be 
d as one that leads to an internal production of entropy that cannot be 
nted for by an exchange of heat with the surroundings. Further, it 
ild be noted that for an irreversible process, d,S and dS can be positive 
fiegative depending upon whether heat is gained or lost by the system, but 
is always positive. 

Vinally, since £ is a property of state, dE is independent of path and fora 
change in state is the same whether the transition is the result of a 


fe (3 a, EAE Boy 
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reversible or an irreversible process. Thus, we can substitute 7 dS for dQ in 
equation (1-2) and obtain 


dE =TdS-PdV+yWde+ > y; dn; (1.5) 


From this equation, it is clear that S is the extensive conjugate of T in the dis- 
placement or flow of thermal energy.* 

Equation (1.5) is often referred to as the “‘Gibbs equation,” and in its most 
general form it includes changes in all possible extensive properties and 
equates dE to the sum of the products of intensive properties and changes in 
their conjugate extensive properties, that is, 


dE = >~ [(Intensive property); X (Change in conjugate extensive property);] 
i 


13. The Gibbs free energy and the concept of “useful work” 


Although the entire fabric of classical thermodynamics can be woven 
out of the first and second laws without introducing any additional properties 
of state, equation (1.5) does not provide the most convenient or illuminating 
framework for the thermodynamic description of processes that take place 
under conditions of constant temperature and pressure. To remedy this short- 
coming, Gibbs introduced the free-energy function G, which is defined as 


G=E+PV-TS (1.6) 


Since £, P, V, T, and S are properties of state, G must also be a property of 
state, so that dG is determined entirely by the initial and final conditions of a 
system and is independent of the intervening path. From equation (1.6) 


dG=dE+PdV+VdP-TdS- SdaT (1.7) 
and when 7 and P are constant 
(dG)r p=dE+PdV-TdS (1.8) 


We now (arbitrarily) divide the total work done by (or on) a system into 
“pressure X volume work” and “useful work” so that 


dW=PdV+dw' (1.9) 


where dW’ represents the useful work and, as before, is positive when the 
system performs work on its surroundings and is negative when work is per- 
formed on the system.® 

Substituting equation (1.9) into equation (1.1), we obtain 


dE =dQ-PdV- dw: (1.10) 
and substituting equation (1.10) into equation (1.8) we obtain 
(dG)r,p=dQ- TdS-dw' (1.11) 
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aS = dS + d,S = (dQ/T) + 4,5, 
T dS =dQ+T4d,S 
1 from equation (1.11) we obtain 


(dG), p=Td,S + dWw' (1.12) 


4) If a system at constant T and P performs useful work on its surround- 
by means of a reversible process, d;S = 0 and the decrease in the Gibbs 
# energy of the system (-dG) is equal to the useful work performed. If the 
is irreversible, d;S > 0 and the decrease in G is greater than the useful 
performed. 
h) If a spontaneous process occurs in a closed system (one that cannot 
ie matter with its environment), dW'=0 and the decrease in free 
of the system is reflected only by an increase in internal entropy; 
1 ls, free energy is completely dissipated. 
» considerations provide some insight into the meanings of internal 
y, free energy, and entropy. Consider for moment an irreversible 
that takes place in a closed system (e.g., the diffusion of a solute from 
on of higher concentration to one of lower concentration). Clearly, the 
] énergy of the system will not be affected since dQ and dW are both 
_ However, the free energy of the system will decrease and reach zero 
equilibrium is achieved. Thus, the potential ability of the system to 
work is dissipated at the expense of the production of entropy; for 
feason free energy functions are frequently referred to as “work func- 
" The same holds for a system that can perform work on its surround- 
Vor any real process the decrease in free energy is greater than the 
unt of useful or recoverable work performed. Thus, the free energy of the 
iiverse has decreased but the internal energy is conserved. 
Why? At the outset it should be stressed that thermodynamics is not based 
‘Wy particular model of matter, and the functions £ and S were intro- 
4 44 properties of state that simply permit a consistent description of the 
tics and direction of natural events. We cannot search for a more pro- 
1 interpretation of these functions in classical thermodynamics. How- 
, let us depart from this secure axiomatic approach and consider a closed 
i subdivided into two compartments by a porous barrier with one com- 
Went containing a monoatomic gas and the other initially void. Clearly, 
gms will diffuse (expand) into the initially empty compartment and equi- 
iin will be achieved when the distribution of atoms is random (i.e., equal 
sentrations in both compartments). Now, the internal energy of this sys- 
4) simply the sum of the kinetic, translational energies of the individual 
and is given by E=n-+m- 07/2 where n is the number of atoms, m is 
wtomic mass and 0? is the mean square velocity of each atom that is only 
ndent on 7, It should be noted that although 0 is a vector, 0? is a scalar, 


F 
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so that while EF embodies the concept of motion, it is divorced from the 
notion of direction. F would ve the same if at any instant every atom were 
moving in the same direction #$S when the movements are directed randomly. 
In contrast, the entropy function is a measure of the randomness of the sys- 
tem. An increase in entropy or an increase in randomness may be defined as 
an increase in the number of p ossible configurations or complexions of a sys- 
tem, all of which are equally likely and accessible; the state of equilibrium 
represents the most random or most probable configuration that the system 
can achieve. Thus, the second Jaw states that the direction of all spontaneous 
processes is from a less random: (less probable) to a more random (more prob- 
able) state. Clearly, with an increase in the randomness of a system, the 
ability of the system to convert its internal energy to a directed effort, or 
useful work, decreases,” 

At the risk of belaboring this point, let us close this discussion with a trivial 
example. When a heavy weight is dropped onto a metal plate in a vacuum, the 
kinetic energy of the falling weight is completely converted into thermal 
energy, which in tum is reflected by an increase in the velocity of random 
motion of particles (entropy)- If the motion were not random, but instead 
if all of the particles could be harnessed to simultaneously move upward 
against the weight, the weight would rise and a spontaneous process will have 
been reversed without outside intervention. The second Jaw infers that the 
probability of this occurrence is so low that this event may be considered 
impossible .® 


1.4 The electrochemical potential 


We now employ the Gibbs free energy function to derive an expres- 
sion for the electrochemical potential of a substance i. As we will see, the 
electrochemical potential is an intensive property of a system, and a differ- 
« ~~ ence in electrochemical potential is the conjugate driving force for the diffu- 

rainael flow of matter. This property will play a central role in all of our 

subsequent considerations so that the importance of understanding its mean- 
ing cannot be overstated. 

Let us start by considering the effect of reversibly transferring matter into 

a system on G. Inasmuch as 


G=E+PV-TS 
dE = dQ - Pav - aw’ 
and for a reversible process 


dQ=TdS 


(1.6) 
(1.10) 


it is a simple matter to show that 


dG =-SdT+Vap- dw’ (1.13) 
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dW’ is the reversible work of transfer of matter into the system given 

y Hy dn, + p de [equation (1.5)]. 

G ts a property of state and, as such, is a function of other state 
§ such as 7, P, and composition, we can express dG as an exact 

tial of these other properties of state as follows: 


dG = (32) dT + (2) dP 
oT P,nj,nj,... oP T,nj,nj,..- 


dG 


+> (— dn; 1.14 
>= iad 


"j,Mj,... Tepresent the number of moles of i,j,... in the system. 
tion of equation (1.14) with equation (1.13) discloses that 
ae) (22) 
— =p; || {-— =V 1.15 
(i P,mysnj,... oP T,nj,nj,... a ) 
aG 
-dw'=>. (32) dn; 1.16 
7 \8m/7,P,nj,... aise 
dW'= > u, dn, + p de 
‘ 
de = 2; F dn, 
#, is the valence of i and F is the Faraday, it follows that 
AW! =X (uy + 21 Fp) dn, 
i 
define the electrochemical potential of i as 
Ay urt2 Fy (1.17) 
dW'e ¥ By dn, (1.18) 
i 


, from equation (1.16) it is clear that @;=(0G;,/dn;) or the partial 
ree energy of i and represents the increase in free energy of the system 

from the addition of one mole of i when T, P, and the composition 
other substances ny .., are constant. In other words, it is the work that 
be performed on the system in order to introduce 1 mole of / when all 
Properties of state are maintained constant. Clearly, Ajj is the con- 
Uriving force for the flow of 1, 
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It should be noted that the electrochemical potential ff; combines the 
partial molal free energy change resulting from a change in the chemical com- 
position of the system (the chemical potential 41;) and (if z; #0) from the 
addition of charge to a system whose electrical potential is y. 

We now wish to define ji; explicitly in terms of measurable parameters. 
Clearly, jt; is an intensive property of state since it is the ratio of two exten- 
sive properties [equation (1.16)] so that dj; can be expressed as an exact 
differential of all other intensive properties such as temperature (7), pressure 
(P), concentrations (c;, ¢;, . . . ), and electrical potential (y). Thus, 

i 2 3H) 
in er PW Cir Cj, --- oP T,W,Cj.Cj.--- 


OW /7,P,¢j,¢),-.- 9¢;/7,P,,cj,.-- 


On; 


dP 


+2 ee) det * (1.19) 
Fi dc; T,P,W,cj,.-. 

Now, from equations (1.15) and (1.16) 

Bi; = (fe -2 (%) --(¥) 2-5 

ta aT \an;)/ on, \aT oe i. 
and 

7H F) (2) a (32) r) 

al ea Soe aoe See se 7 21 

(e ap \an;) "on, \aP/~ on, (ia 


where 5; is the partial molal entropy of i and 0; is the partial molal volume 
of i. 
Further, from equation (1.17) 


ay 
pedlat b =z.F 
(a T,P, Cj, Cj, --- : 


Now we are close to an explicit statement of equation (1.19), and it only 
remains to define the relations between ji; and c;, cj, .... 

The relation between pi; and c; stems from the brilliant insights made by 
some of the pioneers of physical chemistry around the turn of the nineteenth 
century-Henry, Raoult, DeVries, Pfeffer, Hamburger, Arrhenius, and van’t 
Hoff-which led to the notion that the behavior of solutes in dilute solutions 
resembles the behavior of a perfect gas! Thus, in 1886 van’t Hoff wrote: “The 
pressure which a gas exerts when a given number of molecules are distributed 
in a given volume is equally great as the osmotic pressure, which under the 
same conditions would be produced by most solutes when they are dissolved 
in an arbitrary solvent.” 


(1.22) 
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Thus, in analogy with the perfect gas law, we can write 
P,=n;RT/V =c;RT (1.23) 


where P; may be considered the “partial pressure,” “fugacity,” or “escaping 
tendency” (cf. Lewis and Randall, 1961) of a solute in dilute solution. 


Now, 
(HE) - (Et). (% 
0c; OP, i dc; 


and, from equation (1.23) 


(set) = (2 _ (oP. 
aP; ap) \ap, 
From equation (1.21), (8fj;/P) =v, and from Henry’s law P= P; + L;4;P;, 


we that (0P/dP;) = 1. 
Therefore, (0/1,/AP;) = 0;; and since 0; = (8V/8n;) = (RT/P;) = 1/c;, we may 


Nally write 
ie) _ RT _RTdinc; 
dey <5 os VG sige Cy dc; qa 24) 


Now, what about the term 2; x j(01;/8c;) dc;? For an ideal, infinitely dilute 
tyitem (where the perfect gas law holds) the electrochemical potential of i 
ht hot affected by the presence of other particles (7) so that this term can be 
Heplected; thus, we can combine the results of equations (1.19) to (1.24) and 


obtain 


afi; = -5, dT + 0, dP +z; Fdy +RT dine; (1.25) 


An explicit expression for ji, can be obtained by i i i 
y integrating equation (1.25) 
0 T= 0 to T=7, P=0 to P=P, y=Oto W =p, and c;=1 toc; = ¢;, that 


Hy = wy? - 5,7 +0,P +2, Fp +RT inc, (1.26) 


Where the constant of integration, u,°, is the standard electrochemical poten- 
Hal ond is a function of the solvent alone.” 
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Under isothermal conditions, dT = 0, so that integration of equation (1.25) 
yields 


B= (ue?)r + HP + 27;F¥ + RT Ine; (1.27) 
where (u;°)r is a function of the nature of the solvent and T. 
Under isothermal and isobaric conditions, dT = dP = 0, so that 
Be = (ui°)r,p + Fv + RT Inc; (1.28) 


where (u;°)7,p is a function of the solvent, T and P, 

The practical significance of equations (1.25) to (1.28) is that they provide 
the means for calculating the difference in jij across a membrane separating 
two homogeneous compartments, designated by the superscripts 0 and i, 


which have a common solvent. Thus, in general, 


Ait, = (ii)? - (Hi)! = -5(T? - Ti) + HP? - P!) + 2, FV? - w') 
+RT In (2) 
Cj 


Simplifications of equation (1.29) for specific conditions, for example, 
T° = T', P° = P', are self-evident. 

In developing these expressions for jij, we assumed that the composition 
of the system is sufficiently dilute so that the individual particles (molecules, 
atoms, ions) of i do not interact with each other or with those of other com- 
ponents j,k,/,.... Two problems arise if this ideal condition does not 
pertain. 

First, the perfect gas law analogy ¢; = (P;/RT) cannot be applied. However, 
instead we can write 


(1.29) 


a; = ¥40; = Pi/RT 
where a; is the activity of i and 7; is the activity coefficient. Then, 


pu, = uj? + RT Ina; =u; + RTInc;+ RTI ¥; 


aus) _ RT, pp (210% 
0c; Cj 0c; 


The second problem arises from the fact that in a nonideal multicomponent 
solution, u; is influenced by ¢;, Cx, . - - 80 that Z (0u;/dc;) # 0. These effects 
can be included in the activity coefficient y;. Thus, in general, we can write 


and 


du; =RTdina,=RTd inc +RTdiny; 


| 
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where 


- wv (olny 
diny,= > ( 2c, ) de; (1.30) 


Pquation (1.30) includes all of the effects of c;, cj, ¢,, and so forth on ;. 
let us consider the spontaneous (irreversible) diffusion of dn; moles of i 
from compartment 0 to compartment i. Clearly, dG° = -(i;)° dn; and 
WG! « (ii)! dn;. Since G is an extensive property of state, dG = dG° + dG' = 
(Gi)! ~ (i;)°] dn;. According to equation (1.12), when an irreversible pro- 
ees takes place in a closed system, dW’ =0 and dG <0. It follows that in 
ides for this transfer to occur spontaneously, (ji;)° must be greater than 
,)'. Thus, in a system characterized by a single force (Aji;) and a single 
(dn;), the flow must be directed from a region of higher electrochemical 
potential to one of lower electrochemical potential. Further, since the com- 
{ments must be finite, as the flow continues, Aji; > 0 and dG > 0. When 
ji, * 0, then dG = 0 and the system will have reached equilibrium. 
Pinally, according to equation (1.12), when an irreversible process takes 
Pisce in a closed system at constant T and P 


-dG= Ta;S 


From these considerations, 


-dG = [(ii,)° - (i;)"] dn, = Afi; dn, 
Thus, 
Td,S = Ad; - dn; 


Dividing by dt, we obtain T d,S/dt =J;-Aji;, where J; =.dn,/dt and is the flux 
of in moles per unit time. 
Thus, when an irreversible process takes place in a closed system at constant 
7 wd P, the Gibbs free energy is dissipated at the expense of internal entropy 
=, The rate of formation of internal entropy (or dissipation of free 
py) is proportional to the product of the flow (the rate of displacement 
an extensive property) and the conjugate driving force (a difference in the 
Peaugnte intensive property). 


Conditions of equilibrium 


As discussed previously, when a system is undergoing a reversible 
or passing through a series of equilibrium states, dS = dQ/T so that 
» 0, From equation (1.12) we see that in a closed system at constant tem- 
‘Peete and pressure -dG = Td,S/dt. It follows that in such a system the 
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criterion for equilibrium is dG = 0. If the system consists of two well-stirred 
solutions containing n solutes separated by an inert membrane at constant T 
and P, it will be at equilibrium when all diffusional flows cease, that is, when 
Aji; = 0 for all species i that can permeate the membrane. Thus, from equa- 
tion (1.28) we see that any species 7 will be at equilibrium when 


= (Au;°)r, p + 2 FAy +RTAln C= 


Let us distinguish between the two solutions using the superscripts o and i. 
It follows that for an uncharged solute, Afi; = 0 when 


(1.31) 


where, as noted, the standard state chemical potential y,° is a function of the 

nature of the solvent, T and P. Thus, if both solutions have the same solvent, 
Aji; = 0 when (¢,%lc!) = = 1. On the other hand, if the solvents differ and / is 
more soluble in one than in the other, at equilibrium (c;/c;') # 1. The value 
of (c;°/e#) is frequently referred to as the “partition coefficient” (6;); it will 
be considered further in Section 2.1 with respect to the distribution of an 
uncharged solute across the interfaces between a membrane and the adjacent 
solutions. Clearly, if i is more soluble (more stable) in the solvent in compart- 
ment 7 than in the solvent in compartment 0, (u;°)! <(u;°)° so that, at 
equilibrium, ¢;! > ¢,°. 

Ifi is charged, Ai; = 


c?\__ fw Yury) [ss-v")]} 
(<=) =cxe{ Rr " RT 


Thus, if the two solutions have the same solvent properties with respect to /, 
[(ui°)’ - (u,°)?] = 0 and 


(zr) a [ee 2] 


RT 
sao (EE) 


This relation is referred to as the Nernst equilibrium equation and Ay is fre- 
quently referred to as the Nernst potential..It will be discussed further in a 
variety of contexts, but in essence Aw is the electrical potential difference 
necessary to balance (i.e., equivalent to) the concentration ratio of a charged 
species so that the total driving force (Aji;) is zero. 


0 when 


or 


(1.32) 
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Irreversible (nonequilibrium) thermodynamics 


Classical thermodynamics derives its inviolable authority from the 
it is divorced from any considerations of the structure of matter and 
its attention to systems at equilibrium or undergoing idealized, 
changes (which are, in fact, an infinite sequence of equilibrium 
Under these conditions the properties of state are well defined and 
throughout the system. All real processes are described by a set of 
that simply indicate the direction of change, but there is no con- 
the rate of change and, indeed, the discipline does not formally 
the element of time. Although the amount of useful information 
resulted (and will continue to emerge) from the axiomatic applica- 
iwo postulates boggles the mind, it is clear that, because of its limited 
glassical thermodynamics is not ideally suited for the description of 
ms that are displaced from equilibrium and that are characterized 
of matter and energy within the system as well as between the sys- 
its surroundings. 
or nonequilibrium thermodynamics is a relatively recent out- 
of efforts to extend the principles of classical thermodynamics to 
displaced from equilibrium that are characterized by irreversible 
matter and energy.’ It succeeds in replacing the inequalities of 
thermodynamics with equalities and, by explicitly introducing the 
of time, it deals with rates or flows. However, there is a price to be 
uch as the principles of irreversible thermodynamics are not 
valid and the range of validity must be established (where possible) 
, this new discipline does not command the unquestionable au- 
of its forerunner. 


The linear phenomenological equations 


For many years it has been recognized on empirical grounds that in a 
characterized by a single force and a single flow, the flow is a linear 
of the force over a relatively large range of flows and forces. Thus, 
write 


Jy = LyX; (1.33) 


J, is the flow of an extensive property, X; is the conjugate driving force 
nce in the conjugate intensive property), and Z;; is a proportionality 
1 having units of conductance. Ohm’s law of current flow, Fourier’s 
heat flow, Fick’s law of diffusion, and the Poiseuille’s equation describ- 
me flow are familiar examples of equation (1.33). 

alo known that if a system is characterized by several flows and forces, 
may be interactions (coupling) between nonconjugate flows and forces. 
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Further, for sufficiently slow flows in a system that is not displaced too far 
from equilibrium, the dependence of flows on nonconjugate forces is also 
linear. These empirical observations can be described by a set of linear phe- 
nomenological equations as follows: 


J, = Lyi Xi +L, Xj, + Lip Xy oe AB rp es 
Jj = LuiX;, + Li Xj + Lj Xe es i +LinXn 
Jig = Li Xi + Lig Xj + Lig Xu— tt Lin Xn (1.34) 


Ing = Lig Xt LyjXy + Link Xn t °° + Linn Xp 


or 
i 
Jp=LyiXit Yo LyX; 
i 


In these equations the coefficients Lj;, Lj;,L%%, and so on are referred to as 
straight coefficients inasmuch as they relate the flow with its conjugate driv- 
ing force. The Ljj(;4;) are referred to as cross coefficients or coupling coef- 
ficients inasmuch as they relate flows to nonconjugate driving forces. 

It is also possible to express the forces in terms of linear combinations of 
flows, as follows: 


Xj = Ru Ji + Rif Jj + Riek + °° + RinIn 
X= Ri dit RijJ; t+ Ryde Fite *+RindIn 
Xx = Rie + Rug Jj + Rie + °° + RinIn (1.35) 


Xn = Rnidi + RyjJj +RaxJy + 7° +RanIn 


or 
nm 
X;= Rid + py RijJj 
J 


where the R’s are now generalized resistances. 

It should be noted that the L’s and R’s are, in general, functions of the state 
of the system but are independent of the forces. We will encounter specific 
situations later where the L’s or R’s are functions of concentration but are 
independent of concentration or electrochemical potential differences. 

The choice of equations (1.34) or (1.35) for the description of a given sys- 
tem is a matter of preference or convenience. Some investigators prefer the 
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of resistances since (as will be shown) these can be given physical expres- 
| in terms of frictional interactions. It should be noted, however, that we 
always convert one system of coefficients into the other using the relation 


|Lixl 
Rip = — 


|Z |\s the determinant of the matrix of the coefficients L;, and | L;xl is 
“Minor of the matrix corresponding to the term L;,. Thus, for a system 
‘only two flows (J; and J;) and two forces (X; and X;) 


Li -Lj a 
Ra (it): Ram (its Ree Ct) 


The “dissipation” function 


A fundamental assumption that underlies the entire framework of 
ble thermodynamics is that the total change in the entropy of a sys- 
given by the Gibbs equation (1.5) even though the system is not at 

um and is undergoing an irreversible change. (Recall, that the sub- 

of 7 dS for dQ in deriving the Gibbs equation is strictly valid only 
4 1eversible process.) This assumption has a limited range of validity and 
y be employed when the system “is not too far” from equilibrium. 
final analysis, the justification for this approach in a given instance 
only be established by the experimental validation of the theoretical 

Hons, 
the system is not displaced too far from equilibrium so that the total 
in entropy is given by the Gibbs equation, then using the relation 
8 + dS, it can be shown that 


qS\_ 2 
r (7) J > JX; (1.36) 
i 


J, ant X; are the appropriate thermodynamic conjugates (i.e., Table 
Thus, the inequality (d,S > 0) of classical thermodynamics is replaced 
definable equality. The expression 7(d;S/dt) is sometimes symbolized 
and is referred to as the dissipation function. As discussed before, in a 
system at constant T and P, this function is equal to the rate of de- 
(dissipation) of the Gibbs free energy (~-dG/dt); more generally, it can 
be 1elated to the rate of dissipation of free energy or the rate of de- 
in the ability to perform useful work. 
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An intuitive appreciation of the dissipation function can be obtained by 
considering a simple system in which there is only a single flow and a single 
force. Then, X; = RjjJ; so that 


dS 
T (28) = Rj J," 
Clearly, since T(d;S/dT) is always positive, R;; must also be positive and inde- 
pendent of the direction of the flow. Obviously, if in our example J; is the 
flow of electrical current J across a resistor R driven by a potential difference, 


GS) _ jay 
r (42) PR 


which is immediately recognized as the rate at which power is dissipated (heat 
is generated) by such a system. 

An important corollary follows; namely, in a system in which there is only 
a single force and a single flow, all of the free (potential) energy is dissipated. 
Useful work can only be accomplished by coupling one flow to another. 


Onsager’s reciprocal relations 


In two fundamental papers published in 1931, Onsager interwove 
the linear phenomenologic equations (which were generated by empirical 
observations) and the dissipation function (whose derivation involves a ques- 
tionable assumption) in a brilliant synthetic move that immeasurably strength- 
ened the theoretical credibility of the irreversible thermodynamic approach; 
this achievement and others were recognized by the award of the Nobel Prize 
in chemistry in 1968. 

Onsager demonstrated that if the flows and forces of equations (1.34) are 
chosen appropriately (i.e., the correct thermodynamic conjugates) so that 
equation (1.36) is satisfied then 


Ly = Li for all i andj (1.37) 


One important consequence of these “reciprocal relations” is that they per- 
mit a significant reduction in the number of coefficients that must be experi- 
mentally determined. The number of coefficients needed to completely 
describe a system with n flows and n forces is n*, so that the evaluation of 
these coefficients becomes rather formidable even for relatively simple sys- 
tems (i.e., 2 = 3). When Onsager’s relations apply, the number of coefficients 
is reduced to n(n + 1)/2. Another important consequence is that the experi- 
mental demonstration of reciprocity for a given system supports the conten- 
tion that the system is not too far displaced from equilibrium, so that it can 
be described by a combination of linear relations between flows and forces, 
where the sum of the products of conjugate flows and forces is a measure of 


Some basic principles of thermodynamics 19 


of internal entropy production. Equation (1.37) has been experi- 
sified for a number of systems, including diffusion in temary sys- 
., NaCl—KCI—H,0) and electrokinetic phenomena (see Section 
of these systems were charactérized by rather large forces, suggest- 
‘this approach has a wide range of validity and that the admonition 
jyttem cannot be too far displaced from equilibrium may not be 
etive (cf. Dunlop, 1957; Dunlop and Gosting, 1959; Fujita and 
960; Miller, 1960). 

on additional constraint that emerges from these considerations 
ie positive definite nature of T(d;S/dt). Since the dissipation 
never be negative, the straight coefficients must be positive and 
fficlents must satisfy the condition 


SS buly (1.38) 


led flows, energy conversion, and useful work 


, let us briefly consider a subject that will be treated in greater 
ters 5 and 6. Consider a system in which there are two forces 
that can be described as follows: 


LyX) + LyX; 
why Xi LyX; 
‘the Nows are uncoupled and 


(1.39) 


(4) = JX; + JX; = LyX;? + LX; =R;iJj/ +RyJ;’ >0 


, ull of the free energy stored in the forces X; and X; will even- 
Weeipated. 

ity <0 and if (LyX;|>L,X;, then J; <0. This is an example 
foupling, where the nonconjugate force X; can drive J; in a direc- 
lo the direction of its conjugate driving force and thereby per- 
work. The dissipation function under these conditions is given by 


(4) o J)X; + J,X; = Ly X/ + LyX; + 2L;;X;X; > 0 


_ tinee J) <0 and Ly, <0, T(d,S/dt)' is less than T(d;S/dt) and the 
is the rate at which free energy is converted into useful work rather 
winply dissipated. 

etted is one in which energy from a flow in the direction of its 
driving force is converted into useful work by a (negative) coupling 
thal results in the propulsion of a second flow in the direction 
10 His Conjugate driving force (e.g., the flow of matter against a dif- 


20 Basic principles of membrane transport 


ference in electrochemical potential or volume against a difference in pres- 
sure). The efficiency of energy conversion is given by 


ahs power output _ see) =1-T (48) (1.40) 
power input J;X; dt 


For an ideal reversible process, d;S/dt =0 and n= 1. 


2 Isothermal diffusion 


This chapter deals with the diffusion of a solute across a barrier separating 
two homogeneous solutions in the same solvent and at the same temperature. 
The notation that will be employed is illustrated in Figure 2.1. The super- 
‘eripts o and i refer to the two surrounding solutions and the overbars desig- 
fate intramembrane properties. The bidirectional fluxes of a solute 7 are 
_ designated J;°! and J;/°, and the net flux is given by J, = J,2! - Jf.) For pur- 
poses of simplification we will refer to concentrations (c;) and implicitly 
wsume that y;= 1 (or that 7? = v#); however, it is important to stress that, 
Witlctly speaking, activities should be employed rather than concentrations, 
_Winally, we will assume that each solution is perfectly stirred and ignore 
problems that can arise because of unstirred layers,” 


The Nernst-Planck equation 


Consider the system illustrated in Figure 2.2, where 1 cm? of solu- 
containing n solute particles lies adjacent to 1 cm? of membrane surface. 
tly, the number of particles that will cross the membrane in the x direc- 
per unit time is simply 


Jy = cj; 


fe ¥; is the velocity of each solute particle. If v; is expressed in centi- 
f§ per second and c; in moles per cubic centimeter, then J; has units of 

per centimeter squared, second; when uv; = 1 cm/sec all of the particles 
have crossed from left to right in 1 sec. Now, ; is directly proportional 
The force acting on each particle (f;) and is given by vu; =yju; where the 
Wportionality constant u; is defined as the mobility and is simply the 
ty per unit force.* Thus, we can write the general expression 


Jy = cyuyhj (2.1) 


general, force may be defined as the rate of change of energy with dis- 
(dE/dx) so that the force acting on a mole of matter is simply the gradi- 
Of the free energy per mole or the gradient of the electrochemical poten- 
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COMPARTMENT Ayv-yv'-y° COMPARTMENT 
aS ee eee ; 
(e) ! 
wh) 3 MEMBRANE _ | I 
yl yoy 


re 2.2. A 1 cm? volume of solution adjacent to a membrane 
h a surface area of 1 cm?. 


dient. It should be emphasized at this point that equation (2.2) 

t J; is driven entirely by its conjugate driving force and is not 

other flows or forces. In particular, equation (2.2) assumes that 

convective flow of i resulting from interaction with the flow of 

matter is discussed in detail by Helfferich (1962) and will be con- 
Section 4.4. 


AX i 


Figure 2.1. Notations and conventions: c, y, and J designate con- 


centration, electrical potential, and flux, respectively. The subscript) at constant temperature 
(i) designate the chemical species and the superscripts (0, i) desig- r) = +2,F 2.3) 
nate the compartment or interface. Intramembrane properties are © (uo)r + 0,P+ RT inc; +2, Fy ( 


distinguished from those in the surrounding solutions by overbars. 
The definition Ay = W‘- W° was chosen in order to conform with 
the convention accepted by electrophysiologists; namely, that Aw jy 
the electrical potential of the cell interior or inner solution bathing 
epithelial membranes with respect to the extracellular or outer 
solution. 


because 0; is small for most solutes, the term 0;P makes a very small 
to fi, and can often be neglected, particularly in biological sys- 
that 

J, ~ uj, [RTC In &,/dx) + 2,5 (dwax)] (2.4) 
‘te! us apply equation (2.4) to the description of the diffusion of an 
solute across a membrane whose thickness is Ax. Clearly, when 
tion (2.4) is simply 


4 e ~uye;RT\(d In ¢,/dx) = -RTu,(dc,/dx) (2.5) 


We the flux of i in the x direction at some point in the membrane, u; 

ty of dat that point, ¢, is the concentration of iat that point, and 
js the concentration gradient af that point. When the system is in a 
te, J, is constant and must be the same at all points within the mem- 


tial. Thus, we conclude that 
J, = -cju;(di;/dx) (2.2 


Equation (2.2) is referred to as the Nernst-Planck equation and is the most 
common starting point for the formal description of diffusion. The reason for 
the negative sign is that J; is defined as positive when the flow takes place 
from a region of high electrochemical potential to one of lower electrochemi. 
cal potential, so that a positive flow is driven by a negative electrochemical 
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brane. Assuming that u; is also constant through the membrane, equation 
(2.5) can be integrated across the thickness of the membrane from x =o to 


x =i, that is, 


i i 
FA i] dx = -RTu, J (d@,/dx) dx 
Ce] oO 


yielding 
J; = -RTu; Ac;/Ax 


where AZ; =2;' - ¢°. According to the Einstein relation, Dj = RTu;, where 
D, is the diffusion coefficient, so that 


Jj = ~D; Ac;/Ax = -P; AG; 


where P;, the pormontntity coefficient, is simply D,;/Ax. Now, in general, 
G° #c;° and f;' #c;/, but these can be related by “partition coefficients” 
B;. If 8; is independent of concentration® 


Bi = Ge;° = File; 


and equation (2.7) can be written in terms of the concentrations in the ex. 
ternal solution as follows 
J,= -D; 6; Ac; [Ax = -P; Ac; (2. 8) 


where Ac; =c;'- ¢,° and P; now includes the partition coefficient and \y 
simply 8;P;. Equation (2.8) is the statement of Fick’s first law of diffusion 
for a planar barrier under steady-state conditions. 


2.2, The constant-field equation 


When the solute is charged (z; #0) and there is an electrical poten 
tial difference across the membrane, the situation is far more complex. Under 


these conditions 
i ding, dy 
Jy= - uj; [ar ny) + 2;F ¢2 ) 


or 
(2.9) 


Multiplying both sides of equation (2.9) by exp(z;FW/RT) and rearranging 
the right-hand side of the resulting expression yields 


toon (ir) -0 fa (te 


(2.10) 


(2.6) 


(2.7) 
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steady state and that D, is constant through the membrane, equa- 
ean be integrated across the membrane thickness to give 


for(i lao [ too eel» 


RT RT 


(2.11) 


s y! ie v oO 

, this solution is incomplete because the expression on the left- 
@quation (2.11) remains to be integrated. However, in order to 
this, the dependence of W on x within the membrane (i.e., the 
ntial profile) must be known or assumed. The simplest and 
tly used assumption is that suggested by Goldman (1943), 
v is a linear function of x or that dy/dx = Ay/Ax (“the con- 
‘ystumption”) so that p, = = AW(x/Ax), where x is some point in 
between 0 and Ax. The left-hand side of equation (2.11) then 


! Vv T Ax 2; FA 
5, S08) (4 \| (2 ra 
aa —— | ax = J; exp || —=———_— 
J i ( RT Ax FAY RT 
te solution of equation (2.11) is 
D,2,F AW E exp(2; F¥AW/RT) a ae 
RT Ax exp(z; FAW/RT) - 1 
(2.12) is frequently referred to as the Goldman equation or the 
flux equation and, as originally derived, referred only to intra- 
properties. In an effort to relate J; to the measurable properties of 


solutions, Hodgkin and Katz (1949) assumed that the partition 
across the interfaces is concentration independent so that 


© Be,’ = Bye 


(2.12) 


and 2; 


conditions, as will be discussed in Section 3.1, (W° - ¥°)= 
yo that (U! - 9) =(! - y°) = Ay, and 
) - of FAy/RT) 
Piz FAW [4° - cf exp ] (2.13) 
RT 1 - exp(z; FAW/RT) 


* Dip 1dx. 
(2,13) is often referred to as the Goldman-Hodgkin-Katz equa- 


fiws been widely employed for the description of ion transport 
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across biological membranes. However, it should be stressed that, strictly 
speaking, this equation only applies when the electrical potential profile 
across the membrane is linear and when the partition coefficients at the two 
interfaces are concentration independent. These conditions are most likely 
to be approached by thin membranes in which the density of fixed charges 
is low. On the other hand, this treatment is not generally valid for mem- 
branes (or pathways) that are thick and/or are characterized by a high fixed 
charge density. MacGillivray and Hare (1969) have shown that in order for 
the constant field assumption to be valid, a dimensionless value 


a? = (ko kKRT/¥?X Ax?) 


must be large; Ko is the permittivity of free space, k is the dielectric constant 
of the barrier, and X is the concentration of fixed charge in the barrier. Pur- 
suing this analysis, it can be shown that for a 100 A-thick barrier with a 
dielectric constant of 9, the constant-field assumption will be valid (within 
experimental error) if X <3 X 10M. 

Several other approaches have been taken toward the integration of the 
Nemst-Planck equation. These do not appear to offer any a priori advantages 
over the constant field assumption; the solutions are more complex and often 
unwieldy and, as a result, are encountered infrequently. These approaches are 
discussed briefly in Section 3.4 and are described in detail in Helfferich, 
1962; Lakshminarayanaiah, 1969; MacInnes, 1961; and Sten-Knudsen, 1978. 


2.3. Implications of the Goldman-Hodgkin-Katz (GHK) flux equation 


It is of interest to explore some of the implications of the GHK 
equation with the aid of Figure 2.3, where we have plotted J; versus Ay when 
¢° = 100 mM, P; = 0.01 cm/h, z; = +1, and c/ = 10, 50, or 100 mM; J; is ex- 
pressed in microequivalents per centimeter squared, hour. When J; > 0, the 
net flow is directed from compartment o to compartment i; when J; < 0, the 
net flow is in the opposite direction. 

(a) Clearly, when c;? to the relation between J; and Ay is curvilinear 
and the slope of the curve when J; <0 is less than the slope when J; > 0. 
Thus, even when P; is constant, the membrane offers greater resistance to the 
flow of this monovalent cation (per given change in Ay) when it is directed 
from a lower to a higher concentration (compartment i to compartment 0) 
than when it takes place in the opposite direction. This asymmetric behavior 
is referred to as (“Goldman-type”) rectification. The physical basis of this 
phenomenon will be discussed in Section 7.1, but it should be reemphasized 
at this point that it is not the result of a voltage-dependent change in the 
permeability of the barrier. 

(b) According to equation (2.13), when J;=0, the intersections on the 
abscissa correspond to the values 


Aw = (RT/F) In (¢,?/c/') 


= 100mM 
P, = 0.01 em/br 


zie 


o- 
° 


Figure 2.3. The relations between J; and Ay predicted by the GHK equation (2.13) for a mono- 
valent cation when c,° and P; are maintained constant and c;' is varied. The arrows with a single 
asterisk indicate the reversal potentials, with a double asterisk the points where J; = P; Ac;. 
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This is the point where the electrical driving force balances the chemical 
driving force, and for this reason, Ay is referred to as the Nemst equilibrium 
potential (see Section 1.5) of the system or the reversal potential. 

(c) When Ay = 0, equation (2.13) reduces to J; = P;(¢;° - ¢/') or equation 
(2.8) [recall that x/(1 - e*) approaches 1 as x approaches 0], so that if c;,° 
and c} are known, the intersections on the ordinate provide a direct measure 
of P;. Obviously, when Ay =0, i behaves as an uncharged solute regardless 
of the value of z;. 

(d) When ¢;? = ¢;! = c;, equation (2.13) reduces to 


J; = -c;P;2, FAy/RT (2.14) 


so that the flux of i is a linear function of the transmembrane electrical 
potential difference as shown in Figure 2.3. Since 


1, =2;FJ; (2.15) 


where J; is the current of i expressed in amperes per centimeter squared 
(A/cm?), or coulombs per centimeter squared, second (C/cm? sec), and J; is 
expressed in moles per centimeter squared, second (moles/cm? sec), we can 
write 


I, = -(¢,Pjz;? F?/RT) Ay (2.16) 
In analogy with Ohm’s law, J;R; = Ay, we see that 
R; = RT /c,P,2;?F? : (2.17) 


where R; is the electrical resistance of the barrier to the flow of i expressed in 
ohms centimeter squared (2.-cm?). Thus, R; is inversely related to c; and P;. 
The relation between R; and c; will be discussed in greater detail in Section 
Tea. 


Unidirectional fluxes 


When c;! = 0, the net flux from compartment o to compartment / is, 
by definition, the unidirectional flux J;°', which is most often determined by 
labeling solution o with a radioactive tracer of solute i and following the ap- 
pearance of this tracer in solution i, From equation (2.13) 


ee ° 2z,FAy/RT 

sad ae ee ; exp ae 

When Ay > 0, (z; FAw/R7T)/[1 - exp(z; FAw/RT)] > 1 so that 
oJ! = Pic; (2.19) 


where PY hiad is the unidirectional flux of i when Ay =0; this statement 
is analogous to equation (2.8) for the diffusion of an uncharged solute 
when c;' = 0. 


(2.18) 
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Further, it can be shown that when Ay <25 mV, (z; FAY/RT)/[1 - 
(2; FAW/RT)] = exp(-z; FAW/2RT), so that for small values of Ay 
Jj2! = Pic; exp (-z; FAW/2RT) (2.20) 
milarly, when c;° = 0, from equation (2.13) 


rie pot Ss FAw/RT) exp(z; FAW/RT) 
- 1 - exp(z; FAW/RT) 


ew. (-z; FAy/RT) 
p= Pies [= FAW/RT) - 1 


when Ay <25 mV 
Jf? = Pye} exp (2; F4/2RT) 
ining this last equation with equation (2.20), we obtain 


~ Jf? = Pile,’ exp(-z; FAW/2RT) - cf! exp(z; FAW/2RT)] (2.21) 


is a reasonable approximation when Ay <25 mV and is in error by 
15 percent when Ay = +50 mV (see notes).° 


The independence principle 


An important implication of the GHK flux equation is the so-called 
principle, which asserts that the movement of a single ion at a 
within the membrane is influenced only by the electrochemical gradi- 
that point and is independent of the presence of other ions of the same 
tent species. Three consequences of the GHK equation that illustrate 
tepondence principle are as follows. 
, from equation (2.13), we see that 
A. E ~ cf exp(z; FAY/RT) 
Ii Ley?’ - cf!’ exp(z; FAy/RT) 
oy" and c;"' represent a new set of concentrations that differ from c;° 
" Thus, if we double the concentrations of i on both sides of the mem- 


)olding Ay constant, the net flux will double. 
, for two species 7 and j, equation (2.13) predicts that 


a [pale ~ of exp(z; een | 
Ay 1 Py2/[e;? - of exp (@, FAW/RT)] 
#,% 2),0/° = c/° and cj =e! 


(4,/J,) = (Py/P)) 
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nz, = 0,J, =P; Ac; [equation (2.8)], and when c,? = ¢;' so that 
anes 


0, Py 2; FAy/RT (2.26) 
(2,14)]. 


Finally, equations (2.18), (2.20) and (2.21) imply that the unidirectiong 
fluxes of i are separable and influenced judependently. Thus, J;/° is indepen: 
dent of ¢;° and Bhagat is independent of c;’. Or, stated in another way, there aye 
no transeffects of concentration on a unidirectional flux. 

As we will see below, there are a number of diffusional systems that exhibit 
significant departures from the independence relation and, thus, cannot he 
described using the GHK equation or any integrated form of the Nernl- a : 
Planck equation. Nux-tatio equation 
' of the fact that the Nernst-Planck equation cannot be solved 
toy 4, without additional assumptions, an extremely important 
‘ween the bidirectional (“tracer”) fluxes of i across a membrane 
ean be derived from equation (2.4); this relation was derived 
7) and Teorell (1949, 1953) but was rederived and applied 
9) and generally bears the eponym “the Ussing (flux-ratio) 


2.4. The irreversible thermodynamic approach 


For an entirely phenomenologic description of diffusion of a solute 
across a barrier from one solution to another, it is often more convenient |» 
define the driving force as the electrochemical potential difference (Ajl,) 
across the barrier rather than the negative of the electrochemical potenti! 
gradient (-dji;/@x). This approach assumes that the electrochemical potentis! 
of i just within ae membrane a the interface is equal to that in the adjacen! 
solution, that is, #;° = H;° and #7, ay! = ji;'. In essence, this implies that diffusiw 
through the membrane is sufficiently slow so that the distribution of i acrimy 
the interfaces (a rapid process) is essentially at equilibrium. Although dp, lide 
and Ag; have different dimensions, it can be shown that, given these assumyy 
tions, 


/, a tracer of solute i, is placed in solution o and that its con- 
ty) is ¢;°. Since c;! = 0, the (unidirectional) flux of j from o 
4o equation (2.11) is 


' ' sia | 
"1 exp ze dx = D,z? (2.27) 
J RT “ 


tracer of solute i, is placed in solution i so that its concentra- 
is e!, then since c,? = 0, the (unidirectional) flux of k fromi 


i — — 
a uy [ae | (242%) 
J exp ea dx = D,c,, \exp (2.28) 
if RT ar RT 


#yuations (2.27) and (2.28) assume that D; = Dj = Dx, which is 
# Hotlon that / and k are “perfect tracers” for i. 


T(d;S/at) = 2» Ji(-diij/dx) = D> J;(Ai,) 
i 


It should be noted that the “‘switch”’ in sign is not important since the diffe 
ence in fi; can be chosen so that it is positive and J; (Aji;) > 0 (Katchalyhy 
and Curran, 1965, p. 115). 

Using this approach and recalling that Ay = y'- y? 


J; = Lj (44;) =L,(RTA In Cj ~ 2; FAW) (2.235 


When Ac; = cj? ~ ¢7' is small wation (2.27) by equation (2.28), we obtain 
Jj = Ly (RT(Acj/¢;) - 2; FAY] (2.2%) «[% Jfoo( Gag at ) ies 
where ¢; is an “average” or “mean” concentration given by & 
| hae c} . hot c?) » of j and & through the membrane is slow compared to the 
c= In(c,°/ a 2 Qa) aT ee tracers Aigtetbute across the membrane-solution interfaces, 
i? and py! & by! then, as discussed in Section 1.5 
We now define Lj; = ¢;P;/RT, so that 502) 3 
Jj =P; Ac; - P}C;2; FAW/RT (2 1) *o/" exp [4,2 2) 


Thus, when Ac; is small, J; is the sum of a flow caused by the concenin 
tion difference alone (as if the solute is uncharged) and a flow that is a linew 
function of Ay. 


[- = v ee +) 
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Substituting these relations into equation (2.29), we obtain 


J,2# oO -2; FAy 
Lae ~ [25] [ow (=) 030 


If 7 and k are perfect tracers for the movements of i, then (J;7/c;°) = 
(Fi?/e;?) and (Je!/cy!) = (F;'%e;!) so that 


Je!) _[ ei? 2 Fy? - ¥') 
Pane ee i 
Clearly, when J;?! = J/'°, J; = 0 and 


[2] exp (2. 2a%) Ba) 


RT 


which, as pointed out before, is the Nernst equilibrium equation or the rela- 
tion between the electrical potential difference and the concentration ratio 
when there is no net ion flow. 

The original assumptions underlying the flux-ratio equation (2.31) are: (a) 
the solutions surrounding the membrane are perfectly stirred and the solute i 
(or its tracers j and k) are distributed at thermodynamic equilibrium across 
the two solution-membrane interfaces; (b) the chemical state of i within the 
membrane is identical with that in the surrounding solutions (i-e., there are 
no chemical reactions, associations, complex formation, etc. within the 
membrane); (c) the unidirectional fluxes J;°! and J;!° traverse the membrane 
through pathways that have identical properties [if not, the division of equa- 
tion (2.27) by equation (2.28) to obtain equation (2.29) would not be valid] ; 
and (d) the unidirectional fluxes of the tracers j and/or k are not affected by 
interactions with the flow of the parent isotope i or the flows of other solutes 
or solvent (Ussing, 1949). Recently, more general derivations have been 
described by Hoshiko and Lindley (1967), Kedem and Essig (1965), and 
Dawson (1977). Further, the flux-ratio equation has been shown to be valid 
for the cases of substances formed or consumed during passage across the 
membrane (Ussing, 1952), for membranes that are inhomogeneous in all 
three dimensions (Schwartz, 19715), and under non-steady-state conditions 
(Ussing, 1972). 

The power of the flux-ratio equation resides in the fact that it does not de- 
pend upon a detailed knowledge of the properties of the membrane, and, in 
particular, no assumption need be made regarding the electrical potential 
profile within the membrane. Thus, in obtaining the ratio of fluxes, both D; 
and the integral on the left-hand sides of equations (2.27) and (2.28) cancel 
and the final result can be expressed entirely in terms of readily measurable 
external parameters. For this reason, the flux-ratio equation has provided a 
valuable criterion for determining whether the bidirectional and, thus, the net 
flux of an ion across a barrier can be attributed entirely to electrical and 
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themical potential differences. If the bidirectional fluxes of an ion conform 
} @quation (2.31) over a wide range of concentrations and electrical poten- 
differences, we can confidently conclude that transport is caused by 
uimple diffusion. On the other hand, as wili be discussed in Section 5.2, the 
verse is not true; that is, failure of the movements of an ion to conform 
the flux-ratio equation does not exclude the possibility that its net move- 
nt is energized solely by differences in chemical and/or electrical potentials. 
Finally, it should be noted that equation (2.31) can be derived from equa- 
ons (2.18), (2.20) and (2.21), but it should be stressed that these equations 
on the validity of the constant-field assumption. 


Discontinuous diffusion 


The formal descriptions of diffusional processes derived by integrat- 

the Nernst-Pianck equation, which necessarily assume that (dji;/dx) is a 
tinuous function and which embody the independence principle, do not 
juately describe systems in which the diffusing particle interacts with 
sific sites along its path.” As we will see, such systems can display satura- 
transeffects, and competitive behavior that cannot be described by any 
nent in which P;, the permeability coefficient, is a unique property of 
diffusing species and the barrier and is not influenced by the concentra- 
of the diffusing species or the presence of other permeant ions. 

fipproach that has proved useful in analyzing such systems, based on 
theory (Glasstone et al., 1941), was originally proposed by Eyring, 
ty, and Woodbury (1949) and has been recently elegantly extended by 
nn (1972), Lauger (1973), Hille (1975), Sandblom, Eisenman, and 
1 (1977), and others (cf. Hille and Schwarz, 1978). This approach can 
» mathematically formidable and has not been completely developed 
mplex systems. For these reasons, we will simply try to illustrate the 
lying principles employing relatively simple examples. 
Ye fundamental assumption is that the membrane can be viewed as a series 
‘potential energy barriers, depicted as a series of peaks and valleys, that a 
ile must cross in order to pass from the outer (0) to the inner (7) solu- 
Pvery hindrance to advancement (electrostatic or steric) is represented 
§ peak that denotes the minimum amount of energy the particle must 
1o traverse this obstacle. Conversely, every favorable region for 
nee is represented by a valley from which the particle must escape in 

to advance further. The rate constant for advancement is derived from 
iiitieul mechanical reasoning and, in particular, from the “Boltzmann 
Hibution.” Thus, the rate constant for movement of a particle from valley 
valley &, over the peak jx is given by 


kjx = K(KT/h) exp(- 4G,,/RT) (2.33) 
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Figure 2.4. A simple, double energy-barrier model for diffusion 
across a membrane. 


where AG; is the height of the peak from the valley j, k is the Boltzmann 
constant (k= R/N where N is Avogadro’s number) and h is the Planck con- 
stant. K is referred to as the transmission coefficient and represents the 
probability that a particle that successfully reaches the peak will advance to 
valley k rather than “slide back” into valley /; K is generally assigned a value 
of unity. 

The term multiplied by the transmission coefficient is simply the proba- 
bility that a given particle will possess at least the amount of free energy 
AG;x, necessary to reach the peak. 

We now illustrate the application of this reasoning to the simple two-barrier 
representation of a membrane illustrated in Figure 2.4, where the distance 
between energy peaks is designated as A. If c;° is the concentration of diffus- 
ing particles in the outer solution and c;” is the concentration in the mem- 
brane, then the velocity of flow from the outer solution into the membrane is 


Bee = nee (HIN * Kora (he nN (2.34) 


If the area of the membrane is 1 cm? and A is in centimeters, then c,°A is the 
number of molecules in the lamina with thickness A preceding the first energy 
peak; kK,, has units of seconds” so that v;?” is expressed in molecules per 


centimeter squared second. 
Similarly, we can write for the velocity of exit from the membrane into 
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golution i 


07 = Knici™ A - Kim C4 d (2.35) 


When the diffusion of i from o to i achieves a steady state, cj” is constant 
nd vj°" = v/™ =J;. Also, if the barriers have the same amplitudes, Kom = 
and Kno = k,n; So that solving equations (2.34) and (2.35), we obtain 


Ji = (Kom 4/2) (€7? - ¢7') (2.36) 


_ For an uncharged molecule, equation (2.36) has the same form as equation 
(2.8), so that we may write P; = (kom,/2). Since P; = D;6;/Ax and Ax = Ei. 
follows that D;8; = Ko,A?. 
From equation (2.33), we see that Ko, = (kT/h) exp(-AG,,,/RT). Further, 
© can divide AG,,,, into the difference in free energy of i between the outer 
ase and the valley in the membrane phase, AGg, and the difference be- 
n the valley and the internal peaks, AG,,. Thus, we can write 


D;B; =X? (KT/h) exp(-AGg/RT) exp(-AG,,/RT) (2.37) 


ly, the partition coefficient 6; can be identified with the expression 
(-4G,/RT).° Similarly, the diffusion coefficient can be identified with 
expression )? (kT/h) exp(-AG,/RT), where the exponential term 
ents the effect that energy barriers within the membrane have on the 
lity of the diffusing molecule. 

ral points can be generalized from this simple exercise. First, unlike the 
waches that involve integration of the Nernst-Planck equation where uj 
‘#, and £; are phenomenologic coefficients, the rate theory analysis, in 
fviple, outlines a way in which these intramembrane properties can be 
d in terms of the physical and chemical structural details of the diffu- 
pathway. 

md, as pointed out by Eyring, Lumry, and Woodbury (1949), regard- 
the complexity of the system, the only quantities needed to define P, 
amplitudes of the internal peaks relative to the external solutions; the 
of the valleys between peaks do not appear in the final expression. 
the factors affecting partitioning (binding) and mobility appear as a 
jt [equation (2.37)] and need not be dissociated in order to explicitly 
permeability. 

art, given sufficient information regarding the microscopic details of 
jer, this analysis offers an approach to the description of diffusional 
‘wt the elemental level. 

How consider the diffusion of a charged species through a narrow pore 
tains a binding site X with the restriction that, in order to traverse 
“we, ¢ must bind to X and that X can only be occupied by one i at a 
‘The energy profile of this barrier can be represented as shown in Figure 
| Where Ax is the thickness of the membrane and X designates the posi- 
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A (AX-A) 


Figure 2.5. Energy profile for ionic diffusion through a pore that 
contains a fixed site X capable of binding the ion i. 


tion of X; AG represents the free energy profile caused by both an external 
electrical potential (Ay) and internal barriers. As discussed by Hille (1975), 
this system can be described by the equations 


Bee 


of? + [X c. [c;X] 
kxo* 
and (2.38) 
Kyi" ; 
fe x) c+ [4] 
Kix 


where the brackets designate the concentrations of the free and complexed 


forms of X. 

Equations (2.38) are strictly analogous to the equations of chemical kinetics 
describing the reversible reaction between a substrate and an enzyme leading 
to the formation of a product. The rate constants are given by the following 


relations: 
Kox * = Kox exp [-z F(1 - a) Wy /RT] 
kyo * = kyo exp [z Fay ™/RT] 


kgi* = ky, exp [-2F(1 - 8) W*/RT] (2.39) 
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Kix * = Kix exp (2 FBW*/RT) 


Where y°~ is the electrical potential at X with respect to the outer solution 
that Ay = y'- y° = y* + y*!. The k’s (without asterisks) are indepen- 
nt of the external field but are dependent on voltage-independent internal 
nergy barriers as described by equation (2.33). The values of a and § are 
pendent upon the electrical potential profile between o and A and between 
and 7; if these profiles are symmetrical about the midpoints the rate con- 
Jants will be equally and oppositely affected by the applied potential and 
6 =0.5 (Hall, et al., 1973; Stein, 1977; Glasstone et al., 1941). 

‘From equation (2.38) the net flux of i from compartment o to compart- 
nt i under steady-state conditions is given by 


Ip = Kox*c;°[X] ~ kxo* [eX] = kxi* er X] - kix*e'[X] 
her, the total concentration of binding sites per unit area [X], is equal to 


ie concentrations of those sites occupied by i, [c;X], and those which 
ain unoccupied [X], so that 


[X],= [X] + [eX] 
dlving these equations for J;, we obtain 
Kon ki*¢,? = Kae Fey CF 
=[X 
4 (XI. Ge +kyj* + kix*e;! + ko,*c;° 
m equations (2.39) and (2.40) it can be readily shown that J; = 0 when 


ce = es) # Gate) exp (2 = 
cf! fay a RT 
(uation (2.41) conforms to the Nernst equilibrium expression (2.32) only 
WON (Kix Kxo/Koxkxi) = 1, which is one of the constraints on the rate con- 
Win a strictly diffusional system. 


om equation (2 aw, it is clear that the net flux from compartment o to 
ment i when c/ = 


(2.40) 


(2.41) 


= 0 is given by 


\g Kxi*c;°) 
td ae kyi* - | 


Nearly, equation (2.42) has the form of the familiar Michaelis-Menten 
tion of enzyme kinetics. J; saturates with increasing ¢,° (at constant 
the maximum influx under these conditions is [X],k;" and a half- 
influx is achieved when ¢,° = (kyo* + kxj*)/Kox*. 

ersely, when c,° = 0, the net flux from compartment i to compartment 
on by 


J:= (Xe (2.42) 


Kix "kx" of! 


+ Kit + ke (2.43) 


X 
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so that the maximum outflux is [X],xo* and the half-maximum outflux is 
achieved when c/ = (Kxo* + kxi* /kix*- 

Equations (2.42) and (2.43) disclose an extremely interesting property of 
this system; namely, the kinetic properties that describe the relation between 
J; and c;° when c} = 0 differ from those that describe the relation between J; 
and c when c;° = 0. Even when Ay = 0, the membrane will behave sym- 
metrically only when k,o* = ki" (ie., the maximum fluxes in either direc- 
tion are equal) and kj,* = kox* (i¢., the concentrations needed to achieve 
half-maximum fluxes are equal). Otherwise, the membrane will display recti- 
fication; that is, for a given value of c; the flow will differ when 7 is in com- 
partment o alone from that observed when / is in compartment i alone. 

Another interesting feature of this system emerges when we consider the 
bidirectional fluxes determined with tracers. It can be readily shown that if a 
perfect tracer for i, k, is placed in compartment 0, the tracer flux is given by 


Jot = [X] ~ ee 
7m * €,0 * 
* Ko” t Kat? t Kony? + Kix "Ci 


If k is a perfect tracer for i, it follows that (Siler?) = (J;%/c;°) so that 


= * k Ag cf 


2.44) 
[xX] Recs’ tfica * Rae CF kix*e? ( 


P Fad = 
Similarly, the unidirectional flux from compartment i to compartment o is 

given by 
Tene Gis 2, 


lo —__—__ ix x0 1 (2.45) 
a 1. bn 2] ks kage .* kine? 


Thus, the unidirectional fluxes exhibit transeffects; that is, the unidirectional 
flux from o to i is affected by c/ and the unidirectional flux from i to o is 


affected by c;°. ' 
Finally, from equations (2.44) and (2.45) it is clear that the flux ratio is 


given by 
se! r (2) cae 
J, ° ¢;! kix*kxo® 


which from equation (2.41) reduces to 


je he ce ee (4 =) 
gy? “ cf P\ RT 
or the Ussing flux-ratio equation. ; a 
The behavior of this very simple system obviously deviates significantly 


from those which conform to the independence principle, inasmuch as there 
is no unique permeability coefficient P; that relates the net flux or unidire¢ 
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onal fluxes to differences in concentration and/or electrical potential. In 
~ontrast with systems described by the GHK equation (2.13), where P; is in- 
dependent of c; and Ay, in this system P; is explicitly dependent upon these 
eters. In addition, these systems display rectification (even when 
W=0) and transeffects. Finally, although this system conforms to the 
Issing flux-ratio equation, as will be discussed in Section 5.2, pores involving 
ultiple binding sites that permit only “‘single-file diffusion” do not conform 
this relation. 
Finally, as discussed by Hille (1975), and others (Heckman, 1972), the 
jple system illustrated in Figure 2.5 and described by equations (2.38) can 
hibit competitive interactions that will influence the flows of i when 
her species with an affinity for X is present in the outer and/or inner 
ih ion; this should be intuitively obvious. 
$s approach has been successfully applied to the description of ion trans- 
ft across artificial lipid membranes doped with “pore-forming” antibiotics. 
pores appear to possess discrete sites with which a traversing ion inter- 
, and these sites appear to be separated by distances that are large com- 
to the size of the ion. Ionic currents through such systems exhibit 
ion and are subject to competitive inhibition and, thus, do not con- 
to the independence principle. However, it should be noted that if we 
¢ that the membrane is comprised of an infinite number of energy 
iors (i.c., approaching a continuum), the expressions resulting from the 
theory analysis resemble those derivable from the Nernst-Planck equa- 
and the GHK flux equation can be derived if, in addition, we assume 
if all of the barriers are of equal height (Hille, 1975). 
»ther, but similar, approach to systems in which there are discontinui- 
in the chemical potential of the diffusing species i has been described by 
i and Kedem (1966), who view the membrane as a lattice made up of inde- 
{ units each of which can exist in a finite number of discrete states. 
“subsequent paper, Essig, Kedem, and Hill (1966) examine net flows and 
stional (tracer) flows in such lattice models and explore the conditions 
which the flux-ratio equation (2.31) is obeyed. 


A few words about permeability coefficients 


One serious deficiency of all approaches toward a formal description 
te transport across membranes is that there is no generally applicable 
| 10 calculate P; from a knowledge of the properties of i and those of the 
ivy The P; that emerges from the various integrated forms of the Nernst- 
Hk equation is strictly a secondarily derived, phenomenologic (“black 
Coefficient that relates the observed flux to the driving force. The 
‘theory approach outlines, in principle, the energetic factors that deter- 
#,, but the application of this approach to complex membranes presents 
lable if not insurmountable practical problems. 
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Renkin equation can provide an internally consistent body of 
“relative permeabilities,” “equivalent pore radii,” and so on, 
may not be very realistic. 

there is no generally accepted approach to the calculation of the 
wvefficients of charged species from the properties of the ion 
1. Clearly, a formalism that would permit an independent cal- 
from first principles would provide an independent check of 
the flux equations employed to describe diffusional flows. 


Limited success in deriving values of P; for the very simple case of diffus 
of small uncharged molecules through water-filled, right-cylindrical py 
was first achieved by Pappenheimer, Renkin, and their collaborators (Pap 
heimer, 1953). In analogy with equations (2.7) and (2.8), they wrote 


J; =Dj,As Ac,;/Ax 


where D,, is the diffusion coefficient of the solute in an aqueous solut\) 
(rather than in the membrane phase) and A, lumps all of the membrane py 
erties that influence entry into and movement through the pore. Il) 
P; = Dj, As/Ax. Renkin (1954) assumed that the total restriction on diffu 
through the pore, compared to unrestricted (free) aqueous diffusion, cu) 
divided into (a) steric hindrance at the orifice of the pore, which influe 
the probability that a randomly moving particle will enter the pore (jar 
tioning) and (b) hindrances to movement within the pore arising from I) 
tional interactions with the wall (i.e., the difference between D;, and /), 
Pappenheimer, Renkin and Borrero (1951) suggested that the likelihood | 
a particle will enter a pore (i.e., the probability that the particle will 
collide with the rim of the pore) is given by [1 - (r/rp)1’, where 7; is f 
molecular radius of the solute and r, is the radius of the pore. The effec! 
frictional interactions within the pore was described by the Faxen equal! 
which was derived for the case of spheres moving through a fluid cap} 
assuming that the spheres are much larger than the fluid molecules (ie, | 
fluid can be treated as a continuum). Using these assumptions, Renkin (1')) 
showed that 


e 2 3 5 
! - 2.109 (2) + 2.09 ie) - 0.95 a ( 
Pp Pp P 


where A,/A is the fraction of the total membrane area A occupied by jy 
orifices Ap. Thus, if A, Ax, Ap, and r, are known for a given membrwiye, 
can be calculated from the knowledge of 7; and D;, (which can be deter) 
experimentally). 

The Renkin equation (2.46) has been verified in experiments using \W 
branes that have uniform, well-defined thicknesses and that are penetrate! 
pores whose radii are known (Beck and Schultz, 1970; House, 1974), I 
ever, it should be stressed that for the case of biological membranes: (} 
and A, cannot be determined directly; (b) the “pores” (if any) are aly 
certainly tortuous channels that cannot be characterized by a single 1) 
and (c) the pore dimensions are such that the molecules that are likely 
enter are not very much larger than water, so that the assumptions wn 
lying the use of the Faxen equation are unrealistic. Thus, while the ay) 


3 Diffusion potentials 


In Chapter 2 we considered the effect of electrical potential differences = 
the diffusional flow of ions across 4 barrier without reference to the origin a) 

these potentials. In the following sections we will examine how the — 
of ions gives rise to electrical potential differences or diffusion potentials, 
starting with some simple, specific examples and concluding with more gen- 


eral treatments. 


at: Diffusion potential between two solutions of the same salt 
-stirred 


Assume that a homogeneous membrane separates two well , 
an 


solutions (o and i) of a salt that is completely dissociated into a cation Cy 
an anion c_, where Z, =~2Z-. According to the Nernst-Planck equation, 


= dy 

Pera [er (2 ne) rie (2) 
and 2. 
= dv 

s.2-eu_[ar (U8) -2.9 (2) 


must be preserved,’ it follows 
=J;. Equating (3.1) and (3,2) 


(3.1) 


(3.2) 


Since bulk electroneutrality of each solution 
that c,? =c.° =¢j", c+ =C- = cj}, and J, =J- 
and rearranging the terms yields 

dy u.-u- (#7) (2n2) 

ax u,tu_} \Z.F dx 
and integrating across the thickness of the membrane, we obtain 


a 
z? 


of a Oy Sibe at 

bn. gets Be BS ee 
ed ia ape ec 
Assuming that 6; = Gi lef = Gi°/ci? , we obtain? 


1 gp a (22) in(£5) 
Avaal-# Bans ZF ec? 


3.3) 


42 


(3.4) 
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nz, # ~z_, it can be readily shown that 
u, - UL RT cf 
oe Oe (— 
Lu, tza) \ F e;”, 
artz, 1971b). 
stituting equation (3.3) into equation (3.1) and combining terms yields 


= - PARE) 
uy, tu. 


de; 

dx 

J,, = J;, equation (3.5) can be rewritten 

J; = -Ds(de;/dx) 

D,, the diffusion coefficient of the salt, is simply 


(2RTu,u_){(u, + u_-)=RTuz 


important principle to be gleaned from this simple example is central 
Understanding of the origin of diffusion potentials in more complex 
is. Namely, in the absence of an extemally applied current, electro- 
can only be preserved by the equivalent flow of anions and cations 
the membrane. This so-called zero current condition is satisfied in this 
system by the constraint that J, = J_ = J; so that 


[= 1, +1_=2,FJ,+2-FJ_=0 (3.6) 


u_, equation (3.6) can only be satisfied if an electrical potential dif- 
is generated whose magnitude is proportional to the difference in 
and whose orientation is such as to retard the movement of the ion 
greater mobility and accelerate the movement of the ion with the 
mobility. Thus, if the cation has a greater mobility than the anion, the 
‘solution will be electrically positive with respect to the more concen- 
solution and vice versa. The result is that the two oppositely charged 
move with the same velocity, and the diffusion of a dissociated salt 
characterized by a single diffusion coefficient that is given by the 
ie mean of the two mobilities. 

in another way, if u,=u_, then when the same force acts on both 
v, @ uv. and J, =J_; thus, the only force needed is the concentration 
(de, /dx) and, according to equation (3.4), Ay = 0. However, if 
, another force must be present that exerts opposite effects on the 

4 in order that J, =J_; this force is the diffusion potential Ay. 

yw examine the limiting condition where the membrane is ideally 
to one, but not both, of the charged species. According to 
{3,4), ifu,=0, 


» (RF sr) 
Av (£7) in 


(3.5) 


(3.7) 
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If u.=0, then Aw has the same magnitude but the opposite orientation. 
However, if both ions cannot cross the membrane, then neither is permitted 
to cross, otherwise electroneutrality would be violated. Thus, if either u, = 0 
oru_=0, J, =J_=0 and D, =0. Therefore, in this limiting condition, there 
is no movement and the system is in a state of equilibrium. The magnitude 
and orientation of Ay precisely oppose the chemical potential difference of 
the permeant species so that there is no net driving force acting on that 
species and, hence, no flow. 
More explicitly, if u, = 0, we can rearrange equation (3.7) as follows: 


z_ Faw = -RT in(c!/c-°) 


or Afi. = RT In (c,'/c,°) + z_ ¥ Ap = 0. Since Afi_= 0, the anion is at equi- 
librium; Ay is referred to as the Nernst equilibrium potential of the species 
and may be viewed as the electrical equivalent of the concentration difference 
(see Section 1.5). 

We can now appreciate why interfacial potentials cancel [see derivation of 
equations (2.13) and (3.14)] if (a) translocation of an ion through the barrier 
is slow compared to the rate at which the ion can partition across the inter- 
face so that the interfacial distribution can be considered to be an equilibrium 
distribution and (b) the partition coefficient 8; is concentration independent. 
Clearly, under these conditions 


(W° - p°) = (RT/z;F) In 8; = (W! - ') 
so that 
(wi - py?) = (Wi - $9) = Ay 


3.2. The Goldman-Hodgkin-Katz equation 


The problem becomes vastly more complex when we consider a sys- 
tem in which the membrane is surrounded by two solutions containing 
various anions and cations, not simply a single dissociated salt. The simplest 
and most frequently employed approach to such a system was originally 
derived by Goldman (1943) and later rederived and extended by Hodgkin and 
Katz (1949) for the restricted conditions where (a) anions and cations are 
univalent; (b) the permeability coefficient of each ion is constant, concentra. 
tion independent and is given by u;8;R7/Ax; and (c) the electrical field 
across the membrane (dy/dx) is constant. Under these conditions the flow of 
each cation is given by [equation (2.13)] 


_P, Fay [= -e,! op SAPD) (3.8) 


i? ~~ os 1 - exp( FAW/RT) 
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and that of each anion by 


Peay c/° - c! exp(-F oa 
RT 1 - exp(- F Ay/RT) 


Under “zero current conditions” (see Section 3.1) 


I= 9(E4- Ls) =0 (3.10) 


7 (3.9) 


fered, is the sum of the flows of all cations and Z, J_ is the sum of the 
of all anions (recall that a flow is defined as positive when it is directed 

fom o to i). 
Combining equations (3.8), (3.9), and (3.1 0) and solving for AW, we obtain 
Goldman-Hodgkin-Katz (GHK) equation 
Pyc.° + > Pc! 

Bere OT a 
F ¥ Fe, vy foe G.11) 
c a 


many biological systems where the predominant, permeant cations and 
sare Na, K, and Cl, equation (3.11) is simply 


Ay = RT in( Sesh + PxcK? “Facet 


Pyatna’ + Px! + Poiccy” ‘idea 


arly, if an ion / is distributed at equilibrium across the membrane so that 
= 0, then J; =0 and this ion will not participate in the generation of Ay. 
for example, if Afic, = 0, equation (3.12) reduces to 


RT Cra’ tacK? 
Aw = (——} tn { N82 SS K_ 
/ = (ees ae 
re @ = Py /Pya. 


Generalizations of the GHK equation 


As noted, equation (3.11) was originally derived by assuming that 
anions and cations are monovalent,? have constant permeability coeffi- 
ts, and that the electrical potential profile across the membrane is linear 
», 8 Constant electrical field). Subsequently, a number of generalizations 
been introduced that relax some of these restrictions and extend the 
of validity of expressions that have the form of the original GHK equa- 
(3,11), Thus: 

) If the membrane is permeable only to ions having the same sign and 
nee, then (Conti and Eisenman, 1965, 1966) 
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P;¢;° 
RT 2X “ 


i. (3.14) 
ana > Pei! 
i 


Ay = 


Equation (3.14) can be derived directly from the Nernst-Planck equation 
(2.4) without any assumptions regarding the electrical potential profile within 
the membrane or interfacial potentials. ' 

(b) Equation (3.11) is generally valid when the total concentration of uni- 
valent ions (or ions having the same absolute valence) is the same on both 
sides of the membrane; under this condition, it can be shown that the elec- 
trical field within the membrane is constant (Teorell, 1953; Finkelstein and 
Mauro, 1963; Sandblom and Eisenman, 1967). Mg 

(c) Equation (3.11) is valid if the individual permeability coefficients are 
not constant providing that the ratios of permeability coefficients, given by 
the ratio of the products of individual mobilities and partition coefficients, 
are constant (Sandblom and Eisenman, 1967); that is, 


(P;/P;) = (u;8;/u;8;) = constant for all i andj 


(d) The “constant field” assumption is too restrictive. The only constraint 
on the electrical potential profile within the membrane required to derive 
equation (3.11) is that 


i 
J sinh (F W/RT) dx = 0 


This constraint is consistent with the family of electrical potential profiles 
that are antisymmetric (“odd symmetry”) about the midpoint of the mem. 
brane (Figure 3.1) (Mullins, 1961; Barr, 1965; Schwartz, 1971a, b). 

(e) Finally, equation (3.11) has been extended to include conditions when 
interfacial electrical potential differences do not cancel (e.g., the partition 
coefficients are dependent on concentration) and when the electrical poten: 
tial profiles within the membrane are nonlinear; however, these treatments, 
though mathematically sound, are difficult to apply in practice ee Schwartz, 
19714, b; Jacquez, 1971; Barr, 1965; Jacquez and Schultz, 1974). 


3.3. A general expression for a diffusion potential 


We now derive an expression for a diffusion potential, building on 
some of the principles outlined previously, which is generally valid under 
isobaric, isothermal conditions and does not depend upon assumptions re 
garding membrane composition, electrical potential profiles, or ionic mobae 
ties. More rigorous derivations may be found elsewhere (Staverman, 1952) 
Scatchard, 1953; Kirkwood, 1954). 
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MEMBRANE 
MPARTMENT I 
ce) 


COMPARTMENT 
i 


x=0 X= AX/2 X= AX 


Figure 3.1. Examples of the family of antisymmetric electrical 
potential profiles that are consistent with the GHK equation (3.11). 
The constant-field profile is given by line c. The other extreme, 
designated by a, is when the electrical potential difference is the 
result of two interfacial potentials in the same direction; in this case, 
there is no electrical field acting on the ions during passage through 
the membrane and if (as shown) the interfacial potential differences 
are of equal magnitude, J; is given by equation (2.21) (Mullins, 
1961). The curve marked b is an intermediate, antisymmetric profile 
between the extremes a and c. 


ler the diffusion of m charged species across a membrane surrounded 
homogeneous solutions in the same solvent under steady-state condi- 


so that the flows and forces are time independent. According to equa- 
1,36), 


T(d,S/dt) = 5°J;,X; 
i 


, #8 discussed by MacInnes (1961), the barrier may be viewed as con- 
of an indefinite number of layers, or lamina, so that the differences in 
foncentrations of i across each layer can be made infinitesimally small 
ihe transfer of i across each layer can be considered a reversible process. 
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Thus, for each layer we can write 
n 
T(a;S/dt) = > J,X;=0 
i 


or from equation (2.4) 


fafa (288) +a0(@) a0 
or 
3 Js [er(22)| + > zJi (2) =0 (3.16) 


The total current across this layer is given by J = L} 2; FJ;, so that equation 
(3.16) reduces to 


(B+ (2) a(t am 


We now define the transference number t; as° 


t= Gill = Gi FIi/Dacj=0 
Thus, t; is the fraction of the total current across the layer attributable to the 
flow of i when Ac, = 0; clearly, 2; t;= 1. Substituting this expression into 
equation (3.17), we obtain 


{a BE. (ui) (ama) 
ax F 4 2j dx 
Integrating this expression, for a single lamina, across the entire thickness of 
the membrane from solution o to solution i, we obtain 


Ay -- - [> (4) aina| --if 3 (<) FERC) 
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Figure 3.2. System employed for the measurement of transference 
numbers across a membrane. The total current flowing through the 
external circuit is the sum of oppositely directed cationic currents 
and anionic currents. Electroneutrality of the two solutions is main- 
tained by “electrode reactions.” For an excellent description of the 
determination of transference numbers and electrode reactions, see 
MaclInnes (1961, pp. 59-96). ’ 


fer of matter across any layer can be treated asa reversible process 
expression for Ay is generally valid. 
er, in order to integrate the right-hand side of equation (3.18), the 
between ¢;, the transference number of the membrane as a whole 
# must be known. In order to illustrate the problem, let us consider the 
shown in Figure 3.2, where a membrane separates two identical solu- 
eontaining an arbitrary mixture of anions and cations; thus, c,° = c;! for 
If we impose an electrical driving force Ay across the cenhiae by 
of the battery, the total current across the membrane will be made up 
eationic current /,, directed from the anode to the cathode (from com- 
nto to compartment i in the figure), and an anionic current J_, flow- 


i i 


It should be emphasized that the integration includes the boundary cond}. 
tions so that equation (3.18) refers to the properties of the external solutions 
(the initial and final states of the system). A rigorous justification for this 
procedure is given by Kirkwood (1954); in essence it rests on the fact that the 
electrical potential difference across any layer is a property of state (i-e., inde 
pendent of path) so that the overall electrical potential difference (AW) 's 
simply the integral (summation) of these elemental electrical potential diffes 
ences. It also follows that although equation (3.18) was derived assuming thu! 
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ing in the opposite direction. Thus, from equation (2.4), we can write for 
each cation and each anion 


J, = -2,0,U, Fay and J_= ~z_¢_u_ FAW 


(recall that Ay = y' - y?, so that in Figure 3.2, Ay is negative, J, is positive, 
and J_ is negative). Now, we relate the concentrations of i in the membrane 
to those in the outer solutions by means of the partition coefficients 6; so 
that 


J, = -24B,C,U, Fay and J_= -z-$-c_u- FAw 


It follows that the otal cationic current is given by J, : -2, 2478. C,us FAV, 
and the total anionic current is J. = - D4 z-"B-c-u-F* Ay. 
The total current is given by J = J, + /_; thus 


1=-), 2,"B.c.u,F’ dp - > z2p_c_u_F? Ay 
c a 


or 


T= -> 27Bjc;ujF? Ay 
i 


Since for any ion i, J; = -z;7B;¢; uF" Ay, it follows that 


_ [2,7 Bie) (3.19) 


nA (i) . [eee “] 


Clearly, 2; t;= 1. Thus, in general, ¢; is a complex function of concentrations 
(activities), partition coefficients, and mobilities. 

Let us now consider the simple case where the membrane separates two 
solutions containing a single salt at concentrations c;° and cj. Let us further 
assume (as before) that z, = -z_ and that the membrane is “neutral’’ so that 
B, = 8. = 6; and ; is independent of concentration. From equation (3.19) it 
follows that 


inant ieee (3.20) 
bi” ig Brie 009.4 =F ieee, 


Thus, if the mobilities are independent of concentration, ¢, and ¢_ are con- 
stants and equation (3.18) can be written as follows: 


RT P (22) f 
=-\|-——]t, ainc, + |——/t- dinc- 
Ay (=) t J ne, + | x . 


Since the preservation of electroneutrality requires that c, =¢c_=¢; through: 
out the system, it follows that 


i 
av= -(£2) «1 f din ¢ (3.21) 
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since t, + ¢t.= 1, 


i 
Ay = (4 f“i- 2) f ding, (3.22) 


Integrating equation (3.21) across the total membrane, we obtain 


i 


{ substituting the relations between ¢; and u, given by equations (3.20) 
fo equation (3.23) yields 


~ BT (uate) 4 (at) 

al oh (i: + a in( 

ch is identical with equation (3.4). 

In order to apply equation (3.18) to more complex systems, the relations 

ween t; and c; must be determined experimentally. This empirical relation 

‘then substituted into equation (3.18) and the final expression is integrated. 

iis approach has been employed by several groups and excellent agree- 

nt has been found between the predicted and observed values of Ay 
hminarayanaiah, 1969; Gunn and Curran, 1971; MacInnes, 1961). 


The Henderson, Planck, and Schlég) equations 


As discussed above, equation (3.18) is derived using the Nernst- 
ck equation (2.4), which can be solved explicitly under certain restricted 
\ditions. However, in general, the integration of equation (3.18) or equa- 
(2.4) requires some assumption regarding the electrical potential or ionic 
iivity profiles within the membrane. As we have seen, one such assumption 
hat the electrical field within the membrane is constant (i.¢., p is a linear 
on of x), which leads to the GHK equation (3.11). 
nderson (1907, 1908) integrated equation (3.18), assuming constant con- 
ation gradients within the membrane; or, more exactly, assuming that 
fegion separating the two surrounding solutions can be viewed as a con- 
ious series of solutions produced by the gradual admixture of the two 
nal solutions. Thus, at any point x within the membrane 


Bi* = G1? + (Gi - 1°) x/Ax 
from the definition of t;, equation (3.19), 
Gj*ujz? Ax 
| a n n 
(Ax ~ x) >, GP u,2,? +x > Gay2,? 
i t 
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If u; is independent of concentration over the range ¢;° to cj‘ for all ions, 
equation (3.18) can be integrated to yield the Henderson diffusion equation 


all ions are monovalent and C,° = C; 


gui - Ue? a 
; win gpo 
Dd ule - &°) aa 
p< : In (3.24) 
n 
d ui(Gi - 1°) 21 _ vis ue 
i [= V? re yi 
If all of the ions are univalent so that z, = 1 = -z_, equation (3.24) reduce} 
to 
— RT [(u°-v’)-('- a be + oa a, per 
2 = | tp (| 3.25 Cy uy t Ca Ue 
ee” erp "Vole Ga ap. RT 2 Fue + De 


where F 4 pas c,/u,+ D @-?u- 

¢e a 
O xe = 0 = =o 

Ee Fettcs FSS aw identical with the GHK equation (3.11). 

, for the case of a single univalent-univalent salt, where c;° #c;', 

(3.27) reduces to equation (3.4). 

in (3.27) has been extended to:the case where all cations have one 

‘and ali anions have another (Johnson, 1904), and a general formula- 

atbitrary mixture of ions has been derived by Pleijel (1910). 

¥, equation (3.27) is rather unwieldy and difficult to apply except 

the most simple circumstances; the extensions of this expression are 

cumbersome. Planck suggested a graphical method for solving 

(3.27) for — when the compositions of the two bathing solutions and 

mobilities are known. Thus, we let 


(0 ang | (EC%)- Ine (ff) - 
y'-iv? In (C'/C°) +n} \C#- eC? 


(3.20) 
W=yalu,, Vi=y fu 
a 


and Dc is the summation over all cations and Ya is the summation over il] 
anions. For the case of a dissociated univalent-univalent salt, equation (3.25 
reduces to equation (3.4). 

Planck (1890a, b) integrated equation (2.4) for the case where z,=-z_" | 
assuming that diffusion takes place 2cross a “constrained boundary”’ separ! 
ing two well-stirred solutions and that microscopic electroneutrality is clowly 
approximated at every point within the barrier; that is, there is no signifle 
space charge density within the membrane or 22;¢, = 0. [It should be not 
that this assumption differs only slightly from the constant field assumptiw) 
Thus, according to the Poisson equation 


zt) (Ge) (e) Ee" 
aE = | ee z 
(S dx € x i 
where E is the electric field and ¢ is the permittivity of the barrier. Clow 
when £ is constant, 2; z;¢;=0; this is the condition of strict microm, 
electroneutrality throughout the membrane.] The result of Planck’s intoy 
tion® is the transcendental equation 

gui-u? _fin(C'/C)- ne) (S an 

vi-tv? | in(C'/C°)+inéJ \ C'- EC? 
where U°, V°, U', V! are defined as in (3.26); C° and C! are the total (i 
concentrations; and & = exp( FAW/RT). 


and b are evaluated for different values of £ and plotted against the 
‘Yilues of ; the point where a and b intersect defines the value of & 
with equation (3.27). 

, Sehldgl (1954, 1964) has derived a very general solution of the 
k equation for arbitrary mixtures of electrolytes without as- 
fegarding the electrical potential or concentration profiles within 
fine. The derivation is very complex and the final expressions are 
difficult to apply in practice. The reader is referred to the original 
the brief outline by Helfferich (1962; pp. 386-8) for further 
it to say that when both solutions contain equal total concen- 
of ions having the same valence, the predicted electrical potential 
Wear and the final expressions resemble the constant field equation. 


(i) 
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3.5, The Donnan equilibrium and the Teorell-Meyer-Sievers 
(TMS) model 


Up to this point we have limited our considerations to systems 
where the partition or distribution coefficient between the membrane inter. 
face and the adjacent solution 8; is assumed to be independent of concentra, 
tion so that B,;° = 8; for all i. However, in general, this assumption is not 
strictly valid even for so-called neutral membranes and is grossly incorrect for 
the case of charged membranes that contain a high density of fixed (struc. 
tural) charged (dissociated) groups.” In this section we will examine the be. 
havior of charged membranes, where B:° #8;', and develop a general ap: 
proach toward the analysis of transmembrane diffusion potentials. 

Let us start by considering a system where compartment o contains a solu 
tion of a salt that dissociates into two univalent ions at concentrations es 
and c_°, and compartment i contains a solution of the same ions, c,! and c! 
as well as another charged anion X_' with valence Z,,. We assume that the 
intervening membrane is permeable to c, and c_ but impermeable to X ! 
When the system achieves equilibrium, ji,° = ji,' and f° = 4, so that from 
equation (1.32) 


Ay = y' - y° = -(RT/F) In(c,!/c,2) = (RT/F) In(c/e_°) (3.28) 


We now define r = (c,"/c,°) = (c_°/e-'), so that Ay = -(RT/F) In r. Now, 
since both external solutions must be electroneutral, 


=e =e," and c.f =c!+|z,(xX! (3.29) 


Substituting the definition of 7 and rearranging equations (3.29) readily yields 


Cy 


Po? - rlz,X-!|-¢,° =0 


and, solving this quadratic equation, we obtain 


x x 2 1/2 
t= 7 + (2) + 1 (3.30) 


(Note: The positive root is chosen because the square-root term must always 
be larger than the first term and r cannot be negative.) 

This equilibrium system is referred to as the Donnan (or Gibbs-Donnan) 
equilibrium and the value r is often called the Donnan ratio. 

Now, let us consider a membrane that contains a high density of fixed 
negative charges that separates two solutions of a dissociated salt cy, whew 
c;° > cj (Figure 3.3). If translocation of ions across the membrane is slow 
compared to the partition of ions across the interfaces, the interfacial distr. 
butions will approach those predicted by the Donnan equilibrium, where jj, 
is a function of concentration given by equation (3.30) so that B;°?  B/, In 


(DONNAN) 


Vigure 3.3. The origin of electrical potential differences across a 
membrane containing fixed anionic charges according to the é 
‘Teorell-Meyer-Sievers model. The fixed anionic charges are desig- 
nated by ©; mobile cations by @; and mobile anions by ©. 
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this instance 
2° =B°e,°; 7° =c_/6;°; &,'= Bie.’ and 
Clearly, @,° >¢,°, €° <c_°, @,)>c,! andti<ecl. 
As shown in Figure 3.3, the electrical potential difference across this mem 


brane will be the sum of two interfacial (Donnan) potentials and the diff 
sion potential arising within the membrane. Thus, 


Ay = yi - wy? =(v? a w°)+ (pi - w+ (i - y') 


pia esig! 


or 


Ay= |4 (in ar) +(y'- y?)+ 2" (in B| (3.31) 
where (W! - W) is given by the Henderson or Planck equations. 

This model for the diffusion potential across a charged membrane was ptw 
posed simultaneously by Teorell (1935) and Meyer and Sievers (1936) and 
later revised by Teoreil (1953). Although it is particularly relevant fo) 
charged membranes, the general approach applies to all membranes whery 
partition coefficients are concentration dependent, so that interfacial elev. 
trical potential differences do not cancel. 

A particularly interesting special condition arises when the fixed charge 
density of the membrane is much larger than the external salt concentration, 
that is, |z, X¥_'| >> c;. Under this condition equation (3.30) reduces to 


Bi? = (2, X_/c°)>>1 and 6/=@,X_/c')>>1 


Since C_° = c_°/B;° and t! =c_'/B/, as (z, X_) > ©, the concentration of the 
anion in the membrane phase will approach zero (so-called Donnan exclu 
sion). In the limit, diffusion of the salt across the membrane will cease (re. 
call that in the absence of an external current, J, must equal J_) so thu} 
(v' - Y°)= 0. Under these conditions 


Aw =(RT/F) In Bj} - (RT/F) In B;? 
or 
Ay = (RT/F) In(¢;2/¢;') 


which is the Nernst equilibrium potential (equation 3.7) of the system. |) 
this instance Ay is the sum of two equilibrium potentials that arise solely «! 
the interfaces of the membrane. This is the origin of the electrical potentla! 
differences generated by near perfect ion-exchange membranes and jon 
selective materials, for example, pH electrodes (cf. Eisenman, 1967). 
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Osmosis, van’t Hoff’s law, and Staverman’s reflection coefficient 
One of the seminal advances in physical chemistry during the nine- 
century was the recognition that solutes in dilute solution behave as 
ct gas; indeed, the extension of the perfect gas law to dilute solutions 
ts the foundation of the “thermodynamics of solutions.” This point 
illustrated by considering two rigid compartments at the same tem- 

separated by a rigid barrier. Assume that compartment o contains 
s of a gas in a volume V° and that compartment i contains n' moles 
game gas in a volume V!. According to the perfect gas law 


P°V°=RTn? ot P°=RTc? 


¢? = n°/V° and is simply the concentration of the gas in moles per unit 
_ Similarly, for compartment i we can write P* =RTc'. Thus, if the 
is nonporous, 


AP=P° - Pi=RT(c° - c!)=RTAc (4.1) 


is the pressure difference that will be sensed by pressure transducers in- 
into each of the compartments. If the barrier is porous, gas will diffuse 
the compartment with the higher concentration (pressure) to that with 
r concentration (pressure) until AP and Ac are abolished; thus, AP 
considered the “diffusion pressure” of the gas. 
ws now consider the system illustrated in Figure 4.1, where compart- 
» contains an aqueous solution of solute i at a concentration c;° and 
tment i contains a solution of the same solute at a higher concentra- 
¢/; the intervening membrane is assumed to be rigid and strictly im- 
ble to i but freely permeable to water. To the piston sealing compart- 
i, we now apply a pressure of sufficient magnitude to prevent any 
in volume of either compartment. Since in this system a change in 
can only come about as a result of water flow, the pressure applied 
be sufficient to abolish the driving force for water flow so that, from 
15, Aby = My! - Hw° =0, where wy is the chemical potential of 
In other words the state properties of the system must be such that 
is at equilibrium! 


au 
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Figure 4.1. System for studying osmotic pressure and osmotic 
volume flow. 


From equation (1 .27) we can write 
My? =(Myw 2 )° +0, P? + RT Inx,,? (4.2) 


where v,, is the partial molar volume of water (approximately 18 cm?/mole), 
P° is the pressure ira compartment 0, and x,,° is the mole fraction of water 
in compartment 0;* clearly, xy° =ny°/(n;? + n°). Similarly, for compart: 
ment i 


By! = (Uy? ! +5,P!+RT In x, (4.3) 


Equating (4.2) and (4.3), substituting the definitions of x,, and solving, 
assuming that (u,,° )? = (u,,°)', we obtain 


_ ; F nj nj? 
Dy(P! - P?) =RT {In 1+) le+—> (4.4) 
ny ny 
Since, in dilute aqueous solution, n,, >>n;, and since In(1 + x) approaches x 
as x approaches zero, equation (4.4) simplifies to 


rs ‘ > nj} nj? 
WW? Pep = RT SS (4.5) 
ny ny 
Since Dy ny'=V* and dyn,,° = V°, where V is the volume of the com: 
partment, equation (44-5) finally reduces to 


n= (P!- P?) =RT(c} - ¢,°)=RTAG; (4.6) 


which is the van’t Hoff equation for the osmotic pressure (7) across a barrier 
separating two ideal, dilute solutions when the barrier is strictly impermeable 
to the solute. Equation (4.6) is identical with equation (4.1) derived for per. 
fect gases. In this case pressure is applied to the more concentrated solution 
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if solute) which is obviously the side with the lower concentration of water; 
in strict analogy with the diffusion of gases, 7 can be referred to as the 
iffusion pressure of water.” 

| should be emphasized that the system we have just discussed is at equi- 
rium and 7 is the pressure that must be applied to the concentrated solu- 
n to offset, precisely, the difference in water activities across the barrier; 
this sense equation (4.6) is analogous to the Nernst equation (2.31; 3.7), 
re the electrical potential difference across the barrier precisely offsets 
difference in activities of a permeant ion. In both instances the equilib- 
m condition is achieved by a balancing of different but equivalent driving 
ees, 

ithe piston illustrated in Figure 4.1 is released so that (P! - P°) = 0 (ie., 
lh chambers are now open to the atmosphere), water will flow from o to 
ind the initial rate of water flow (i.e., before there are any changes in the 
“wentrations) is given by 


Dy Jw =Jy=Lyt (4.7) 


re L, is the hydraulic conductivity of the membrane; for a barrier perme- 
by right cylindrical pores whose diameter is much larger than that of 
ler, Lp is given by the Poiseuille equation. Equation (4.7) is an example of 
® linear relation between flows and forces discussed in general terms in 
pter 1. 
More than a half-century elapsed between the formulation of van’t Hoff’s 
and the next major advance in the theory of osmotic pressure and os- 
lic water flow; namely, the generalization of equation (4.6) by Staverman 
include membranes that are permeable to the solute (Staverman, 1951). 
Staverman reasoned as follows: If compartments o and i in Figure 4.1 are 
ed with the same aqueous solution with a concentration c;, then the appli- 
ion of a pressure to the piston will force volume from i to o, and the 
Wirafiltration permeability” of i can be calculated from the change in V° 
i) c;°. For example, if the membrane is impermeable to i, then only pure 
iter will be forced into compartment o and after a given period, the initial 
will decline to a value given by (c¢,°)’ =(c;°V°)/(V° + AV), where 
Y” is the increment in volume. Conversely, if the membrane is equally 
meable to solute and water, the solution in compartment i will simply be 
iisferred, unaltered, to compartment o and c¢;° will not change. Obviously, 
the membrane is permeable to i but less so than to water, the change in 
will fall between these extremes. The ability of the barrier to discriminate 
tween the solute and water can be described by a reflection coefficient o 
lined as 


0; = i fq, Rite), Mitra (4.8) 


te for small volume flows, c!**"4 is simply ¢;! and c;*™*° is the cal- 


lated composition of the solution that was forced into compartment o. 
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Clearly, if the membrane is impermeable to the solute 0, = 1 and, if the 
membrane is equally permeable to solute and water, o; = 0; for all intermedi: 
ate cases 0< 0; <1. Clearly, o; is a function of both the solute i and the 
membrane. The term “reflection coefficient” was chosen because it indicates 
the extent to which the solute particles are “reflected” by the membrane in 
comparison with the solvent. 

Staverman then showed that if cj >c,° the initial pressure that must be 


the case of the solute i, from equation (1.27), 
Ap; =0; AP+ RTA Inc; (4.13) 


Aline; = In c,? - Ine; = In(e;7/e#). 
order to linearize equation (4.13) (ie., eliminate the logarithmic term), 
ul K introduced the notion of a mean concentration ¢;, defined by 


applied to the piston to prevent volume flow is given by “a ( e? - of a ( Ac; mgs 
Test = 0; = 0,RTAc; (4.9) 7 \in (/c')/ — \In (c;?/c#') 

and the initial rate of volume flow from o to i when AP = 0 is g this definition, equation (4.13) simplifies to 
verano (4.10) A, = 0; AP + RTAc;/é;= 0; AP + 1/6; (4.15) 


where Teg is the effective osmotic pressure across a nonideal membrane. 

Equations (4.9) and (4.10) are completely generalized descriptions of os 
motic pressure and osmotic volume flow. 

As we will see, 0; is a “cross coefficient” that is the measure of the coupling 
(interaction) between the flows of solute and solvent across a membrane; and 
o; together with the two “straight coefficients” P; and Lp are sufficient to 
completely describe the flows of solute and solvent across a given membrane. 


, what is the meaning of ¢;? When is this important simplification 
Using the series expansion of In(c;?/c), it can be readily shown that 
the concentration difference is small so that Ac; > 0 and (c°/e?) > 1, 
6,7 (c;° + c')/2; in other words, when the concentration difference is 
. 6; is the arithmetic mean of the concentrations in compartments o 


i, é' 
stituting equations (4.12) and (4.15), into equation (4.11), we obtain 


4.2. Interactions between solute and solvent flows: T d;S/dt = J,(0,AP + n/¢;) + Jy 0y(AP - 7) (4.16) 
The Kedem-Katchalsky equations BI eacrangerment yields 
T d,S/dt = (J;0; + Jy0y) AP + [(Ji/¢:) - Jw8wI 8 (4.17) 


, (J;0; + JwOy) is simply the flow of volume J,. Further, since in dilute 
tion é,,0,, = 1, equation (4.17) can be written as follows 


T d;S/dt=J,P+Jpn (4.18) 


Jp = (Vie) - Jw/éw)| and is simply the difference between the 
ties of the oppositely directed flows of solute and solvent. Jp is thus a 
ire of the velocity of exchange of matter across the barrier. 
us, equation (4.11) has been transformed into a new expression of flows 
forces where AP is the conjugate driving force for the flow of volume Jee 
n or RTAc; is the driving force for Jp. 

* can now write the linear phenomenologic equations that describe this 
em as follows: 


Jy = Lp AP + Lpp RTAG; 


As implied in Section 4.1, when a membrane that is permeable to 
both solute and solvent separates two solutions having different concentra 
tions, there will be oppositely directed flows of solute and solvent, and as 
discussed in Section 1.6, in general, these two flows can interact. In 1958 
Kedem and Katchalsky (K and K) analyzed such systems using the formalism 
of irreversible thermodynamics; in one fell swoop they interjected this rela 
tively young discipline into the biological literature, and their contribution 
has had a profound influence on the conceptualization and investigation of 
membrane transport processes. 

As discussed in Section 1.6, in a system characterized by two flows and 
two forces, the “dissipation function” is given by 


@=Td;S/dt=J;X,;+JyXy (4.11) 


where, in the system we will consider, J; is the flow of an uncharged solute 
and J,, is the flow. of water. Therefore, X; is simply Ay; and X,, is Au,,. We 
will define the flows as positive when directed from compartment o to com 
partment i and Ay; = y,;° - wi! and Apy = By? - Hy’: 

From equations (4.2) to (4.6) it is a simple matter to show that for a dilute 
solution where n;/n,, > 0 and n,vy > V 


Au, =0,,(P? - P!)- 5, RT (ci? - cj!) = Dy(AP - 1) (4.12) 


(4.19) 
Jp = LppAP +LpRT Ac; 


fe, according to the Onsager relation, Lpp = Lpp- 


——= 
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tuting the definitions of Jp and J,, we obtain 


(euien- Ulé)  a-o 


Ji 0; i Jwiw 


Although equations (4.19) seem abstract, they, in fact, more faithfully 
describe the actual experimental system than does equation (4.16). Thus, 
when Ac; = 0 and a pressure is applied to compartment i, we measure the 
rate of increase in the volume of compartment o or J, (not J; or J,,); the 
ratio (J,/AP),=9 is the hydraulic conductivity L,. We can also measure the 
change in composition of compartment o due to the fact that the applied 
pressure will, in general, force a fluid across the membrane that does not have 
the same composition as the solution in compartment i. This differential 
effect of an applied pressure on the flows of solute and solvent is commonly 
referred to as ultrafiltration; Jp is a measure of this differential flow and the 
ratio (Jp/AP),7=9 =Lpp is a measure of the ultrafiltration properties of the 
membrane. 

Conversely, if AP = 0 and c;° # c/, there will be a displacement of volume 
because of the flow of water (osmosis) from the compartment with a lower 
concentration of i to the other and an oppositely directed flow of volume 
caused by the diffusion of i down its concentration gradient. What is actually 
measured is the total displacement of volume, and the ratio (J,/7)a p=o =Lpp 
is a measure of the relative contributions of solvent flow and solute flow ty 
the total displacement of volume when AP = 0. 

If the system is not too far displaced from equilibrium (i.c., if the flows and 
forces are small) Lpp = Lpp, so that 


(ol ™)ap=o = Up/AP)n=0 (4.20) 


Thus, osmotic volume flow in the absence of a hydrostatic pressure is equal 
to the “ultrafiltration effect” in the absence of a concentration difference 
This is precisely the reasoning behind the approach employed by Staverman 
to determine the reflection coefficient o;. Further, from equation (4.19), we 
see that when Ac, #0, then J, = 0 when 


(AP);, =o = ~(Lpp/Lp) n (4.21) 


Comparison of equation (4.21) with equation (4.9) reveals that 9, © 
~UZpp!Lp). 

Finally, when AP=0 and Ac; #0, there will be an exchange flow of sol. 
vent and solute given by Jp; the membrane property that describes this ex. 
change is (Jp/T)ap-o =Lp. 

We now examine the relations that prevail if the membrane is ideally sem). 
permeable, that is, permeable to solvent (water) but not solute (i). Clearly, 
from equation (4.9), 6; = 1 so that -Lpp = Lp. In addition, under these con 
ditions Jp = -(J,/¢,,) = -(Jy d,) = -J,. Since (Jp/Tapeo = Lp, it follows 
that Lp = Ly. Thus, for a semipermeable membrane, -Lpp = Lp = Lp. 

What if the membrane cannot discriminate between solute and solvent’ 
From equations (4.19) we see that when Ac; = 0 


(Jp/Jv) =o os polLp 0) (4.22) 


ce in dilute solutions Jy0,, >>J;0;, equation (4.22) reduces to 


pete) 1-4, (4.23) 
Jié, = v; and J,,/¢,, = Yy, where vy and v,, are the velocities of solute 
lvent molecules, respectively, (see Section 2.1). Thus, 


(0;/0w)n=0 = 1- % 

ly, if the barrier cannot discriminate between solute and solvent, there 
be no “ultrafiltration” in response to an applied pressure so that u; = vy, 
= 0. (Needless to say, if the barrier is ideally semipermeable, uv; = 0 and 


Be ations developed here lead directly to two expressions that com- 
ly describe the interactions between solute and solvent flow across a 

neous membrane that separates two solutions. The introduction of 
finition of o; into equation (4.19) yields 


Jy = Lp(AP - 0;RTAC;)) (4.24) 


at J, is directly proportional to the difference between an applied hydro- 
pressure and an oppositely directed effective osmotic pressure. In gen- 
, when there are 7 solutes 


n 
Jy =Lp (ar- > oRTAc;) 


her, solving equation (4.24) for AP and substituting this expression into 
‘dquation for Jp (4.19) yields 


J, = (1 - 9;) CjJyt w,RTAc; (4.25) 


w= 6(LpLp ~ Lpp” /Lp. Since LpLp is always greater than Lpp* 
jon 1.6), w; is always positive. Further, when J, F: 0, Vi)s,=0 = 
TAc;. Thus, «9; is a measure of solute permeability and is related to the 
familiar permeability coefficient by P; = w,RT. : 

en Ac; =0 then J; =(1 - 0;) ¢;J, or from equation (4.24) Ji =(1 - 9;) 
L, SP. Thus, if o;# 1, there will be a flow of solute that is directly pro- 
rtional to the flow of volume or, alternatively, the applied hydrostatic 
ure (AP). This component of solute flow is frequently referred to as 
Ivent drag,” implying that it arises from an entrainment of solute particles 
4 flowing stream of solvent. However, as we will see, this interpretation 
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implies that solute and solvent share a common pathway through the mem, 
brane (e.g., a pore) and need not be correct, 
When AP=0, -J, =L,0;RT Ac; and, substituting this expression int 


6; is the partition coefficient described earlier (Sections 1.5 and 2.2) 


equation (4.25), we obtain (c) o=1- Qi _ Wifiw AX (4.29) 
: Yen Pw 


Ji= (-C ~ 9) oj Lyd; + wJRT Ac; (4.26) 
expressions for L, and w, have relatively straightforward interpreta- 
and are intuitively satisfying (however, see Section 4.5). The expres- 
for 0;, however, provides important insight into the meaning of this 
coefficient and warrants some discussion. 

der a lipid, porous membrane that separates two aqueous solutions 
solute 7. Assume that water can cross the membrane only through the 
and that i, a large lipophilic molecule, is excluded from the pores but 
diffuse across the membrane through the lipid matrix. Under these con- 
fiw = 0 and equation (4.29) reduces to 


0; = 1 - (w,0;/Lp) (4.30) 
AP=0, then J, = -Lpo;RT Ac; [equation (4.24)] so that 


43. Physical interpretations of Lp,w, and o 


It should be clear from equation (4.24) that L, is a straight coeffi 
cient that relates the rate of volume flow (which in a dilute aqueous solution 
is essentially equal to the rate of water flow) to the conjugate driving force 
namely, the effective pressure difference across the membrane. Similarly, 
from equation (4.25), we see that w), is alsoa straight coefficient that relates 
the diffusional flow of an uncharged solute to its conjugate driving force A¢, 
In contrast o; is a cross coefficient that provides a measure of the interaction 
between the flows of solute and solvent and nonconjugate forces. It should lye 
emphasized that these coefficients are strictly phenomenologic and are ind. 
pendent of underlying mechanism; that is, they are not based on any partiew 
lar model of the membrane or the flow processes. Although this assures 
general validity, it is intellectually unsatisfying and, over the years, sever! 
attempts have been made to provide these coefficients with intuitively satiy. 
factory physical interpretations. The most successful approach is based on 
the assumption that the generalized phenomenological resistances or condu¢ 
tances can be translated into mechanical frictional coefficients. Thus, it is 
argued that when a particle achieves a steady-state velocity (i.e., it is not 
accelerating) the total force acting on the particle is zero and the thermo. 
dynamic driving force X; is precisely balanced by opposing frictional forces 
Thus, we can write X;=-2, F,;, where the summation includes all possibly 
frictional forces encountered by the moving particle. The frictional forces 
are given by Fi; = -f,;(v; - 7), where Ji; is a frictional coefficient and (v; - ¥) 
is the difference between the velocities of the solute 7 and some Eonigonl / 

in the system with which it interacts. For the case of the flows of an un 

charged solute and water across a membrane, the membrane is generally con 

ee to rs the stationary reference for the flows so that Up», = 0 and the 
our frictional coefficients that describe the system are ‘ , 

Swm; Where fj, is the frictional coefficient bce re ee eRe p fn” cnihamnilt eke 
forth. Using this approach, Kedem and Katchalsky (1961) have derived the , fecall that Ay, =0,4P+RTA Inc; [equation (4.13)] so that when 
following relations: © 0, Au;=0, 4P. Thus, equation (4.32) simply states that the flow of i 
product of a mobility term, a concentration term and a driving force 
Section 2.1), where the latter is derived solely from the applied pressure. 
rly, under these conditions J; is directly proportional to AP and J,, 
tla that are often applied to infer the presence of “solvent drag.”” How- 
1, obviously, there is no entrainment of solute in a solvent stream so that — 
inference is entirely unjustified. 


-Jy = L),RT Ac; = 0,0;RT Ac; 
(4.31) 
= J. = L,RT Ac; tt 0,5; 


, the first term on the right of equations (4.31) is simply the volume dis- 
ment caused by osmotic water flow and the second term is simply the 
of solute that diffuses across the membrane driven by Ac;. Obviously, 
or P;=0, equation (4.31) reduces to the van’t Hoff expression for an 
ly semipermeable membrane. When w;>0, the deviation from ideal 
ior is solely because of the fact that the ideal osmotic volume flow is 
d by the volume flow of the permeant solute in the opposite direction. 
insider the situation when Ac;=0 and a hydrostatic pressure AP is 
across the membrane. Under these conditions 


Jy =L)AP and J;=(1- 0;) cjJy 
g equation (4.30) for Z, and substituting, we obtain 


(a) Lp = Swi /fiym OX (4.27) 


‘a Yw is the volume fraction of water in a membrane whose thickness |y 
we 


(b) w= Bi/(Siw + fim) Ax (4.28) 
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Referring once more to equation (4.29), we see that when the permeant 
solute and water traverse the membrane through a common path, there are 
two factors that contribute to the decrease in volume flow from that pre- 
dicted by the van’t Hoff equation for an ideal semipermeable membrane. The 
first is the oppositely directed volume flow of solute and the second arises 
from the frictional interaction between the oppositely directed flows of sol- 
vent and solute. Since o,;, w;, Lp, and 0; can be determined experimentally, 
equation (4.29) provides a criterion for determining whether there are inter- 
actions between solute and solvent flows across the barrier. Thus, if 


0; <= (2,0; /Lp) 


then f;,, > 0 and a common path may be inferred. 

Finally, according to equation (4.29), o; decreases with increasing w; and 
may assume negative values when w, is sufficiently large. This situation is 
referred to as a negative or anomalous osmosis, inasmuch as J, will be directed 
from the compartment with the higher solute concentration to the compart- 
ment with the lower solute concentration and if cf >c,° (Figure 4.1), a 
negative pressure would have to be applied to the piston to prevent volume 
flow; this apparent anomaly, which has been observed, is accounted for by 
the fact that volume flow of the diffusing solute as well as that resulting from 
the frictional interaction between solute and solvent overweighs the nomnally 
directed osmotic flow (Grim and Sollner, 1957, 1960; Talen and Staverman, 
1965; Lakshminarayanaiah, 1969). 

In summary, this frictional analysis permits us to express phenomenological 
coefficients in terms that provide some intuitive feeling for their bases. The 
principal shortcoming, however, is that only f;,, and f,,; can be evaluated in- 
dependently and, then, only under special conditions. Because f;,, is the 
frictional coefficient between a mole of solute and an infinite amount of 
water, it can be estimated from the diffusion coefficient of i in free solution 
using the Einstein relation: 


D? = RTu; = RT[fiw? 


(4.33) 


where the superscript o denotes that the process is taking place in free solu. 
tion. If we assume that the frictional interaction between solute and water in 
the membrane is equal to that in free solution (e.g., if i diffuses through pores 
that are large compared to the size of the solute and solvent molecules), then 
fiw? =Siw. Further, as discussed by Kedem and Katchalsky (1961), f,,; (the 
frictional coefficient between 1 mole of water and an infinite amount of 
solute) is equal to ¢;f;,/¢y. Thus, under the unlikely condition that water. 
solute and solute-water interactions in the membrane do not differ markedly 
from those in free solution, fj and f,,; can be approximated. 

- On the other hand, there are no generally acceptable approaches that per. 
mit independent estimates of fj, (the frictional coefficient between 1 mole 
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of i and a membrane of infinite thickness) or f,,,, from first principles; con- 
pequently, although the frictional analysis seems intuitively reasonable, it 
femains phenomenologic. 


The Renkin-Durbin model 


Renkin and Durbin have suggested a somewhat different approach 
for the physical interpretation of o that is also based on hydrodynamic rea- 
foning but is couched in terms that are physically more meaningful than 
those of the frictional approach. Renkin (1954) argued that the ultrafiltra- 
tion of an uncharged particle through a membrane permeated by right cylin- 
drical pores could be described by the combined expressions of Ferry (1936) 
and Faxen (1923) and derived the expression 


(Aig/Ap) = 2[1 - Gi/rp)]? - (1 - (ri/rp)\* 
+ [1 - 2.104(7;/rp) + 2.09(7;/rp)® - 0.95 (ri/rp)°] 


where (A;7/A,) is the effective area for ultrafiltration divided by the total 
pore area, and 7; and 7, are the radii of i and of the pore, respectively. As 
uussed in Section 2.7, Aj includes the restrictions upon the entry of i into 
the pore as well as the hindrance to the flow of i through the pore. 
A similar expression can be written for the filtration of water (A ,,,/A,) by 
‘wmply substituting 7,, for 7;. 
Durbin (1960) later argued that 


o; = 1 - (Ajz/Awy) 


‘and, therefore, can be defined entirely in terms of 7;,7,,, and 7p. Studies using 
‘uttificial membranes with pore radii estimated as ranging between 15 A and 
‘$2 A are reasonably consistent with this analysis (Durbin, 1960). However, it 
should be emphasized that the smallest pore radius is approximately 10 times 
that of water (1.5 A), and it is not clear whether these hydrodynamic models 
apply when the pore radius approaches the dimensions of the permeating 


yolecules. 


Electrokinetic phenomena 


44. 

In Section 4.3 we considered the interactions between the flow of 
solvent (water) and that of an uncharged solute i across a homogeneous mem- 
rane. We now briefly consider the case when i is a completely dissociable 
salt, the approach will follow that of Kedem and Katchalsky (1963a), but for 
the sake of simplicity we will only consider a salt that dissociates into a 
Monovalent cation (+) and a monovalent anion (-) [see also House (1974) 


for a similar approach]. 
The experimental system we will refer to is illustrated in Figure 4.2. The 


68 Basic principles of membrane transport 


PISTON 


LEACEE EES, | PTE, ETE LLL TOD 


PRESSURE 
TRANSDUCER 


Figure 4.2. System for studying electrokinetic phenomena. The 
electrodes are reversible to the anion so that the voltmeter & records 
the electromotive force of the anion given by equation (4.38). 
Current J can be passed across the membrane from the external 
electromotive force (battery) B. 


apparatus permits measurement of J, and AP. In addition, electrodes inserted 
into each compartment connected to a high resistance electrometer (volt 
meter) permit measurement of the electrical potential difference across the 
membrane Ay under conditions of “essentially zero” current flow. And, 
identical electrodes connected to a variable electromotive force (battery) are 
used to pass current J across the membrane. Thus, we can measure AY when 
T= (as well as the “external current” needed to “clamp” Ay to zero. For the 
sake of convenience and in order to adhere to other treatments of this subject 
in the literature, in this section we define Ay = ~° - w! and continue to 
define Ac; = ¢,° - c/, AP=P° - P!, and Au, = u;° - pj. 

The dissipation function for this tertiary system (i.e., three species, three 
flows, and three forces) is 


T d,S/dt =J, Afi, +J-Of-+ Jy Spy (4.34) 
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Aj, =0,4P+RTAInc,+ F Ay or, using the linearizing approxima- 
discussed in Section 4.3 [equation (4.14)] 


Afi, =0,4P + (RTAc,/é,)+ FAW (4.35) 
rly, 
Oj_=0_MP+(RTAc_/¢.)- FAy (4.36) 


ce bulk electroneutrality must be maintained in both compartments, it 


0.2 =c°=c? and c=ct=e/ 
t the difference in chemical potential of the 1:1 salt is 
Ay; = Api, + Ap_= u,AP+ 2(RT Ac;/¢;) (4.37) 


0; = 0, + 0_. Equation (4.37) simply states that the chemical potential: 
4 dissociable salt is related to the concentration of individual particles 
ing from dissociation (i.e., a colligative property of solutions). 
ow, the expressions for Ap, and Aj. given by equations (4.35) and 
) could be substituted into equation (4.34) to give an explicit statement 
the dissipation function. However, it is more convenient to define the 
romotive force of the system, &, which can be measured by inserting 
tical electrodes reversible to one of the ions in each solution and measur- 
the electrical potential difference between these electrodes potentio- 
trically (ie., using a very high resistance voltmeter (electrometer) that 
very little current). The observed electrical potential difference is the 
ntial driving force (energy) for the flow of that ion through the external 
t and, if the electrodes are reversible for the anion, is given by (Katchal- 
and Curran, 1965, pp. 149-50) 


& = -(Ai_/F) = ~(RT/F) Ain c_+ AW (4.38) 


units of volts. Clearly, in any system under steady-state conditions, & bears 
instant relation to Ay. 

Now, if the electrodes used to measure the electrical potential difference 
to pass current across the membrane are strictly reversible to the anion, 
cation is neither withdrawn from nor released into the bathing solutions 
the electrode reactions,? so that J, across the membrane can be equated 
the net displacement of the undissociated salt i, and J, = Jj.- 

Further, the current that must be passed across the membrane in order to 
p & to zero is the difference between the transmembrane flows of 
tion and anion which is equal to the flow of anion (electrons) from one 
Mpartment to the other through the external circuit.2 Thus, the current 
ugh the membrane is given by 


I= F(J, - J.) (4.39) 
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Substituting equations (4.37), (4.38), and (4.39) into equation (4.34) yields 
a new expression for the dissipation function given by 
T d,S{dt=Jy Ap, + J; Au; + 1& (4.40) 


It follows that we may describe this system by means of the following three 
linear phenomenological equations: 


Jw = Law Ab + Ly Oyj t Ly & 
Jy =Liw Atty + Ly Au; + Liz 
T= Ly Spy + Ly; Ou; + Ly& 


(4.41) 


These equations contain nine phenomenological coefficients; but if On- 
sager’s reciprocal relations are obeyed, Ly; =Ljy, Lyz = Lyy, and Lj = Lz; so 
that the number of independent coefficients is reduced to six. Thus, in gen- 
eral, the total water flow is comprised of a hydraulic component, an osmotic 
component, and a component coupled to the flow of current (or, alterna- 
tively, driven by the electromotive force), which is referred to as electro- 
osmosis. Similarly, the total flow of salt is made up of a fraction resulting 
from ultrafiltration, a fraction resulting from diffusion, and a fraction driven 
by & that is referred to as electrophoresis. Finally, the total current is made 
up of a streaming current coupled to the flow of water, a diffusion current 
driven by the concentration difference of the salt and the (electrical) current 
actually driven by the electromotive force. 

Although equations (4.41) completely describe this system, as discussed in 
Section 4.2, it is frequently more convenient to deal with volume flow, 
which is actually measured, rather than water flow. Recall that Auw = 0,, ° 


(AP - RTAc;) and Ay; =; AP + RT Ac;/¢;, so that equation (4.41) can be 
written 


wae =J,(AP - RTAc,;) + a +J& (4.42) 
and since in dilute solution ¢;0; << 1, equation (4.42) reduces to 

Ae =J,(AP - RTAc;) + ane +1& (4.43) 

The corresponding phenomenological equations are 

J, =Ly (AP - RTAc;) + Ly (RT Ac;/¢;) + Ly & (4.44) 

J, = Lj OP - RTAc;) + Ly(RT Ac;{e;) + Ly & (4.45) 
and 

1=Ly(AP- RTAc;,) + Ly(RT Ac; /¢;) + Ly& (4.46) 


where, once again, Ly = Ljy, Lys = Lyy, and Ly = Ly. 
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e problem now is to translate the straight and cross coefficients in equa- 
(4.44) to (4.46) into experimentally useful or practical terms that 
¢ the properties of the system. This can be accomplished by examining 
flows under experimental conditions when two of the three driving forces 
lows are set equal to zero. Clearly: 

When Ac; = 0 and & =0 (so that Ay = 0), equation (4.44) reduces to 
Ly AP, so that Lyy is obviously equal to L,. Under the same conditions 
» Liy AP, so that from equations (4.24) and (4.25) Liv = (1 - 6;) Clip = Lu 
) When (AP - RTAc;)=0 and & = 0 but Ac; #0, then 


J, = Ly (RT Ac;/¢i) 


) When AP=0 and Ac; = 0, then /= L7& = Lz Ay so that Lzis simply 
“electrical conductance of the membrane Gy. Under the same condi- 
Jy = Ly & = LyzllGm.- Thus, we can define a new cross coefficient 
© (Ly/Gm), which is the volume flow associated with unit current flow 
n AP=0 and Ac; = 0 or the electro-osmotic permeability. Finally, under 
conditions J; = Liz& = Liz/G,,. Since the current carried by J,, which 
recall is equal to J;, is simply FJ,, we can define the transference mae 
., the fraction of the total current carried by J, when AP = 0 and Ac; = 0) 


t, = (FI) = F Lisl Gm 

‘Thus, we have converted the three straight and three cross coefficients into 
new set of experimentally meaningful coefficients, where Lp, «;, and Gy, 
the corresponding straight coefficients and 0;, ¢,, and 8 are the cross 


fficients. 
We can now rewrite equations (4.44) to (4.46) as follows: 


Jy = Lp AP - ojLpRT Ac; + we (4.47) 
J, = Lp(1 - 9;) &(AP - RTAci) + w, RT Ac; + te/F (4.48) 
1 = B'Gm(OP - RTAG) + GmteRT ACFE; + Gm& (4.49) 


i hat formidable, 
Now, although equations (4.47) to (4.49) appear somewhat f 
they provide insight into a number of well-described electrokinetic phenom- 
ena. For example, if both bathing solutions have identical compositions so 
that Ac; = 0 and & = Ay, the equations describing volume and current flow 
are simply 


Jy=LpOP+6'l (4.50) 


and 
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Thus, when-a current is passed across the membrane, a pressure must be 
applied by the piston to prevent electro-osmotic volume flow; this pressure, 
referred to as the electro-osmotic pressure, is simply. 


Pe =-B'l/Ly (4.52) 


If a pressure is applied across the membrane, producing a volume flow 
under zero current conditions (i.e., in the absence of an external current 
source), a streaming potential will be generated given by 


-Ay = B'AP= 6'J,/Ly (4.53) 


and the streaming current (i.e., the external current necessary to clamp Ay to 
zero in the presence of volume flow) is 


1=B'G, OP (4.54) 


Finally, by simply manipulating equations (4.50) and (4.51) the following 
identities emerge: 


(a) (/AV)ap-o = (U/AP)ay =o (4.55) 


so that the relation between electro-osmotic volume flow and Aw when 
there is no applied hydrostatic pressure is equal to the relation between the 
(streaming) current necessary to reduce Ay to zero and an applied pressure. 


(b) UMv)ay -0 = -(AP/AY)z, =o (4.56) 


so that the relation between the streaming current and a pressure-induced 
volume flow is equal to the ratio of the electro-osmotic pressure necessary to 
prevent volume flow and the applied electrical potential difference. 

These relations have been experimentally confirmed repeatedly and sup- 
port the validity of applying Onsager’s reciprocal relations to systems in 
which the flows of solvent and electrolytes are coupled (Miller, 1960). 

Finally, it should be noted that in systems characterized by at least two 
flows and two forces, the straight and cross coefficients can be experimen- 
tally defined in several ways. For example, the electrical conductance of the 
membrane G,, can be determined under a variety of experimental conditions. 
Obviously, one desired condition is when Ac; =0 so that & = Ay. However, 
even under this condition, we can determine Gm, or 0J/8(Ay), either when 
J, = 0 or when AP=0. The values obtained will, in general, not be equal 
since 


(Gm)ap=o =(Gm)y,=0/(1 + Pah’) 


where Pg =(AP/AwW) when J, and Ac; are zero; Pg is, therefore, a coeffi- 
cient of electro-osmotic pressure. Thus, there is no “unique” electrical con- 
ductance, and the choice of values is a matter of experimental necessity or 
convenience. 
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Figure 4.3. A composite membrane made up of different regions 
arranged in parallel. 


The same principle applies to all of the phenomenological coefficients. The 
important point is that they are defined for a given set of experimental con- 
ditions, and their values can only be applied under the same set of conditions. 


45. Transport of solute and solvent across composite membranes 


Up to now we have been concerned with interactions between the 
‘flows of solute and solvent across homogeneous membranes that have uni- 
form properties throughout and can be characterized by a single set of ele- 
ntal straight and cross coefficients. We will now consider inhomogeneous 
or composite membranes that, because of interactions between flows, can 
exhibit properties not seen in simpler systems. We will limit our formal con- 
siderations to the flows of an uncharged solute and solvent and simply point 
Out, where appropriate, the behavior that can be observed with charged 
wpecies. The approach will closely follow that of Kedem and Katchalsky 
(19635, c) and House (1974). 

Let us first consider the membrane illustrated in Figure 4.3, which is made 
up of two regions a and 6 that have different properties and are arranged | 
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parallel to one another. In such a system the forces acting across each region 
are the same, and the total flow of any substance is simply the sum of the 
parallel flows across each region. Thus, we can write for the jth flow 


Jp= yf +d (4.57) 
where the y’s are fractional areas, and 
a 
“> Ljx'X, and Jj? = > Lj? Xy (4.58) 
k= 


k=j 


where aj is expressed per unit area occupied by regions a and ); clearly, 
y7+7° =1. The set of elemental straight and cross phenomenologic coeffi- 
cients Lik’ and Lix? relate the flows across each region to the sane set of 
driving forces; if the two regions have different properties, they must differ 
in one or more of these coefficients. 

Now, the total observable flow is given by 


n 
Jj = > LieXn (4.59) 
k=j 
where the set of coefficients designated by Lj, are the overall observable 
coefficients and are some combination of the elemental coefficients. Sub. 
stituting equations (4.58) and (4.59) into equation (4.57), we see that each 
overall coefficient is simply the area-weighted sum of the borepmponding 
elemental coefficients, or ; 


Lie = Lj? + YL in? (4.60) 


This is not surprising in view of.the fact that the L’s are generalized con. 
ductances, so that strict additivity is analogous to Kirchhoff’s law of conduc- 
tance in a parallel electrical circuit. 

Now, the interacting flows of solute i and solvent w across region a can be 
described by equations (4.19) as follows: 


=Lp*AP + Lpp*RT Ac; (4.61) 
and 
Jp? = Lyp* AP + Lp*RT Ac; (4.62) 
and a similar set of equations can be written for region b. 
It follows that 
Lp = Lp’ + y°Lp” (4.63) 
Lpp = ¥*Lpp* + y’Lpp” (4.64) 


and 


Lp =y*Lp* + y°Lp” 


(4.65) 
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Equation (4.63) simply states that the overall Lp is the area-weighted sum 
of the elemental Z,’s; in view of the fact that L, is a pure conductivity, this 
nclusion is not surprising. 

Since 0; = -(Lpp/Lp), it follows that 


0; = -(YLpp* + Y’Lpp’ Lp 


o;'= ={y* pp '|Lp) /. (y°Lpp Ee) 
plying the first term on the right by (L,4/L,°) and the second term by 
ip IL”), we obtain 


0; = 0;°(y"Lp[Lp) + 04°(y"Ly /Lp) (4.66) 
Now, when Ac; = 0 anda presuiee difference AP is applied across the mem- 
ane, J, = L, AP, J, =L,° AP, andJ,° =L,” AP. Thus, 
(1e") (“") r (zie) ee 
|] = | an = 
L, Jy Lp Jy 
1 follows that the overall o; is the sum of the elemental reflection coeffi- 
ants of the two regions weighted by the fractional volume flow through 
ese regions in response to an applied pressure. 
Kedem and Katchalsky (1963b) have shown that for electroyltes there is 
| additional term in the expression for o; that may either be positive or 
wgative, so that the overall reflection coefficient may be greater than or 
jaller than either of the elemental o’s. 
discussed in Section 4.3, w;=¢;(LpLp - Lpp”)/Lp, so that using equa- 
(4.63) to (4.65) it can be readily shown that 
a,b; ay Bb 
16k b 
= 7707 + y?u9? + (0,4 - 0;°)? (Pte Ci 
P 


(4.67) 


Wiis, unlike the cases for L, and o;, w, is not simply an approg nately 
ighted combination of he elemental w’s except when 0,7 = 0,2; if this 
ity does not hold, the overall ww, will always be greater than the area- 
ighted sum of the elemental w’s. 
What is the origin of this significant departure from the behavior of a 
jomogeneous membrane? Recall that w, is defined as (J;/RTAc;) when 
|, © 0. Thus, in order to determine w;, we must apply a hydrostatic pressure 
© compartment i of Figure 4.1 sufficient to abolish J, (for the sake of this 
iration, we will continue to assume that ¢/ > c;). The required pressure 


7°L a yL b 
AP = 0,RT Ac, = ly, of! + >) a |araq (4.68) 


P P 
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Under these conditions 


Jy = Sy? + JS,” =0 (4.69) 
Now, from equation (4.24) 
"Jo" = 11, (AP - of RTAc;) (4.70) , 
c! 
and i 
yy? = 7? Ly (AP - o,°RTAc;) (4.71) 
Substituting equations (4.63) and (4.68) into equation (4.70) and (4.71), we 
obtain 
ay a,b; b 
Lb 
oJ a; (Ae) (0? - o*)RT Ac; (4.72) 
Lp 
and 
aya by b 
rat os = € p Y Lp ) (o? - o*)RT Ac; (4.73) 
ip 


Thus, when J,=0 and o,7=0;", there will be no volume flow across 
regions a and b. But when o;,* # 0,”, the condition of zero overall volume 
flow will be achieved by the presence of equal and oppositely directed vol: 
ume flow across the two different regions; in other words, a fluid circuit will 
be established. Further, each of these fluid streams may be associated with a 
coupled solute flow so that there may also be a recirculation of solute across 
the membrane. The direction and relative solute content of each stream can 
be deduced intuitively. If o,? is greater than o,7, the fluid stream through 
region b will be directed toward compartment i (i.e., toward the compart: 
ment with the higher concentration), but will be relatively poor with respect 
to solute content. On the other hand, a fluid stream of equal volume will 
flow from compartment i to compartment o through region a, which will be 
richer in solute content. Thus, the difference in solute contents of the oppo: 
sitely directed fluid streams contributes to the flow of i from compartment i 
to compartment o (which is assumed to be strictly diffusional when J, = 0), 
The total effect is to augment the rate of appearance of ¢ in compartment ( 
and, thereby, to exaggerate the apparent permeability of the membrane to /. 

This phenomenon can be described formally by writing the appropriate 
equations for solute flow, that is, 


Figure 4.4. A composite barrier made up of two different regions 
arranged in series separated by an intermediate zone that can be 
considered infinitesimally thin. 


bstituting y°J,° =-y7J," into equation (4.75) and summing equations 
74) and (4.75) yields 


J = 7,207? - 0/2) 6; + (yw ;? + yw) RT Ac; (4.76) 


arly, if 0,2 >0;*, J,” is negative (directed from compartment i to com- 
tment 0) and J; is augmented by the differential solute flow.* Alterna- 
ly, if o,? <o;7, then J, is positive (i.e., it is directed from compartment 
to compartment /) and, once again, J; is augmented. From equation (4.76) 
4s clear that (when J, = 0) 


(Jy /RT Mey) = 2; = (7707 + y?w;”) 


when 0,7 = g,?! 


Let us now tum our attention to a composite membrane made up of two 
erent barriers arranged in series separated by an infinitely thin intermedi- 
region m as shown in Figure 4.4. The relations between flows and forces 
i define such an array are 


x; =X/ +X)? (4.77) 


If = 77S - 6;7) 6; + y7wART Ac; (4.74) 
and 


Si? = yPI°(1 - 0;°) 6, + y?exPRT Ac, (4.75) 
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where X* = X;° - Xj" and X;? =X;" - X},and 
n n 
Tp> Yo LEX * =) Lin? Xa? p*) 
kaj k=j 


Equation (4.77) asserts that the sotal jth force across the system is the sum 
of the forces across barrier a and barrier b and is simply a statement of the 
principle that the difference in intensive properties of state between two 
points is independent of path. Equation (4.78) asserts that under steady-state 
conditions the flows across each barrier arrarmged in series are equal; this is 
simply a statement of the conservation of matter. 

By definition, w, is given by (J;/RT Ac;) when J, =0. Under these condi- 
tions, according to equation (4.78), 


J, = w'RT Ac? = wPRT de; = w,RT Ac; (4.79) 


where Ac,“ =c;° - cj" and Ac? =¢/" - c/ so that Ac; = Ac;* + Ac;?. It 
follows that 


Ji{(1/wj?) + (1 /eo;”)] = RT AG; 
so that 
(1/9;) = (1/erg*) + (1/0,”) (4.80) 


Since (1/w;) is simply the overall resistance of the barrier to solute diffusion, 
equation (4.80) is a statement of the strict additivity of resistances arranged 
in series. 

Similarly, when J, = O, AP=0;,RT Ac; and the analogous relations across 
the two membranes are AP* =0,RT Ac;? and AP? =0;RT Ac,?. From 
these equalities and those given in equation (4.79) it follows that 


G; = 0f£(w;/1;%) + 0?(c24/@;”) (4.81) 


Thus, the overall reflection coefficient is the surm of the two elemental coeffi. 
cients weighted by the relative resistances of the two barriers to the diffu- 
sional flow of solute. 

Equations (4.80) and (4.81) have several interesting implications. 

First, if either w,;? or Kaa? is dependent upon concentration, the system 
illustrated in Figure 4.4 will exhibit an asymmetry with respect to “perme. 
ability.” That is, the «2; measured with concen tration c;° in compartment o 
and ¢;' in compartment i will differ from that measured when the solutions 
are reversed. Further, under these conditions o; will also display a polarity 
{equation (4.81)], so that the effective osmotic pressure across such a barrier 
for a given Ac; will depend upon the orientation of the solutions. 

But, perhaps the most interesting and unexpected property of the series 
array illustrated in Figure 4.4 is the relation between the overall L, and the 
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slemental L,’s. Since Lp is given by (J,/AP) when Ac; = 0, we can write 

Jy = Lp (AP* - oT Ac;*) = Lp*(AP? - o;°RT dc;”) (4.82) 
re, according to equation (4.77) 

AP=AP*+ AP? (4.83) 
ce when Ac; = 0 
RT Ac# = -RT Ac;? (4.84) 


b 
—~_ = AP - (o;* - o;°)RT Ac;* (4.85) 


Thus, in spite of the fact that the L,’s are conductances, the overall L, is 
f simply equal to [e,*b AL + Ly”)| as would be expected for the. sum- 
tion of electrical conductances in series unless 0,7 = 0,” or Ac;* = 0. 

Now, Ac;* in equation (4.85) can be evaluated from the following relations: 


Jy= (1 - 62) Jy 6j + wART Aci? = (1 - 0/°) Jyé? + wPRT Ac; (4.86) 
re ¢;" and ¢,? are the average concentrations across barriers a and b, 
ctively, given by equation (4.14). Clearly, when Ac; = 0, ¢;* = ¢,? ; thus, 
m equations (4.84) and (4.86) it follows that 
b 
<< Js (4.87 
RT Ac; Cj ie + wy? v ) 


Substituting equation (4.87) into (4.85) and rearranging, we obtain 


Lp" + Lp” erty - 4.88 
? {( £*Ey” ) t fae +w? v Ri ree 
apAesf Ae\ ofa | oot Var 4.89 
re) = (ca) + Gs) + ret Gj (4.89) 


According to equation (4.89), the resistance to volume flow across this 
jamposite membrane will always exceed the sum of the resistances of the two 
Hes barriers unless 0,7 = 0;”. 

The physical reason for this intuitively unexpected property arises from 

behavior of the central compartment. Thus, let us assume that Ac; = 0 
di that a positive pressure is suddenly applied to compartment i. If o,? <a;*, 
fluid stream crossing barrier 6 from compartment i to compartment m 
contain more solute than the stream that flows from compartment m to 
partment o across barrier 4. Consequently, solute will accumulate in the 
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intermediate compartment. In time, a steady state will be achieved in which 
cj" >c;° and the effective pressure generated by the concentration differ- 
ence (Ac;’) across the less permeable (the larger o;) barrier a will oppose AP 
and retard volume flow across the entire system. Alternatively, if o;? >0/*, 
the sudden application of a hydrostatic pressure to compartment i will 
deplete solute from the central compartment, since the solute content of the 
fluid stream crossing membrane 6 will be less than that contained by the 
stream crossing membrane a. Thus, when a steady state is achieved, c/'”" < c 
and the effective osmotic pressure across membrane 5 will oppose volume 
flow in response to the externally applied hydrostatic pressure. In short, if 
o;7  o;”, there will be an accumulation or a depletion of solute in the central 
compartment that, in either instance, will lead to the development of an “‘in- 
tramembranous effective osmotic pressure” that reduces the effect of AP on 
J,. Thus, the observed overall resistance to volume flow will always exceed 
the sum of the two elemental resistances. Further, as discussed by Kedem 
and Katchalsky (1963c) the degree of solute depletion or accumulation in the 
central compartment is a function of the direction and magnitude of J, so 
that, in general, a plot of J, versus AP will be nonlinear and will display 
rectification (i.e., for a given pressure difference the resistance to flow will be 
greater in one direction than the other); this point is illustrated by House 
(1974, p. 74). 


4.6. Summary 


Although the irreversible (or nonequilibrium) thermodynamic ap- 
proach does not enjoy the unassailable status of classical thermodynamics, it 
has proved to be remarkably successful for the description of interactions 
among flows and forces in relatively well-defined membrane systems. The 
applicability of the linear phenomenological equations and Onsager’s relations 
have been repeatedly demonstrated for the cases of solute-solvent and solute- 
solute interactions even in systems that are significantly displaced from equi- 
librium. The descriptions are phenomenologic and do not rest on assumptions 
regarding mechanism. Thus, although the approach can provide an internally 
consistent and experimentally useful set of descriptive parameters, it sheds 
little insight into underlying mechanism. As noted, some efforts have been 
made to translate the phenomenological coefficients into more mechanisti- 
cally meaningful terms, but it is not always clear whether these represent real 
advances or whether they simply substitute one set of phenomenologic state 
ments with another. For this reason, the approach has been subjected to some 
mixtures of skepticism and criticism. Clearly, a great deal more must be 
learned about the microscopic properties of solvent, solutes, and membranes 
before realistic, mechanistic models of transport processes can be formulated; 
but, until then, a useful and internally consistent phenomenologic treatment 
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ds better than none at all. Certainly, the contributions of classical thermo- 
fynamics, which rank among the most noble achievements of the human 
Intellect, cannot be criticized as being devoid of meaning simply because 
many antedated our current concepts regarding the molecular structure of 
matter. 
Finally, there are a number of important practical messages to be gleaned 
from the phenomenologic treatment of composite membranes. 
First, we see that asymmetric behavior (e.g., rectification of solute perme- 
ability coefficients, reflection coefficients, and pressure-volume flow rela- 
tions) which, on intuitive grounds, is not generally expected of an inert 
barrier can arise simply from the presence of heterogeneities along the 
eadth or thickness of the barrier. Clearly, an inert membrane penetrated by 
form pores having different radii and/or by nonuniform pores that are 
Petrow at one end and wider at the other would be expected to display such 
behavior and no more complex explanations need be invoked. 
Second, it is clear from these simple considerations that overall straight (L 
and «,;) and cross coefficients (0;) need not bear a clear-cut relation to the 
fespective elemental coefficients (i.e., they need not simply be area-weighted 
averages). Thus @équivalen “/ pore dimensions of complex membranes cal- 
culated from values of Lp using Poiseuille’s law or from values of 0; using the 
Renkin—Durbin approach (see Section 4.3) may be entirely unrealistic. 
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5 Active transport 
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5.1. A definition of active transport 


As discussed briefly in Chapter 1, in a system characterized by a 
single force and a single flow, the flow must be directed from a region of 
higher electrochemical potential to one of lower electrochemi¢al potential so 
that T d;S/dt =J;X;>0. In other words the flow and the force must have 
the same sign. However, if there are two or more flows and forces, negative 
coupling between flows can drive a substance in a direction opposite to its 
conjugate force; we referred to this phenomenon as energy conversion (see 
Section 1.6). Such phenomena are, of course, characteristic of biological 
membranes, but can also be observed in artificial membrane systems. For 
example, as discussed in Section 4.3, if a membrane is more selective for 
solute than solvent (0; <0), solvent flow will be directed from a region of 
high solute concentration to the region of lower solute concentration, a 
direction opposite to the osmotic driving force (anomalous osmosis). Further, 
from equation (4.25) we see that J; can be negative if J, is negative so that 
interactions between the flows of i and solvent can drive i from the compart. 
ment of lower concentration to that of higher concentration (e.g., solvent 
drag). 

In 1948 Rosenberg proposed that the term active transport should be 
applied to flows that are directed against their conjugate driving forces, 
that is, J; is “active” if J;X;<0. Conversely, if J;X;>0, J; is considered 
“passive.” This is an unambiguous definition based entirely on classical 
thermodynamic grounds that can be readily applied experimentally. 

An alternate definition was proposed by Kedem (1961), who simply re. 
arranged equation (1.35) as follows: 


Ji=(Xi/Ri) - LL Ri Fj/Rit) - Rr F/Ris) (5.1) 
i 


where J, is the flow of a chemical reaction and R;, is the cross coefficient 
that quantitates the coupling between J; and J,. The J; represents all other 
flows. Thus, in an isothermal system J; is determined by its conjugate driving 
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force, interactions with the flows of solvent and/or other solutes (“‘osmo- 
‘@imotic coupling”) and coupling to the flow of a chemical reaction (“chemi- 
Osmotic coupling’). Kedem suggested that the term “‘active transport” be re- 
‘served for those flows of i where R;, # 0. That is, J; is “active” if it is coupled 
p a chemical reaction; if not, it is referred to as “passive.” 

The definition proposed by Kedem has a number of advantages over that 
proposed by Rosenberg (Curran and Schultz; 1968%Schultz, 1 

First, it serves to distinguish between flows that are directed against their 
Wnjugate driving forces because Of interactions with the flows of other 
plutes or solvent and those that are driven by a chemical reaction. Some 
fer to the former as “secondary active (osmo-osmotic) transport” and to 
the latter as “primary active (chemi-osmotic) transport” (Mitchell 1967, 
970). Since a major purpose of any categorization of transport processes is 
0 direct future research, a definition that implies coupling to a biochemical 
faction is more useful than one based exclusively on the direction of the 
"W@nsport process. 
Second, the Kedem definition does not depend upon the direction of flow. 
e direction of J; may be the same as that of X; but it still will be con- 
dered active if it is accelerated or slowed by “positive” or “negative” cou- 
pling to J,. In contrast, the Rosenberg definition, which hinges on the direc- 
on of J;, can lead to an awkward situation. For example, in the absence of 
jupling to the flows of J;, equation (5.1) reduces to 


Ji = (X¢/Rii) - (RirJ/Ris) (5.2) 


Now, if Rj; is very large (e.g., the barrier is essentially impermeable to 7), 
® -(K;,J,/R;;) and is essentially independent of X;. Thus, in principle we 
pould adjust X; so that in one instance the flow of i is directed in the same 
ection as its conjugate force and in another the same flow is opposed to its 
iniving force. According to the Rosenberg definition, in the first instance J; 
auld be considered passive and in the other it is considered active in spite of 
ihe fact that J; is not influenced by X; and in both instances it is driven by 
ie Same mechanism! 
Thus, the Kedem definition departs from the classical thermodynamic 
sunds of the Rosenberg definition and introduces considerations of coupled 
dws and flow rates. However, although it is unambiguous and intuitively 
ye satisfying than the Rosenberg definition, it is far more difficult to apply 
Ni practice. In principle, in order to unequivocally establish positive or nega- 
' coupling between J; and a chemical (metabolic) reaction, either (a) the 
Wuction must be identified and a nonzero R;, must be demonstrated or (b) 
1 possible J; and Rj; must be determined in addition to R;; and X;; clearly, 
i complex biological system this could pose formidable experimental 
ybutacles. 
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a2) Experimental criteria for categorizing transport processes 


As discussed in Section 2.5, if the movement of a substance across 
a membrane is the result of simple diffusion and can be described by an inte- 
grated form of the Nemst-Planck equation, the relation between the two 
unidirectional fluxes is given by the Ussing flux-ratio equation 


TPF? = (ciPlez!) exp(-z; FAY/RT) (5.3) 


Since the advent of commercially available radioactive isotopes for many 
biologically important ions, this relation has been applied extensively and 
has provided the most important criterion for identifying the simplest of all 
transport processes. Thus, if the unidirectional fluxes conform to equation 
(5.3) over a wide range of values of c,°, c#, and Ay, it is safe to conclude 
that the flow of i is driven solely by differences in chemical and electrical 
potentials and traverses the membrane in its free form without interactions, 
complex formation, and so on (see Section 2.5 for the assumptions under. 
lying the flux-ratio equation). 

It is an all-too-common misapprehension that the flux-ratio equation can be 
employed to distinguish between “active” and “passive” transport perhaps, 
in part, because of the title of this landmark contribution, “The distinction 
by means of tracers between active transport and diffusion” (Ussing, 1949) 
It must be emphasized that failure to conform with equation (5.3) does not 
exclude the possibility that the flows are driven entirely by their conjugate 
driving forces and most certainly does not imply the presence of coupling 
between the flow of i and a chemical reaction. In short, equation (5.3) 
simply describes a sufficient but not necessary criterion for the conclusion 
that the movements of a given substance are attributable solely to differences 
in electrical and chemical potentials. 

It is instructive to consider, briefly, processes that can lead to deviations 
from the flux-ratio equation inasmuch as some of these departures can be 
distinctive. For the sake of simplicity we will consider an uncharged species 
or an ion when Ay is clamped to zero. Under these conditions equation (5,3) 
teduces to 


(5;°1F°) = (e:/e;') 


Single-file diffusion 


Suppose i crosses the membrane through narrow pores that can only 
accommodate a single row of two or more particles. Under these conditions 
tracer and parent species will not cross the barrier randomly. Insiead, the 
tracer movement will interact with nontracer movement, since it must pass 
through the pore by a series of jumps caused by the impact of collisions ut 
the two ends of the pore. 

It can be readily shown (Hodgkin and Keynes, 1955) that this will accelet 
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the movement of the tracer in the direction of its electrochemical poten- 
difference and retard the flow in the opposite direction. The flux-ratio 
juation takes the form 


Ji21F'° = [(eAlef) exp(-z; FAW/RT))" (5.4) 


here n is the number of particles in the single file. When z;=0 or Ay = 0, 
quation (5.4) reduces to (J;°!/J/?) = (c;°/e')". Clearly, equation (5.2) will 
inform with the flux-ratio equation (5.3) only when 


[(ei2/c;') exp (-z; FAW/RT)} = 1 


or J; = 0. (For other formal descriptions of single-file diffusion see Sjodin, 
965; Hladky, 1965; Heckmann, 1972; and Hille and Schwarz, 1978.) 


Membrane interactions 


If the movement of an uncharged solute across a membrane is the 
gult of complex formation with a mobile membrane component, or “‘car- 
jer,” but is not affected by any nonconjugate driving force, the process is 
ferred to as facilitated transfer. This “carrier-mediated” process will be 
ssed in greater detail in Section 6.2. Suffice it to say at this point that 
flux ratio is given by 


(71TH?) = [eK + ef NN le'(Ke + 67°) (5.5) 


here K, may be a complex parameter and is defined by equation (6.10). 
‘early, equation (5.5) will only conform to the Ussing flux-ratio equation 
inder the trivial condition when c,° = c;! and J; = 0. 

Finally, significant deviations from the Ussing flux-ratio equation can result 
if, in addition to simple diffusion, the molecule interacts with a membrane 
#Mponent that mediates an obligatory one-for-one exchange of the solute 
sross the barrier. This type of transfer process is termed exchange diffusion 
ind was first proposed by Levi and Ussing (1948). It will be discussed in 
wivater detail in Section 6.2; suffice it to say at this point that it closely 
embles facilitated transfer, except that it cannot give rise to a net flux but 
simply increments both unidirectional fluxes equally. If simple diffusion is 
jmplicated by exchange diffusion, the flux ratio for an uncharged solute is 

n by 


Ie 6 + SP) (5.6) 
Jf of +H 1P) 


Where J/” is the maximum velocity of exchange diffusion (see Section 6.2). 
learly, the numerator and denominator of the right-hand expression are 
incremented equally by a term that weights the magnitude of the diffusional 
process and the contribution from the exchange process. When J; >> P,, 
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then (J;°'/J;°') approaches unity; when P; >>J,, the contribution of the 
exchange process is small and (J;°'/J;'°) > (€;°/c;'). 


Interaction between tracer and parent species 


One of the assumptions underlying the derivation of the flux-ratio 
equation is that the flows of the tracer and the parent species do not interact 
(i.e., the cross coefficients are zero, so that each flow is driven exclusively 
by its conjugate driving force). In general, this may not be true, and Kedem 
and Essig (1965) and Essig (1966) have described the effects of positive and 
negative coupling between these flows on the flux ratio. An important con- 
sequence of their analysis is that these interactions can give rise to departures 
from the predicted ratio that resemble those attributed to single-file diffusion 
and exchange diffusion. These results caution against the inference of specific 
membrane properties or transport mechanisms on the basis of departures 
from the predicted flux ratios. 

It should be stressed that in all of the conditions described here, which can 
lead to significant deviations from the flux-ratio equation, J; is driven entirely 
by X; and thus must be considered “passive” by any criterion. Let us now 
consider situations where J; is influenced by nonconjugate forces. 


Solvent drag 


As we did in Section 2.5, let us assume that compartment o contains 
a tracer j and that c =0. According to equation (4.25), the net flux of / 
from compartment 0 to compartment i, which is equivalent to the unidirec. 
tional flux, is given by 


Jj°! = (¢;°/2) [2wjRT + (1 - 7) Jy] (5.7) 
Similarly , if compartment i contains a tracer k and c,° = 
Jie? = (Cx! |2) [200 RT - (1 ~ 0%) Jv) (5.8) 


Assuming that j and k are perfect tracers for i, the flux ratio is given by 
ee “i. 2 RT + (A= 01) Jo (539) 
J#?— \c#]) | 2a,;RT- (1 - 0) Jy ; 


Thus, unless J, = 0 or o; = 1, one unidirectional flux will be augmented and 
the other will be retarded by volume flow. The magnitude of this effec! 
clearly depends on the relative magnitudes of w;RT and (1 - 9;) Jy. 

A similar expression can be written for possible solute-solute interactions 
(solute drag), which takes the form 


Nig i. (<4) [a Se] (5.10) 
"Fa c;! 2w,RT - Lid; 
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where J; is the net flow of a species j and Lj; is a measure of the extent of 
coupling between J; and the flows of i. 


Active transport 


The effect of superimposing an active transport process on a simple 
diffusional process depends on the characteristics of the former. In the sim- 
plest case if i is actively transported from compartment o to compartment i 
{see Section 6.5) 


Jit = Prey? + [cI (K, + ¢(°)) (5.11) 
‘and J;'° = P;c;'. Thus, the flux ratio is given by 
yt (4) ha 
= 1 sa * 
a Nee P(K, + c;°) ihe 
Summary 


Figure 5.1 shows the relations between (J,°/J,!°) and (c;%c/) cal- 
‘ulated from equations (5.3)-(5.6), (5.9), and (5.12). The value of ¢;! was 
fixed at 100 mM and c,° was varied from 20 to 200 mM; other transport 
fameters were selected arbitrarily for scaling purposes. We see that the 
jations from the line of identity are in some instances quite distinctive. It 
particularly important to note that for any value of (c;%c;/) other than 
ity all of the systems we have considered yield flux ratios that could be 
| terpreted as indicative of active transport. Thus, flux ratios must be 
Wnterpreted with caution and can be particularly misleading when examined 
yer a narrow range of concentrations. On the other hand, data obtained over 
wide range of concentration ratios may provide clues to possible underlying 
sport processes. 


ie. An irreversible thermodynamic model for active transport 


The Kedem definition of active transport implies that the “flow of a 
hemical reaction” J, can be treated as is any other flow in a system and that 
Wupling between this flow and the flow of matter can be viewed as the 
dus operandi of an active transport process. Accordingly, in a system in 
Ich i is actively transported as a result of coupling toJ/,, we may write the 
lipation function 


T(d,S/dt) =J,;X;,+J,X, (5.13) 


here X, is the conjugate driving force for the flow J,. 
Now, what is the conjugate driving force for the flow of a chemical reac- 
tion? Consider the reaction 


Wg A tugpBt'**——~ vol + vpDt::: (5.14) 
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Figure 5.1. Relations between flux ratio and concentration ratio for 
an uncharged solute, determined as discussed in the text: a, simple 
diffusion ; b, single-file diffusion; c, facilitated transfer; d, simple 
diffusion plus solvent drag; e, simple diffusion plus exchange 
diffusion ; f, simple diffusion plus active transport. 
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where the v’s are stoichiometric coefficients of the reaction in which the 
feactants A, B, . . . are converted to the products C, D, . . . Assuming that 
only reactants are present at the beginning of the reaction, then, at any time 
thereafter the progress, extent, or degree of advancement of the reaction df 
is given by the number of moles of reactant that have disappeared and/or the 
number of moles of product that have appeared. Accordingly, ? 


df = -(dny |v4) = -(dng/vg) = (dnclvc) = (dnp/vp): * (5.15) 
Now, at constant T and P, the change in free energy of the reaction is 
(dG) 7, p = Ma dng + Up dng + uc dnc + Up dnp t+: (5.16) 
(dG)r,p = 2. aj dn; (5.17) 
j 


where dn, is negative if j is a reactant and positive if j is a product. Substitut- 
equation (5.15) into (5.17), we obtain 


(dG)r, p= dS Do vu (5.18) 
J 


d differentiating by time (t), we obtain 
(dG/dt)r, p = (dS/dt) 20 yu, (5.19) 
ii 


where, now, the stoichiometric coefficients are negative for reactants and 
tive for products. Clearly, df/dt is the rate of advancement of the reac- 
or the flow J,. Thus, equation (5.19) can be rewritten 


(dGlat)p, p=Jp- 2 vay (5.20) 
J 


As discussed in Section 1.3, the free energy of a system decreases during 
y spontaneous process (i.e., dG <0). Further, as illustrated by equation 
1,12), in a system at constant T and P that does not perform “useful work” 
dG/dt) = T(d;S/dt). Thus, equation (5.20) reduces to 


T d,S/dt=J, >) (-vju;) (5.21) 
i 


it follows that the conjugate driving force of the reaction is given by Z(-»,u;) 
h (at constant T and P) is the decrease in Gibbs free energy of the reac- 
per unit rate of advancement. This relation was formulated in 1936 by 
pioneering work of DeDonder, who defined the driving force for the flow 
a chemical reaction as the affinity A. 

Assuming that flows and forces are linearly related, we can describe this 
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two-force, two-flow system by the following equations (see equations 1.35): 


X,= Any = RJ; + Rir Jy 
and (5.22) 
A=RyJ, + RrrJ, 


where the R’s are generalized resistances and, assuming that the reciprocal 
relations apply, R;, = R,;- 

As discussed in Section 1.6, in principle, the linear phenomenological equa- 
tions and Onsager’s relations are expected to apply only when a system is not 
displaced too far from equilibrium, and in the final analysis the validity of 
these relations must be established empirically. Whether biological transport 
processes can be described by means of linear relations among J;,J,, Aji;, and 
A is an unsettled question. Katchalsky and Curran (1965) have argued, on the 
basis of kinetic analysis, that a linear relation between J, and A is likely to 
hold only over a very limited range where the reaction is very close to equi 
librium. Prigogine (1961) has also argued that linearity is to be expected only 
when A <RT (approximately 600 cal/mole), a value that is considerably 
smaller than the free energy changes associated with most exergonic bi 
chemical reactions. More recently, Prigogine and Lefever (1975) have shown 
that if an overall chemical reaction is the result of a number of intermediate 
reactions, each of which is not very far displaced from equilibrium, the over 
all reaction rate may be a linear function of the overall affinity even though 
the latter may be much larger than RT. 

Experimental data bearing on this issue are limited. Blumenthal, Caplan, 
and Kedem (1967) have shown that the linear phenomenological equations 
and Onsager’s reciprocal relations can be used to describe coupling between 4 
chemical reaction and current flow (ion transport) in an artificial membrane 
system where A was as high as 3000 cal/mole. It is of interest that in thel 
system (as in most enzymatically catalyzed reactions) J, was a decidedly 
nonlinear function of A; nevertheless, the linear formalism yielded the cori! 
data. j 

In short, the applicability of equations (5.22) to the description of hh» 
logical transport processes is a moot point that cannot be settled on the busy 
of the available data. Nonetheless, it is instructive to carry through an analysis 
of these equations inasmuch as the results have bearing on energy conver |eyy 
in general and may provide some insight into biological transport process 
The approach we will employ closely follows that of Kedem and Caplay 
(1965) and Essig and Caplan (1968). 

Rearranging equation (5.22), we obtain 


is) [Are 
fin (OM) Ute (5,28) 
: = Ry 
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We will focus mainly on the situation where coupling to J, retards the flow of 
J; in the direction of its conjugate driving force or actually drives J; in the 
Opposite direction so that J; Au; <0. This is referred to as negative coupling 
dif we define the directions of J, and Aji; as positive, then R,, > 0. 

We start by defining two combinations of the phenomenological coefficients 


= fee we = “Rir 
2 (Fe) at a Ga roe 
Since R;,? <RjR,,, it follows that -1<q<1; q is referred to as the 
e of coupling for the following reasons: 
(a) If q=+1, then (J;/J,) = +Z. Thus, the ratio of the flows is independent 
the ratio of the forces, and if one flow is fixed the other is uniquely deter- 
ed; in other words, the two flows are completely (tightly) coupled. It 
ows that if the driving flow J, is fixed then J; is independent of Ag; and 
system operates as a “constant current source.” This is, of course, an ideal 
treme. Alternatively, if J; is affected by varying Ag,, then J, would be 
ected proportionally. 
{b) If q = 0, then R;, = 0, and 


(34) 27 hn; _ (AM IRu) _ Li Adi (5.25) 
J, A (A/R;y) L,,A : 


, when q=0, the ratio of flows is proportional to the ratio of forces; 
flow is influenced only by its conjugate driving force so that the flows 
entirely uncoupled. 
©) For values of q between 0 and -1 (negative coupling) the relation be- 
n the flow ratio and the force ratio varies as illustrated in Figure 5.2. 
© quadrant delineated by the origin (0), the value (J;/J,Z) = -1, and the 
(Afi;Z/A) = +1 is referred to as the driving region since in this region 
, > Au;. Thus, J; is driven in a direction opposite to Afi; (i.e., J; Al; <0). 
ily, if we define J,A > 0 and Aji; as positive, J; is negative throughout 
tégion.) 

quadrant immediately above this one, delineated by the origin, the 
(J;/J,Z) = 1, and the value (Agi;Z/A)=1, may be referred to as the 
ing region since in this region J; Aux; > 0 but J; decreases with increas- 
#/ or decreasing force ratio (i.e., when A increases at constant Ag; or when 
(lecreases at constant A). 

the interface between the retarding region and the driving region J; = 0. 
Aji; = 0 only when q = 0; for all nonzero values of q, Afi; > 0. Thus, 
my to an exergonic chemical reaction establishes and maintains an 
ochemical potential difference for i in the absence of a net flow of i; 
eater the degree of coupling, the greater will be Aji; at fixed A. These 
ts of intersection with the abscissa are referred to as static heads. The 
Hiditure of metabolic energy to maintain differences in electrochemical 
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Jj/InZ to a source of energy, but in neither instance is useful, or osmotic, work 
performed; thus, the efficiency of energy conversion [see equation (1.40)] is 


zero. Between these boundaries chemical energy is converted into osmotic 
work and the efficiency of the process is given by 


RETARDING eS 
REGION n= -(J;Mi;/J,A) (5.26) 
Substituting equations (5.24) into equations (5.22), we obtain 
Au; + (q/Z) A 
Pe ee ae 5.27 
Ry(1 - 4’) 27) 
(q/Z) Sui; + (1/Z)" A 
ee oe 5.28 
‘ Ry(l - ¢’) G29) 


From these expressions it follows that the static head electrochemical 
jotential difference is given by 


(Adi;)o = (Ad;)J,=0 = ~qA/Z =, - qA(Rii/R,,)"!? (5.29) 


‘ 
DRIVING d the rate of energy expenditure (power) needed to maintain this static 
REGION dis : 
(Aiii)o” (1 
SA) s.<9 = 5 lice = ; 
Gr Alyn = Re (aro (5.30) 


‘Thus, the rate of energy investment needed to maintain a given static head 
ases if the passive resistance to the flow of i decreases and/or the degree 
eoupling decreases. 

‘The level flow is given by 


ee A(q/Z) 
ay ee Oe een eae ! 
(Ji) ap;=0 e na 25 (5.31) 
the rate of energy expenditure needed to sustain this flow is 


(J,A)apj=0 = Ri (JiJan,-0/9? (5.32) 
, the rate of energy investment needed to maintain a given level flow 
ases if the resistance to flow increases and/or the degree of coupling 
ases, 
the efficiency of energy conversion is zero at each of the two bound- 
of the driving region, it follows that n must pass through a maximum 
en these boundaries. It can be readily shown that this maximum effi- 
is solely a function of the degree of coupling and is given by 


max = q7/[1 + (1 - q?)"/?}? (5.33) 


+ Nmax = 1 when q = +1 but falls very quickly with decreasing values of | 
‘Wor example, when q= 0.95, nmax = 0.52, so that extremely tight cou- 
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Figure 5.2. Relations between flow ratio and force ratio for differ 
ent degrees of negative coupling (modified fram Kedem and Caplan, 


1965). 


potentials in the absence of net flows is characteristic of all symmetrical cells 
under steady-state conditions (i.e., when cell composition is constant). 

In addition, when Au; = 0, J; = 0 only when q = 0; for all nonzero values «if 
q, there is a finite flow of J; in the direction opposite to J, that increases will) 
increasing q. A net flow of i when Ay; =0 is referred to as level flow, \nu 
much as energy is invested to bring about a flow of matter in the absence f 
an electrochemical potential difference and, thus, without the performan 
of osmotic work. A number of epithelial tissues such as gallbladder, tony! 
proximal tubule, and small intestine invest metabolic energy in ordey {) 
absorb Na and Cl against minimal electrochemical potential differences wn 
operate close to level flow. 

Thus, the two boundaries of the driving region are the conditions, /, * 
when Aj; #0 and J; #0 when Aji; = 0. Both are dependent upon coupling 
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pling would seem to be a requisite property of efficient energy conversion 
processes. 

Finally, these relations suggest experimental approaches for determining A 
if it is constant. From equations (5.27), (5.28), and (5.29) it can be shown 
that 


A = -Z(Ati;)o/q = (8J;/0J,) (Alii) (5.34) 


Thus, determination of the relation between J; and J, under conditions close 
to static head provides a measure of A. Alternatively, from equations (5.27) 
and (5,28), it can be shown that 


A =~ [0(Au;)/8J,] i FI))az,-0 (5.35) 


Thus, determination of the relation between J, and Agi; under conditions 
close to level flow provides another measure of A. It should be noted that 
the determination of A using equations (5.34) or (5.35) involve independent 
measurements. Thus, obtaining the same value of A using both methods 
Suggests that A is constant throughout the driving region and supports the 
internal consistancy of the model. 

In concluding this section, it should be reemphasized that the model pro- 
posed by Essig and Caplan rests on several simplifying assumptions that are 
Open to question. These investigators have summarized the body of experi- 
mental data, largely drawn from studies on isolated frog skin and toad urinary 
bladder, that are consistent with this linear phenomenologic approach and 
have derived values for A in these systems that appear to be quite reasonable 
(cf. Caplan and Essig, 1977). However, further study is necessary to sub- 
Stantiate this simple model and to establish its range of applicability. In any 
event the approach employed by Kedem, Caplan, and Essig provides a formal 
foundation for the appreciation of energy conversion resulting from the inter- 
actions of flows and forces. 


6  Carrier-mediated transport 


6.1. A definition of the term carrier 


The notion that biological membranes possess components capable 
of interacting with a solute at one side, conveying it across the barrier and 
teleasing it at the other side, was introduced in order to account for a number 
of characteristics of many biological transport processes that seemingly could 
not be attributed to simple diffusion. These characteristics include structural 
Specificity (e.g., chemical and stereoisomeric specificity), saturation kinetics, 
and competitive inhibition.! However, as discussed in Section 2.6, diffusion 
through pores or channels that contain specific binding sites can also exhibit 
many of these characteristics. 

The essential difference between a pore and a carrier, as these terms are 
generally employed, is that a pore has openings on both sides of the mem- 
brane so that molecules in the two surrounding solutions have simultaneous 
access to the transport pathway. In contrast, the binding site of a carrier is 
accessible from only one side of the membrane at a time and alternates be- 
tween the two sides. Thus, a carrier necessarily involves a mobile or transi- 
tional element, whereas all of the structural properties of a pore can, in 
principle, be stationary or fixed. The precise nature of this “mobile” or 
“transitional” property is unknown; a variety of models that have been sug- 
gested are described in the excellent review by LeFevre (1975). In any event, 
the restriction that solutes at both sides of the membrane do not have simut- 
taneous access to the binding site is the sine qua non of the carrier hypothesis 
that distinguishes this transport mechanism from specific pores or channels. 

The kinetic behavior of carriers is generally described using schemes such as 
that illustrated in Figure 6.1. It is assumed that the carrier X is (a) confined 
to the membrane; (b) can exist in the free or bound (iX) forms; and (c) that 
the total amount of carrier per unit area is constant. In addition, it is fre- 
quently assumed that the association-dissociation reactions occur only at the 
membrane interfaces.” Carriers facing the outer compartment are designated 
X° or iX° and those facing the inner compartment are designated X? or iX!. 
The k, and k, are association rate constants for the formation of the com- 
plex iX; k_, and k_, are the corresponding dissociation rate constants. P, and 
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Figure 6.1. General model of carrier-mediated t i i 
1 : ~ tansport involving » 
carrier XY that can form a binary complex with the solute i. 


P, are the rate constants for the translocation of iX and X across the barrie, 
translocation is assumed to be a first-order process, but its mechanism noed 
not be specified. 

Although this model can be solved generally using the familiar approach of 
chemical kinetics, for the sake of simplicity we will assume that translocation 
through the membrane is slow compared to the association-dissociation reac 
tions at the interfaces (i.e., the translocation Step is rate limiting). Unde; 


these conditions the reactions at the interfaces are close to equilibrium ani 
can be described by 


Ky = ¢,° [X°] /[iX?] (6.1) 
and 


Ky = ef (X4] /[ix"] (6.2) 
where K, and K, are dissociation constants given by (k_,/k,) and (k./k,), 


respectively, and the terms in brackets are concentrations of free or bound 
Carrier per unit area of membrane. 


The assumption that carrier is confined to the membrane and that the 
amount is fixed is given by the relation 


[Xe] = [X?] + [X?] + [x4] + x4] (6.4) 


where [X;] is the total concentration of all forms of X and is constant, 
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Y, under steady-state conditions, the distribution of X at the two 
aces is constant (time independent) so that 


P,([iX°] - [iX"]) = P(X") - (X71) (6.4) 


the rate of flow of bound carrier from o to i must be equal to the rate 
tm of free carrier from? to 0). 

four independent equations permit us to solve for the concentrations 
four possible carrier forms (X°, iX°, X', and iX*) in terms of ¢,°, ¢}, 
| 2 P 1, and P ae 

¢ unidirectional flux of i from compartment o to compartment / is then 
ly J,°! = P,[iX°]; the unidirectional flux in the opposite direction is 
© P, (iX"]; and the next flux is given by J; = J;2! - J,9. 

will now consider some of the classical categories of carrier-mediated 
ort in terms of this kinetic approach. We will first consider the trans- 
of uncharged solutes and later deal with charged species. 


Facilitated transfer 


The term facilitated transfer was introduced to describe carrier- 
ted transport processes that are not directly coupled to nonconjugate 
of energy and, thus, cannot lead to transport against an electrochemi- 
potential difference; they can only lead to equilibration of the transported 
f@ across the membrane and are, therefore, sometimes referred to as 
tated diffusion. These processes can be described by the model illus- 
d in Figure 6.1 and Equations (6.1) to (6.4) when K, = K, =X. Under 
condition it can readily be shown that 


P, (X}e;°P2K + Pyc;') 


0 ee ee nahi‘) 

Cj (P,K+P,K+ 2P,¢; )+ K(2P,K + P2¢; + Pye; ) 

jo Pi [X,]c/(P2K + P,c;°) 

J; Fy o o oO (6.6) 

Cj (P,K +P,K + 2P2C; )+K(2P2K +Pic; + P2¢; ) 

ce the denominators in equations (6.5) and (6.6) are identical, 

Ore tn! 

Blo = J,= [X}(? - ¢7') PiP2K (6.7) 


¢)°(P2K + P,K + 2P,c/) + K(2P2K + Pyc} + Pyc}) 


Thus, when ¢,° =c;', J;=0. As will be shown (Section 6.5), this limiting 
dition is entirely because K,=K, so that the carrier process is sym- 

nical. 

‘According to equation (6.5), when ¢;! is constant, J;°! will be a saturable 

tion of ¢,° given by 


Jyh SMe; (Ky + 01°) (6.8) 
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where J/” is the maximum velocity and is given by 


E P,[X,] @2K +Py¢;') 
PK +P,K + DP Cy 


and K; is the value of ¢,? at which J,?! = J/"/2, given by 


_ K(2P,K + Pic + Pyc/') 
P,K +P,K + 2P,¢/ 


When P, = P, =P, equation (6.9) reduces to J/” = P[X;] /2; equation (6.10) 
reduces to K, = K; and the unidirectional fluxes are given by 


(ORE) ae 
K+ ci? 


Because of the similarity between equation (6.8) and the familiar Michaelis» 
Menten formulation of enzyme kinetics, the K,’s determined experimentally 
from plots of J,°! versus c;° or J;!° versus c;' are frequently identified with 
the equilibrium constant of the interaction between the transported solutt 
and the carrier (K). It should be clear from equation (6.10) that this is true 
only when P, =P,; when P, #P,, the experimentally determined K, is no! 
equal to K, but is a rather complex function of K, the P’s, and either c,° wy 
Cr. 
Thus, before an experimentally determined K, can be interpreted as a mew 
sure of the association or dissociation constant of the solute-carrier interay, 
tion, we must establish the identity between P, and P3. 

How can this identity be established? 

It is clear from a comparison of equations (6.5) and (6.11) that when 
P, =P, J;°! is dependent on c;° but is not influenced by c/; similarly, J;!? is 
independent of ¢,°. In contrast, when P, #P2, J;7! is influenced by c;! ani 
a” is influenced by c,°; these influences are referred to as transeffects lw 
cause the unidirectional flux is affected by the concentration of solute in the 
(trans) compartment toward which it is directed as well as the (cis) compus} 
ment from which it originates. 

When P, >> P,, equation (6.5) reduces to 


J,°1 = Pi [Xi] ei" 
. 7? [2 + (K/e#)] +K 


thus, when ¢;° is constant, J;°! increases with increasing c;'. This phenome 
non is referred to as transstimulation, inasmuch as increasing the concent!» 
tion in compartment / stimulates the unidirectional flux from compartment \) 

Further, from equation (6.7) we see that J;=0 when P, =0. Under this 
limiting condition the bound carrier iX simply shuttles back and forth acim 
the membrane and mediates the obligatory one-for-one exchange of solute ij 


J” (6.9) 


K; (6.10) 


cf, ae 
Kt c} 


Jet= and Jj? = (6.11) 


(6.12) 
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partment o for solute in compartment i without net transport; this pro- 
is referred to as exchange diffusion (see page 85). 

When P, >> P, equation (6.5) reduces to 


Jel = Py [Xe] ci? 
t 


= 6.13 
QR +e? +e} ( ) 


s, when ¢;° is constant, J;°! decreases with increasing c}. This phenome- 
is referred to as transinhibition and occurs when the translocation of 
nd carrier is rate limiting so that an increase in ¢;! leading to an increase 
iX' reduces the amount of free carrier available for recycling across the 


brane. 
Thus, in principle, experiments can be designed to determine whether an 
ilibrating system exhibits transeffects and the absence of such effects 
t be demonstrated before we can consider K, a measure of the “affinity” — 
the carrier for the solute. 


Competitive interactions and countertransport 


Consider the system illustrated in Figure 6.2, where a carrier is ca- 
of interacting with two uncharged solutes i andj, with K, being the dis- 
tion constant of iX and K, the dissociation constant of jX. The six equa- 
that describe this system are 


¢;° [X°] - Ki[iX°] =0 
ej? [X°] - K2[iX°] =0 
cf [X4] - K, [iX#] =0 
oj [X"] - K,[jX‘] = 0 
[X?] + [4] + [Xe] + fax!) + 7X?) + LX] = 1%] 


(6.14) 


P,{iX°] - P, [iX!] + P,[X?] - P2[X"] + PaliX°] - PaliX"] = 0 
With some patience, these six equations can be solved for the six possible 
rt forms in terms of P’s, K’s and concentrations, and expressions for the 
ectional fluxes of i and j can be derived as above. 
an be readily shown that when P, =P, = P3 =P 
s Bf ia ¢ in 

K,[1 + (¢/7K2)] + 1° 
is precisely analogous to the expression for classical competitive in- 
Yon in enzyme kinetics. Clearly, J;°! decreases with increasing c,° and 


dP! (6.15) 
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Figure 6.2. Model of competitive inhibition in a system capable 
only of facilitated transfer. 


when c;°=K, the concentration of i needed to achieve a half-maximal 
velocity is 2K, (i.e., it has been doubled). 

Let us now consider the case in which c;° =0 but c} > 0. From equation 
(6.15) it follows that 


J, re c;? 


ote 
ue K, +¢/? (6.14) 
and 
mai 
Jf? = ae ee (6.17) 
K,[1 + (c//K2)] +¢; 
Thus, 
oe i 
Jie! - Ff = J, = 5" tes ee 
i i i i K, +¢,° Ki {i + (c)/K2)] +o! (6.18) 
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Ne now inquire: When will J; cease? By setting equation (6.18) equal to zero 
d rearranging, it can be shown that J; = 0 when 


(ci/c;°) = 1 + (¢/K2) 


Thus, in spite of the fact that this is an equilibrating system, it is capable 
f transporting i from compartment o to compartment i against a concentra- 
ion difference as long as c} > 0. In general, when c;? # 0, J; ceases when 


(ci /c;°) = (Ko +e )(K2 +¢;°) (6.20) 


‘that the system is capable of transporting i against a chemical potential 
ference, providing c} ey. 

s phenomenon, which was described formally by Wilbrandt and Rosen- 
(1957, 1961) is referred to as countertransport. 

The mechanism underlying countertransport is intuitively obvious! In the 
iple case where c;° = 0, J,°! is unimpeded, whereas J;!° is retarded by com- 
ition between i and j in compartment 7 for free carrier. This emerges 
patly by substituting the relation K, = cj [X‘]/[/X‘] into equation (6.19), 
gives 


(cf/e;°) = 1 + (UX'1/[X"]) 

se the total amount of carrier is fixed, as the ratio of [jX“] to [X*] in- 
eases, the amount of free carrier available to mediate the movement of i 
im compartment i to compartment o decreases so that i becomes “trapped” 
compartment i. Thus, when cj! >c;?, J;°! will equal J;° only when 
> ¢,? : 

jowever, it should be clear from our previous considerations that coupling 
4 Nonconjugate source of energy is required for any mechanism to trans- 
matter against an electrochemical potential difference; what is the source 
‘this requisite energy input in the case of countertransport? This issue will 
dealt with more completely later (Section 6.4); suffice it to say at this 
int that the energy for the flow of i against its concentration difference is 
ed from coupling to the flow of j from a higher to a lower concentration 
the coupling mechanism is the common carrier system. In short, be- 
j and i share a common carrier, the downhill flow of 7 from compairt- 
i 7 to compartment o can provide the energy for the uphill flow of i in 
Opposite direction. 

iis phenomenon can be illustrated more vividly by considering the con- 
om where P, = 0, Clearly, then, in the absence of j, J; = 0 and the system 
capable of mediating exchange diffusion. When / is present in one or 
h solutions, in general, there will be a net flow of i, which under steady- 
conditions must equal the net flow of j in the opposite direction; in 
, this system will carry out an obligatory one-for-one exchange of i and 


(6.19) 
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j. When P, = P3 =P and P, = 0, the solutions of equations (6.14) are 


Jie! = P[X,] Kr0i°(K.¢j + Ki¢)/D (6.21) 


and 


Jil? = PLX,] K20j\(K20;° + K1cj°)/D (6.22) 


where 
D= 2KiK2¢,°c;! + 2K17¢;°c/ +K,K2.0;2c} + K,K2? (cj + ¢,°) 
+ Ky? K2(¢/? + ¢/) + 2K27¢;2c (6.23) 


af and gfe are given by equations (6.21) and (6.22) after replacing c;° by 
éf", cj! by a. and K, by K,; the denominator remains the same. 
It follows from equations (6.21) and (6.22) that J; = 0 when 


(c;{/ei?) = (c/¢;°) (6.24) 


and under these conditions J; is also zero. This is the static head state for an 
obligatory, one-for-one, countertransport system. Since Ap; = RT In(c;'/c;") 
and Ay;=RT In(e}/c;?), it is obvious that the condition of zero flows is 
reached when Ay; = Ayj. If (¢//¢;°) > (ci/ei), there will be a flow of / 
from compartment o to compartment i against a concentration difference 
driven by an equal flow of j in the opposite direction from a higher to » 
lower concentration. Conversely, if (c¢;!/¢;?) > (c//e¢;?), there will be a flow of 
j from o to i against a concentration difference driven by an oppositely 
directed flow of i. 

It is instructive to compare equation (6.20), which applies when P, ¥ 0), 
with equation (6.24), which applies when P, = 0. Clearly, for any finite valuy 
of K, when ¢;! > cj? 

(Ky + ¢j/)(Ko + 6°) < C//e;°) 
Thus, when P, # 0, the static head value of Ay; will be less than Ay;. That 
is, because X can cross the barrier in the free form, the flow of i is not tightly 
coupled to the flow of j and the efficiency of energy conversion is less thy 
the ideal of unity. Thus, a high efficiency of energy conversion in a countey 
transport system requires that the free carrier be immobile. 

Finally, there is an obvious point that, nonetheless, should be emphasizi 
In the “real world” where the efficiency of energy conversion is less thy 
unity, some input of metabolic energy is necessary in order to maintaly 
steady-state differences in chemical potential; the countertransport system 
illustrated in Figure 6.2 cannot accomplish this by itself! For example, let ws 
suppose that at the outset of an experiment, c;? =c;! and c/ >c/?. Then, 
there will be a flow of i from compartment 0 to compartment / so that we 
will achieve the condition (c;!/c;°) > 1. However, in the real world the cow 
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dition c > c;° cannot be sustained indefinitely unless some other mechanism 
tervenes; otherwise, this concentration difference will ultimately vanish, 
however slowly. Consequently, the final steady-state condition of this iso- 
lated system is (c;//c;°) = (c;/c;°) = 1. The only way a true steady state can be 
stained, where (c;/c;°) and (c//c;°) are greater than unity, is if there is 
me other mechanism capable of replenishing j in compartment i. For 
example, j may be the product of a chemical reaction in compartment / or 
there could be a mechanism capable of actively transporting j from compart- 
Ment o to compartment 7. In any event, suffice it to say that a constant 
pply of energy is necessary to maintain any system in a steady state that is 
splaced from equilibrium. 


4. —_ Cotransport 


Countertransport (sometimes referred to as “antiport”) is one 
néchanism for carrier-mediated coupling between flows of solutes that can 
ud to transport against a concentration gradient. Another mechanism that 
prevalent in biological systems (cf. Schultz and Curran, 1970) is cotrans- 
(or “‘symport”) where the net movements of two solutes take place in 
same direction. Consider the general scheme illustrated in Figure 6.3, 
yhere i and j can combine with the same carrier X to form the binary com- 
es iX and jX and the terniary complex #/X. In the most general case, 
order of binding is random and all forms of the carrier can translocate 
ss the barrier. We can, of course, impose a preferred order of binding by 
sting the values of the K’s and can restrict the carrier forms capable of 
pinslocation by setting one or more of the P’s equal to zero. The resulting 

pressions for the fluxes of i and j can become cumbersome and their 
ivations are rather tedious, but the underlying approach is the same as that 
Witlined previously for simpler systems. The general solutions of the model 
1a number of selected conditions are given by Schultz and Curran (1970). 
¥ order to illustrate the principal features of cotransport, let us consider 
vo simple examples. 

Example I. Assume that i must bind to X first to form the binary 
plex iX and that only then can j be bound to form ijX and that the 
ry complex dissociates in the reverse order (i.e., j leaves first); this 
tion is achieved by setting the dissociation constant K4 =© so that 
» 0. The system then reduces to that enclosed by the dashed rectangle in 
re 6.3. Let us further assume that P, =P, = P3 =P. 

‘The solution of this model for the bidirectional fluxes of i and j are 


J, Pi ;? 


a! hd 6.25 
- [K, K2/(K2 + ¢;°)] + ¢/° : 
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and 
ipo a. (6,28) 


1 "KK, +¢{(Ki +c!) Ka +e; 


where J;” = [X],P. (It should be noted that there are no transeffects he. 


cause the P’s are equal.) 
Thus, the maximum velocity of i in either direction is independent of / but 


the concentration needed to achieve a half-maximal velocity is a function of /, 
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clearly, when c;=0, K,= XK, as is to be expected. On the other hand, the 
maximum velocity of 7 in either direction is a direct function of the con- 
centration of / in the cis solution. Further, from equations (6.27) and (6.28) 
we see that 


Jj ot te Jj? c 


—_— = — d —_— = 
Jet (Kr + ¢/?) Sr ecg ee +¢/) win 


The expressions on the right-hand sides of these equations are referred to as 
the coupling coefficients between the flows of i andj. 

It follows from equations (6.25) and (6.26) that J;=0 when (ci'/c:?) = 
(K, + of MK ain ¢/'); the resemblance between this expression and equation 
(6.20) for the countertransport system should be noted. Thus, (¢,!/c;°) > 1 
when c;° > c ; that is, the system is capable of transporting i against a con- 
centration difference and the required energy is derived from carriey-mediated 
Oupling to the flow of j in the same direction from a higher to a lower 
centration. 

From equations (6.27) and (6.28) it can be shown that J; = 0 when 


of _ (Ki +7) 

Gj? eK, +c”) 
us, (c;'/c;°) > 1 when c,? >c;, so that j can be transported against a con- 
tration difference energized by coupling to the flow of i down a concen- 
jon difference. 

Example II. Now let us consider the same sequence of binding and 
jation but with the restriction that P, = 0 so only the free carrier and 
terniary complex can translocate. Under these conditions the bidirec- 

fluxes of i and j are complicated by transconcentration effects; the 
it solutions can be found in the review by Schultz and Curran (1970). 
important point is that under these conditions the coupling coefficient 
unity so that J;°!=J,!, J/° =J/° and J,=J;. Further, J; and J; will both 
when (c;'/c;°) = (¢;°/c}); again, this expression should be compared 
equation (6.24) for the analogous countertransport system. 
Since for any finite value of K, when c;° >c/! 


[(K2 + ¢/?)(K, +¢/)] < (ccf) 


follows that the static head value of (¢;/c;°) will be less when P, #0 than 
P,=0. 

again, when P, = 0, the flows of i and j are tightly coupleq and the 
derived from the flow of one of these species from a higher to a lower 
ntration can be completely utilized to drive the flow of the other 
from a lower to a higher concentration. When P, #0, the flows are 
tightly coupled and the efficiency of energy conversion is reduced. 


(6.30) 
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In short, the principles that underlie carrier-mediated cotransport processes 
are strictly analogous to those that govern the behavior of countertransport 
processes. In both instances the energy derived from the flow of one species 
down an electrochemical potential difference can be employed to propel the 
flow of the other species against an electrochemical potential difference; 
the coupling mechanism responsible for this energy conversion (or transduc- 
tion) is a common carrier system, and the efficiency of energy conversion is, 
as discussed in Section 5.3, determined by the tightness (or degree) of cou- 
pling between the flows.? 


6.5. “Primary” active transport 


All of the carrier models we have discussed to this point are inher- 
ently symmetrical. We have seen that in the cases of cotransport and counter- 
transport, a solute i can be propelled from a lower to a higher concentration 
when there is an asymmetric distribution of the “coupled” species j. We 
referred to these processes as secondary active transport because of solute- 
solute (or osmo-osmotic) coupling. In the case of primary active transport the 
carrier system is inherently asymmetric and this essential asymmetry arises 
from coupling the system to an exergonic chemical (metabolic) reaction. In 
order to illustrate this point, let us consider the hypothetical model shown 
in Figure 6.4.4 We assume that the total amount of free carrier can exist in 
either of two mobile forms; a form designated by X which (as before) can 
combine with the solute i, and a form designated by Y which is incapable of 
binding i. We further assume that the distributions of total free carrier be. 
tween X and Y at each interface are given by the equilibrium constants 
K, =((¥°]/[X°]) and K3 =([Y¥*]/[X4]). The solution of the equations that 
describe this system for J; is simply 


i 
=70i_7io=7;m cf? Biiddd| ieee Ute ae iv 
Sy al oe irae nee en aa el 


Thus, both unidirectional fluxes conform to Michaelis-Menten kinetics, 
where the K; for J;?! is K,(1 + Kz) and that for J? is K,(1 + K3). Clearly, 
if K, #K, the system is inherently asymmetric, inasmuch as the concentra: 
tion of i in the outer solution needed to achieve a half-maximal flow from 
compartment o to compartment 7 will differ from the concentration of j in 
the inner solution needed to achieve a half-maximal flow in the oppoalty 
direction. 

From equation,(6.31) we see that when J; = 0 


(cilci?) = (1 + K3)/( + Ka) 


Thus, when K; >K2, the system is capable of transporting 7 from a lowes 
concentration in compartment o toa higher concentration in compartment /, 
if K3 <K-,i will be actively transported in the opposite direction. 


(6.32) 
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Figure 6.4. Highly simplified model of a system capable of active 
transport. 


Although this conclusion follows directly from kinetic considerations, what 
bout the energetics? Some insight into the energetics of this simple, primary 
ive transport system can be gained by describing the static head condition 
g the formalism of irreversible thermodynamics. 

In general, the dissipation function for this system is given by 


T d,S/dt =Jix Xjx +JxXx t+ JyXy 


here, as before, J’s are flows, the X’s are the conjugate forces and the sub- 
Nipts designate the species. When J; = 0, the static head condition, J;y =0 
that this expression reduces to 


d,S 


= P([X°] - [X4]) RT In rae) +P([Y°] - [¥*]) RT In ( vl) 


(6.33) 
1 steady-state conditions when J; = 0, 
P([X°] - (X*])= -PCY*) - [¥']) 
that equation (6.33) reduces to 
Tas = - 07 . t areal 7 
re P([Y°?] - [Y‘]) RT in [Y?] (6.34) 
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Since [Y°] =K,[X°] and [Y‘] = K3[X"], equation (6.34) may be written 


T da;S 
dt 


K 
= -P((¥?] - [Y‘])RT In GF) (6.35) 
2 
Now let us assume, as illustrated in Figure 6.4, that a chemical reaction is 
responsible for converting X‘ to Y‘ and that Y‘ then translocates from the 
inner interface to the outer interface where it spontaneously re-forms X° 
Then, under steady-state conditions 


-J,=P((Y?] - [¥']) 


where J, (as before) is the rate of the chemical reaction. Further, it should be 
clear that the expression 


ere i) 


is simply the difference between the chemical potentials of the products, 
X° and Y‘, and the reactants, X‘ and Y°, of the chemical reactions at the 
two interfaces or the change in free energy of the system per unit rate of 
reaction. Thus, we may write 


T d;S/dt =J,A 


where A, the affinity, is given by the sum of the chemical potentials (grad), 
ents) of the two carrier forms. 

In short, when K3 #K3, a static head will be achieved where c,° #c! and 
free energy must be dissipated by the system to maintain this asymmetry 
The energy is invested to maintain gradients of free carrier and the rate of 
energy expenditure is determined by the magnitudes of these gradients ay 
the rate at which they tend to be dissipated by recycling of X and Y. 

Finally, it should be noted that the maximum osmotic work that can bw 
performed by this system is given by 


i 
+ 
RT In a =RT in 22) 
Cj ] +K, 


and when K; >K, this is less than the free energy of the driving react\iw) 
(A), which is RT \In(K3/K,). The reason for this is that P, #0, so that 
gradients of X established by the chemical reaction can be dissipated by the 


movement of X through the barrier. If P, = 0, the system is tightly (ideully) 
coupled and the static head that will be achieved is given by (c;'/¢,")* 
(K3/K2). At this point all flows will cease inasmuch as Ay, is now equul anil 


opposite to A. If we now perturb the system so that (¢//c,°) > (K3/K,), the 
chemical reaction may be driven in the reverse direction (Mitchell, 1970) 
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6.6. Rheogenic carrier-mediated transport 


Up to this point we have limited our considerations to carrier- 
Mediated transport of uncharged solutes. We will now briefly consider the 
bases illustrated in Figure 6.5, where a monovalent cation i is transported 
ross a membrane after forming either (a) a positively charged complex 
(iX,.) with a neutral carrier (X) or (b) a neutral complex (iX) with a mono- 
falent anionic carrier (X_). We will use these examples to illustrate an ap- 
roach toward the formal description of such systems as well as several 
joints bearing on the distinction between carrier-mediated versus diffusional 
sport of charged species. Since these mechanisms result in a carrier- 
hediated transfer of charge, we refer to them as “rheogenic” or “current 
herating” (from the Greek, rheo = current or flow). 
We assume that the presence of an electrical potential difference across the 
brane only affects the translocation rate of the charged carrier form(s) 
d has no effect on the association-dissociation reactions at the interfaces. 
thus, the equations that describe the association-dissociation reactions in 


e systems are 
K[iX°] - c;° [X°] =0 

K[iX"] - ¢/[x‘] =0 
ii the conservation of carrier is given by 

[Xe] + [X?] + 2X4] + P] = 14] 

re [X,] is constant; these equations are identical to those that describe 
rged carrier systems. 
if the first-order translocation rate constant of the neutral carrier forms is 
n by P, then the rate constant for translocation of the charged carrier 
ns must be modified by the factor &% or {@-) where & = exp( FAW/RT). 
value of a must lie between 0 and 1 and is a function of the electrical 
with the membrane. If the electrical potential barrier is symmetrical 
out the midpoint of the membrane, the rate constants for translocation 
wn Oo to i andi to o will be equally and oppositely affected so that a = 0,5 

Geck and Heinz, 1976; Stein, 1977; Hall et al., 1973). Under these 


ditions: 
{a) When the carrier is neutral, the steady state is described by 


P{[X?] - [x4] + €-Y? [ix,7] - EY [iX,'1} = 0 
the unidirectional flux of i from compartment o to compartment 7 is 
P¢~'/2(ix,°]; the unidirectional flux of i in the opposite direction 


n by PE*? [ix,/]. 
{h) if the carrier is a monovalent anion so that the complex iX is neutral, 
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the steady state is described by 


“ye x 
afr e Pie! [x2] - 9x4] + eX?) - EX} =0 
o 
2) ff = 3 Solving these equations, we obtain: 
lin) x< =f (a) When the free carrier is neutral, 
= 
: Jit = {PUXi) o(°(KEN? + HD (6.36) 
i") 
= z Sf? = {P[X,] ¢((KE"? + ¢4°)}1D (6.37) 
s {l/s = eee 
bead fan —— x Or-1/2 _ , igt/2 
aj; a =. Jp= {P(X 1] K(eiPE' - cf §" HD (6.38) 
5 3g here 
i ¥ y Bs =2K? +Ke(1 a g712) +Ke/( iy £12) + Pc f(E'? + gti?) (6.39) 
: a as ae. Lhe 
e = ES (b) When the free carrier is a monovalent anion J;°', J;'°, and J; are also 
5s n by equations (6.36) to (6.38), but D is given by 
— & 
o-- <5 BS 2ejPc! + Kote?) + Ket e7)+ KEP +E?) (6.40) 
i Se Several points should be noted: 
a3 First, in both cases J;=0 when (c;%c;/) =exp(FAW/RT), that is, the 
s . mst equilibrium distribution, as must be the case for any equilibrating 
nara 2 2 tem. 
oO 36 Becond, when Ay =0, equations (6.36) to (6.38) reduce to equations 
“+ 3 s }.11) for facilitated transfer of an uncharged solute. 
Pie Third, when c;? =c;!, the relation between J; and Ay is described by an 
- Mids By symmetric curve about the origin. 
. s& Finally, in Figure 6.6 the relations between J; and Ay calculated from 
x< eS ations (6.38) to (6.40) are compared with that predicted by the “con- 
2 Sb t field” flux equation (2.13). In the example shown in Figure 6.6, c;° and 
=i Ea were fixed at 100 mM and 10 mM, respectively; P was set equal to 0.01 
Oo 83 & co har i/h (as in Figure 2.3); K was set equal to 1 mM and X, was adjusted so that 
tap a? eo oe values calculated from equations (6.38) to (6.40) when Ay =0 were 
oO 2 F al to the value given by equation (2.13) when Ay =0. Clearly, over the 
24 Aw =+100 mV the three curves are quite distinctive, but over the 
: < oa oli Aw =0 to Ay = 58 mV (the equilibrium or reversal potential) they are 
< <———— os tinguishable. Further, it can be readily shown that as c;' approaches c;° 
ox 3 Ey curves become increasingly linear, the intersections approach the origin, 
iz > one can choose parameters that expand the range over which they are 
%: tially indistinguishable. 
re are important practical messages to be gained from this exercise. As 
a be discussed further (Section 7.1), experimentally determined relations 
ween J; (or J;) and Ay (so-called current-voltage or /-V relations) can 


Ide a considerable amount of information regarding the properties of a 


ay 
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80=— —|00 


“—< 


mae 


GHK 
——-— Uncharged Carrier 


—-—-— Charged Carrier 


=100 mM 
60 


c!=1om™M 


° 


C* 
-mediated transport of a monovalent cation by an uncharged or 


The arrow indicates the “reversal potential” (58 mV) when the concentration ratio 


J; and Aw predicted by the GHK equation for simple diffusion of a monovalent 


— 


sa=0 


estion and by equations (6.38) to (6.40) for carrier 


Figure 6.6. Relations between 


checgeed carrer © 


Carrier-mediated transport 113 


system and the underlying transport mechanism. Thus, if only c;° is known 
and there is reason to believe that J; is not influenced by nonconjugate driv- 
ing forces, the value of Ay when J;= 0 (the reversal potential) permits us to 
calculate c;! from the Nernst relation; on the other hand, if both c,? and c/ 
are known and (Ay)j,-o is not equal to (RT/F) In(c,7e?), the transport 
process must be influenced by a nonconjugate driving force. This is true 
regardless of the underlying transport mechanism! 

Further, as shown in Figure 6.6, the relations between J; and Ay predicted 
by the constant-field equation and by the two rheogenic carrier mechanisms 
are quite distinctive providing studies are carried out over a sufficiently wide 
range of Aw’s. As discussed in Section 2.3, if J; is given by the constant-field 
equation, the value of J; when Ay =0 is a measure of P; so that knowing 
(AY)z,<0 and (J;)ay =o (ie., ¢°, cj, and P;) the entire curve can be con- 
structed. Significant departures from this predicted curve may suggest the 
involvement of a neutral ion carrier or a charged carrier. However, it must be 
reemphasized that the relations between J; and Ay predicted by these three 
quite different transport mechanisms may be indistinguishable over a wide 
Tange, so that they must be interpreted cautiously .* 


6.7. Distinction between carriers and pores 


By this point, it should be clear that pore-mediated transport and 
€arrier-mediated transport may display a number of common features, such 
as saturation kinetics, competitive inhibition, and transeffects. Further, as 
discussed in Section 6.6, the relations between J; and Ay for charged solutes 
may not provide clear grounds for distinguishing between simple diffusion 
and carrier-mediated transport. How, then, can one distinguish between these 
two different transport mechanisms? 

Considerable insight into this problem has emerged from studies of artificial 
lipid bilayers “doped” with antibiotics that are capable of rendering these 
barriers permeable to ions; in the absence of these agents, lipid bilayers 
‘posses extremely high electrical resistances (>10® Q - cm?) indicating that 
they are essentially impermeable to hydrophilic charged species. Two groups 
‘of such agents have been identified: (a) those that serve as mobile carriers or 
“fonophores” (true “ferry boats” in the classical sense) such as nonactin, 
nomycin, and monensin; and (b) those that form water-filled pores that 
n the entire thickness of the barrier such as gramacidin A and amphoteri- 
B. Much of the experimental evidence that ion transport through pores 

exhibit saturation kinetics, competitive inhibition and transeffects is 
tived from studies on artificial lipid bilayers treated with gramacidin A 
f, Eisenman, 1973; Urban et al., 1978; Hladky, 1974). 

Studies on nonactin-doped lipid membranes indicate that each carrier 
ecule can transport, at most, about 10* ions per second across a bilayer 
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approximately 100 A thick (Lauger, 1972). On the other hand, studies on 
lipid membranes doped with pore-forming agents (as well as theoretical con- 
siderations) indicate that each pore can conduct as many as 10° to 10” ions 
per second, a value two to three orders of magnitude greater than the maxi- 
mum “turnover” number of Anown carriers (Lauger, 1972; Hladky, 1974). 

If these arguments can be extrapolated to biological membranes, measure- 
ments of single-site conductance could distinguish between pores and mobile 
carriers. Estimates of the conductances of single Na entry sites in a variety 
of excitable membranes indicate that each site is capable of transporting 
approximately 10° ions/sec, suggesting that passive Na entry into these cells 
is mediated by pores or channels rather than carriers (Armstrong, 1975). 

However, it must be emphasized that the notion that the turnover number 

of carriers is low compared to that of pores is based on experiments with 
artificial lipid membranes, approximately 100 A thick, doped with iono- 
phores whose diameters are approximately 10 to 15 A (cf. Simon and Morf, 
1973). Thus, these carriers, in fact, are small ferries that must diffuse a dis 
tance at least seven times their size through a viscous milieu in the course 
of shuttling one ion across the barrier and then must retum for the next 
passenger; this is a time-consuming voyage and any lateral movements in the 
plane of the membrane (side trips) will increase its duration. It is reasonable 
to question whether these artificial carrier systems are realistic models of 
purported carriers in biological membranes. This question cannot be answered 
ai this time and, in the absence of evidence to the contrary, we can certainly 
conceive of carrier models that do not involve the shuttling of a binding site 
back and forth across the entire thickness of the membrane. For example, 
pore with a recognition site that alternately faces one side of the membrane 
and then the other would, according to our definition, be considered » 
“carrier” and might exhibit turnover times that equal or rival those observe! 
in “true” pores. Thus, the use of “turnover number” to distinguish betwee 
pores and carriers in biological membranes may not be generally valid. 

As already discussed, the essential property of carriers, as they are gener 
ally conceptualized, is that their recognition or binding sites are capable uf 
undergoing transitions that shift their accessibility from one side of the mew) 
brane to the other. In contrast, we view pores as being simultaneously acces 
sible to solutes from both sides. Criteria for distinguishing between pores ail 
carriczs should be direct consequences of this fundamental difference. \n th 
sense, the phenomenon of “uphill” counter-transport is the most widely 
accepted positive criterion for the implication of a carrier rather than 4 pore 

(Widdas, 1952; Rosenberg and Wilbrandt, 1957; Wilbrandt and Rosenberg. 

1961; LeFevre, 1974). It is difficult to see how, in general, the downhill flaw 

of a solute through a pore in one direction can energize the flow of anolle 

solute in the opposite direction; yet, this phenomenon has been widely ol) 

served experimentally (LeFevre, 1974) and can be readily accounted for us)ig 


simple carrier models. 
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In addition, it is difficult to see how some well-documented examples of 
cotransport, where the downhill flow of one solute energizes the uphill flow 
of another in the same direction, can be accommodated by a pore model. 
Certainly, transport through pores can exhibit positive coupling. An excellent 
example of this phenomenon is solvent drag (Andersen and Ussing, 1957). In 
addition, there are several examples of uphill cotransport that appear to be 
attributable to solute-solute interactions (solute drag) (Ussing, 1966; Franz 
and Van Bruggen, 1967; Biber and Curran, 1968), which are presumably the 
result of frictional interactions (momentum transfer) that take place in 
water-filled channels and do not involve carriers. However, in all of these 
instances of cotransport caused by purely mechanical interactions, the 
toichiometry or coupling coefficient is highly variable. On the other hand, 
there are a number of well-documented cotransport processes where the 
coupling coefficients appear to be fixed and independent of the concentra- 
tions of either of the cotransported species or their respective conjugate 
driving forces (Schultz and Curran, 1970). These phenomena can be readily 
described using carrier models that are simple extensions of well-established 
mzyme reactions; it is far more difficult to see how they can be the results 
of transport through pores.® 


8. Summary 


The carrier hypothesis has led a precarious existence since its incep- 
in 1933-35 (Osterhout, 1933, 1935; Lundegardh, 1935); nonetheless, 
has survived and now appears to be flourishing. There is a conglomerate of 
etic evidence consistent with this notion (heroically summarized by 
Fevre, 1974), and there are no equally acceptable alternatives for many 
these observations. During the past fifteen years this area has embarked on 
“new direction aimed at isolating purported carriers and characterizing their 
avior after reconstitution in artificial lipid membranes. It is hoped that as 
fesult of these efforts, abstractions that were created to satisfy kinetic 
rvations will acquire indisputable biochemical identities. 

In this chapter we introduced the types of transport phenomena generally 
ributed to carriers using a kinetic approach. In some instances these for- 
$ are qualitatively correct but do not hold up quantitatively under 
us testing; such inconsistencies, or “anomalies,” and possible explana- 
are summarized by LeFevre (1974). Suffice it to say, that the present 
of development of this area resembles that of enzymology before 
mes were finally isolated, purified, and characterized with respect to 
ture and mechanism of action; many of these problems will not be 
ved until this area of membranology achieves the current state of the art 


enzymology. 


7 Some principles of 
electrophysiology 


The underlying theme of this monograph is that the formal description of all 
transport phenomena centers about the determinations of flows and forces 
that are then related by means of generalized conductances or resistances. Up 
to this point we have focused on the relations between fluxes (determined 
using chemical or isotopic techniques), electrochemical potential differences, 
and permeabilities, that is, the kinetic approach. In this final chapter we will 
briefly consider the electrophysiologic approach toward the same end, for two 
reasons. First, whereas many students of the field are comfortable with the 
concepts of fluxes, concentration differences, electrical potential differences, 
and permeabilities, they often react with horror when faced with equivalen|! 
electrical circuits that involve currents, resistors, and electromotive forces 
(batteries); our aim will be to show that these are in fact two ways of express 
ing the very same concepts. Second, whereas the kinetic and electrophysi«. 
logic approaches share the same theoretical base and differ only in terminol. 
ogy, the latter enjoys a very advanced technology and, as a result, a number 
of practical advantages. For example, currents as low as 10°’? A can be 
measured accurately with almost instantaneous time resolution; this is equiva: 
lent to the net flux of approximately 1 X 107!” equivalents of an ion and 
there is no chemical or isotopic technique that can come close to matching 
this level of resolution. There are a number of other technical advantages 
offered by the electrophysiologic approach that need not be belabored, 
suffice it to say that the brilliant studies of Hodgkin, Huxley, Katz, Cole, ani 
others that led to the unraveling of many of the mysteries of the axon (Cole, 
1968) could not have been accomplished using other approaches.’ 


Cake Electrical circuit analogs of flows and forces 


Let us begin by assuming that the flow of an ion i is strictly diff, 
sional so that, in general, it can be described by the Nernst-Planck equation 


J, = -uj;t;2; F((RT/z;F) (@ |n G/dx) + (d/dx)] (7.1) 
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As discussed earlier, the electrical current resulting from this flow is given 
by J; = z;FJ;, so that equation (7.1) can be rearranged to give 


ar ee (=) (4%) . (<2) 
uj;2z;? FE; 2,F dx ax G2) 


Assuming that J; is constant (steady state), equation (7.2) can be integrated 
over the thickness of the membrane yielding 


oe ee ek 
uf eames n( 2) ae (7.3) 


where, as before, Ay = y! - y?. It should be noted that the right-hand side 
of equation (7.3) refers to properties of the external solutions; as discussed 
previously, this is valid providing Ay includes boundary (interfacial) electrical 
potential differences. 

Now, following Finkelstein and Mauro (1963), we define the integral resis- 
tance of the membrane to i as 


i 
Ri = f (dx /u;z;? F 7Z;) (7.4) 
Oo 
and the electromotive force of i as the Nernst potential, 
E; = (RT/z;F) in(c;%e;) = ~ (Au; /ziF) (7.5) 
Substituting equations (7.4) and (7.5) into equation (7.3), we obtain 
I,R;=E;- Aw (7.6) 


Equation (7.6) can be readily recognized as a statement of Ohm’s law and, 
more generally, is an example of the linear phenomenological equation that 
felates flows and forces (1.33). In essence, this equation states that the 
Nernst-Planck equation for the flow of an ion driven by its electrochemical 
potential difference can be represented (symbolized) by a current J; flowing 
fcross a resistor R,; that is connected in series with a battery E£; as illustrated 
in Figure 7.1(a); the measured electrical potential difference across the 
harrier is obviously given by 


yi - y° = Ay =E;- 1,R; 


Clearly, equation (7.4) cannot be integrated without some assumptions 
fegarding u,; and the concentration profile of i across the membrane (i.e., the 
felation between @; and x). However, insight into the meaning of R; can be 
pained by considering the case where u; is constant and ¢,° = c;', so that ¢c; in 


#quation (7.5) can be replaced with an average concentration ¢, that is equal 


(b) 


Figure 7.1. (a) Equivalent electrical circuit representation of the 
Nernst-Planck equation for ionic diffusion. (b) Equivalent electrical 
circuit for the diffusion of n ions across a membrane. (Note: As dis: 
cussed in Section 7.2, the electromotive force E can be generalized 
to include flows that are not accounted for by simple diffusion.) 
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to [8;(c;° + c;)/2]. Then, 


R; = [Ax/u;2z;?F7é,] my) 
Since, P; = 8; RTu;/Ax, the relation between R; and P, is 
Rj = [RT|P,2;? F263] ) ") 


The conductance of the membrane to / is, simply, G; = (1/R,) so that 
G; = [Pj2z;’ F°¢,/RT] (7.9) 


Two important points should be noted: 

First, the electrical resistance of the membrane to / is directly related to its 
thickness and inversely related to the local concentrations of i integrated over 
the membrane thickness, or the content of i per unit area.” 

Second, although G, is directly related to P;, these terms are not synony- 
mous inasmuch as G; is also dependent upon ¢;. This point, which will be 
illustrated below, must be emphasized since the terms conductance and 
permeability are often used interchangeably in the literature. 

Now, if we assume that the electric field within the membrane is constant, 
equation (7.1) can be integrated directly, yielding 


Ee. (Zats’a¥) [e - oj exp(z; Zable) 


RT 1 - exp(z; ¥ Ay/RT) rte) 


Further, from equation (7.6) we see that when £; is constant (i.e., when c;° 
and c/ are constant), R; = -[8(Aw)/d/;], so that 


. --[ (RT) ee] 


P,z;°F? Aw 7? - cf é i 


where & = exp(z; FAy/RT). 

Clearly, multiplying equation (7.10) by equation (7.11) yields equation 
(7.6). ' 

When c¢;? = c; then J, = -[(P;z;°?¥? Ayc;)/RT] and R; approaches the ex- 
pression given by equation (7.8). f 

Further, when Ay = 0, then J; = P;z;F(c;° ~ c;') and 


R= ( at ) fasscter (7.12) 


Piz7°F*?} (ci? - ¢?) 


As discussed previously (Section 4.2), the term [Ae;/In(c;?/e;')] is the loga- 
rithmic mean concentration of i, which approaches the arithmetic mean ¢; 
when Ac; > 0, so that equation (7.12) approaches equation (7.8). 

We now return to an issue that was introduced earlier (Section 2.3) and 
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_ deferred for further discussion, namely, “Goldman-type” rectification. Refer- 
’ ring to Figure 2.3, we see that according to the constant field equation a plot 
of J; versus Ay is linear only when c,? = c;; when ¢;? # c;, the plot is curvi- 
linear in spite of the fact that P;, c;° and c;! are constant. Now, the ordinate 
labeled J; can be relabeled J; (=z;¥J;) so that the curves actually are current- 
voltage or “J-V” plots. The slopes of these curves, at any value of Ay are 
given by [0/;/0(Ay)] (the partial differential is employed to indicate that P;, 
¢;°, and c;! are maintained constant) and are equal to G; (or 1/R;). Clearly, 
when ¢,? #c;, the slope conductance G; is a function of Ay and is lower 
when J; is directed from the compartment with the lower concentration to 
that with the higher concentration (in the example shown, from compart- 
ment i to compartment o so that J; is negative) than when J; is in the opposite 
direction. When c,° = c;', G; is independent of Ay and the relation between J; 
and Ay is linear. 

Figure 7.2 is a plot of R; versus Ay for the conditions illustrated in Figure 
2.3, calculated using equation (7.11); these curves clearly confirm the im- 
pressions gained from inspection of Figure 2.3. Thus, when c;° = 100 mM 
and c;! = 10 mM, R; decreases as Ay is varied from +100 to -100 mV; when 
cf = 50 mM, R; is always less than when c;' = 10 mM for any finite value of 
Ay but still decreases when Ay is varied from +100 to - 100 mV; when ¢c,° * 
c} = 100 mM, R; is still lower and is now independent of Ay [see equation 
(7.8)]. 

Figure 7.3 provides the explanation for these phenomena. This figure illus 
trates the concentration profiles of i within the membrane under several con 
ditions calculated according to Arndt and Roper (1972). 

When c;° = 100 mM, c;/ = 10 mM and Ay = +100 mV, the concentration 
profile is concave upward; since J; is directed from compartment i to com 
partment o, this profile can be viewed as the result of a convective (fluid) 
flow of current, interacting with i. Conversely, when Ay =-100 mV, J; \s 
directed from compartment o to compartment i and the concentration profile 
is convex upward. When Ay =0, the concentration profile is linear. It |» 
obvious that the amount of i within the membrane (per unit area) is lens! 
when Ay =100 mV, most when Ay =-100 mV and intermediate whey 
Aw = 0. From equation (7.4) we see that R; is inversely related to the amoup| 
of (mobile) i within the membrane. It follows that R, will be greater when 
Aw = 100 mV than when Ay = -100 mV and that an intermediate value of 
R; will be observed when Ay = 0. 

When ¢;? =c;! = 100 mM, the concentration of i within the membrane \§ 
independent of position and Aw; consequently, R, is independent of Ay 
Further, it is obvious that the amount of i within the membrane under (hii 
condition (for any finite value of Ay) is always greater than that whew 
c; = 10 mM so that R; will always be lower (see Figure 7.2). 
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Figure 7.2. R; as a function of Ay for a monovalent cation when c;° and P; are maintained constant and c; is varied. 
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Figure 7.3. Concentration profiles within the membrane when 
cj? = 100 mM and c/ = 10 mA | for Ay = +100 mV, 0 mV, and 
-100 mV (solid lines). When c,? = c;! = 100 mM, ¢; = 100 mM 
throughout the membrane a is independent of Aw (dashed line), 
(Note: For the purposes of illustrating this point, we have assumed 
that BY? = B; = 
dade The “overall” resistance and electromotive force 


If there are parallel flows of n ions across a membrane, equation 
(7.6) can be written as follows: 


n n 

ISA 

i i 

Now, in this parallel system, the total conductance of the membrane G,, !¥ 
equal to Z7 G;, so that equation (7.13) reduces to 


C, (7.13) 


es es 


m=, (GiEiIGm)~ OY (7.14) 
I 
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where R,,, is the total electrical resistance of the membrane and is given by 


“ -1 
Rm = b ciry| = (1/Gin) 


G; is referred to as the partial ionic conductance of i, inasmuch as it repre- 
sents the contribution of J; to the total membrane conductance. 

Clearly, from equations (7.6) and (7.13) it follows that G; = [a/;/0(-Ay)] 
and G,,, = [d//8(-Aw)] so that (G;/G,,,) = (01;/07). In other words, (G;/Gyn) 
is the ratio between the change in J; to the change in total current / in re- 
sponse to an imposed Ay when all E, (i.¢., c;° and c;') are maintained con- 
stant; as discussed earlier, this is the transference number t; of the total 


membrane to the species i under the given set of conditions.* Thus, equation 
(7.14) can be rewritten as follows: 
IRm =E - Ay (7.15) 


where £ is the total electromotive force across the membrane given by 


(7.16) 


If the ionic flows are strictly diffusional, £; is given by equation (75) 
which, when substituted into equations (7.15) and (7.16), yields 


-RT S (ti\ (ee\ - 
Rm = 5 u (4) m(2) Ay (7.17) 
us, when J = 0 
ry. 1 &/4),,. 
a FE (Hale) 5 SG )am os 


Equation (7.18) is identical to equation (3.18) which describes the multi- 
ic diffusion potential under “zero current conditions.” In other words, E 
the electrical potential difference across the membrane under zero current 
ditions or the electrical potential difference at which the membrane must 
clamped in order to abolish current flow. 

Clearly, equation (7.17) can be represented by the equivalent electrical cir- 
jt illustrated in Figure 7.1(b), where, in general, R;, Rj, Rx,...,R, and 
&), Ex, ..., En ate given by equations (7.4) and (7.5); and Jj, Jj, 14,--+1n 
given by equation (7.6). In the particular case where the electric field 
thin the membrane is constant, R;, Rj, Ry, . . .» Ry are given by equation 
WN) It Ips Tes «++ Tn ate given by equation (7.10), and Ej, Ej, Ey,...,En 
given by equation (7.5). 
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Finally, the expression for the overall electromotive force can be general. 
ized to include the possible effects of nonconjugate driving forces on J;. In 
short, any process that can generate a flow (current) of i will contribute to 
E; and E. These include the coupling of J; to solvent flow (solvent drag); the 
coupling of J, to the flow of another solute, j (solute-solute cotransport, 
solute-drag); the coupling of J; to a metabolic reaction (active rheogenic 


transport); and $0 ON. Thus, in general, 


(OE 


1. 42. 5 
+ ea z;F(1 - 0;) éJy (Streaming potential) 
mi 


(Diffusion potential) 


(7.19) 


n 
Z s Zs zjsvd; (Cotransport potential) 


mi 

loan 
+G > 2i FyJ,+°-- (Active transport potential) 

im i 
where v is the coupling coefficient between J; and J; and ,J; is the flow of / 
resulting from active transport. ; 

Thus, the equations of flow that have been discussed in previous sections 
can be translated into measurable electrical parameters such as currents, 
resistances, batteries (electromotive forces), and electrical potential diffey. 
ences and equations (7.15) and (7.19) are simply the Thévenin equivalents 
of a system characterized by the flow of one or more ions; Figure 7.1(b) \s 
the electrical circuit representation of such a system. 


73. The “short-circuit technique” 


Although the emphasis in this monograph has been on principles 
rather than techniques, the short-circuit technique has such instructive 4» 
well as practical value that a would be remiss if it were not discussed in somw 
detail. This technique was introduced by Ussing and Zerahn (1951) for the 
study of Na transport by isolated frog skin but, whereas it has been employe! 
almost exclusively for the investigation of ion transport by epithelial tissues, 

in principle, its applicable to any membrane. 
Let us first consider the theoretical basis of this technique and then illuy. 

i jcation. 

aon is bathed on both surfaces by a solution of a salt that diy. 
sociates into 4 cation j with valence z, and an anion j with a valence 2), the 
system can b¢ represented by the equivalent circuit illustrated in Figure 7.4 
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Figure 7.4. Equivalent electrical circuit for the transport of a cation 
iand anion j across a membrane under open-circuit and short-circuit 
conditions. 


As shown, provisions are made to measure the electrical potential difference 
across the membrane using the voltmeter V (which is assumed to draw no 
current) and to pass an external current J across the barrier driven by a 
variable battery B. 

The general solution of this circuit is 


IRm =E - AW 


where Ry, = R/R;/(R;+R,) and E = Rp, ((E;/R;) + (Ej/R;)] or alternatively 
[equation (7.16)] 


E = (G;/Gy,) E; + (G;/Gm) Ej 
When J = 0 (i.e., the switch S is open) 
Ay = £= Ry, {(E;/R;) + (E/R,)I (7.20) 


Thus, assuming that R; and R; are finite, there will be an electrical potential 
difference across the membrane if E, and/or E; are nonzero; this is referred to 
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as the open circuit potential. This electrical potential difference can be 
abolished by passing an Opposing current across the membrane such that 
IRm = - Aw. Clearly, the magnitude of this current is given by 


()ay=o0 = E;/R;) + (E;/R;) (7.21) 


an when both bathing solutions have identical compositions so that 
Cy = ey =c; and ¢° = c; = cj, equations (7.19) and (7.20) yield (when J = 0) 


AV = Rm [Ei/Ri) + Ej/R,)] = RmziF(Si)acj=0 + Rmzj)F(Jj)acj=0 


or 


Ay =Rm Ci)ac;=0 1 G)acj=o] (7.22) 
so that 


A 
(ac, ayeo= Re = [Aaciao + Dacjao (7.23) 


Thus, when both bathing solutions have identical compositions, the current 
necessary to abolish the transmembrane electrical potential difference is equal 
to the algebraic sum of the currents of i and j that are driven by nonconjugate 
forces; this current is referred to as the short-circuit current, which we denote 
1 oa 

In other words, if i and j can only cross the membrane by diffusion, then, 
when both solutions have identical compositions, there will be no flows, 
£=0 (7.19), Ay =0 (7.18), and obviously Tse =O (7.21). However, if / 
and/or j are driven by nonconjugate forces, then, when Ac; = Ac; = 0, in 
general, F #0, Ay #0, and the J,. is given by the algebraic sum of the cur. 
rents of i and j. 

In order to illustrate these points further, let us consider the system shown 
in Figure 7.4, assuming that the anion j can only cross the membrane by 
simple diffusion and that the cation i is actively transported from compart: 
ment o to compartment i. When c;? = cf and cj” = cj , then £; = 0, and from 
equation (7.22) under open-circuit conditions (i.e., when / = 0), 


Ay =E=R,,E;/R; (7:24) 


Since E = Ry, Fz; J; [equation (7.19)], it follows that Ay > 0; that is, the 
inner solution will be electrically positive with respect to the outer solution, 

Now, what are the individual currents through the two parallel pathways 
under open-circuit conditions? These may be described as follows 


1,R;=E;,- Ay (7.25) 
and 


IR; = -Ay (7.26) 
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Thus, the flow of j from the outer to the inner solution is driven solely by the 
electrical potential difference established by the active transport of i. Clearly, 
in order that electroneutrality be maintained under open-circuit conditions, 
the restriction on the currents is that 7;=-J;. Using this restriction and 
equations (7.25) and (7.26), we obtain equation (7.24). Thus, the physical 
system is consistent with the equivalent electrical circuit model; the reason 
for the emphasis on “consistent with” will be discussed in Section 7.4. 

Now, the amount of an (opposing) current that must be passed across the 
membrane from the external source B in order to clamp Ay to zero is simply 


Tse = (E/Rm) = (E;/Ri) = 21F 05; (7.27) 


Thus, the short-circuit current is equal to the rate at which i is actively 
transported from compartment o to compartment i, when 4J; (in moles per 
centimeter squared, second ) is converted into amperes per centimeter squared 
by the factor z;F. 

Finally, it should be pointed out that when Ay = 0, there is no current flow 
through R; so that the entire J,, passes through the pathway responsible for 
active transport of 7 but in the opposite direction. Thus, the J,, may be 
viewed as the external current that must be passed through the membrane in 
order to “buck out’’ the current generated within the membrane by the active 
transport mechanism and thereby render Ay = 0.* 

Now let us descend from theory to practice and examine the application of 
the short-circuit technique to a “real” membrane where Na is actively trans- 
ported from the outer to the inner solution and the movements of Cl are 
entirely diffusional; this system is illustrated in Figure 7.5. As already dis- 
cussed, when both surfaces of the membrane are bathed with identical NaCl 
solutions, under open-circuit conditions the inner solution will be electrically 
positive with respect to the outer solution and 


AW =Rm FadNa 


Now, since electroneutrality must be maintained, there must be a flow 
of Cl from the outer solution to the inner solution equal to gJyag. Since 
[C1]? = [C1j', this flow must be driven entirely by Ay; thus, the electrical 
potential difference secondary to the active transport of Na is the driving 
force for the accompanying flow of Cl]. The magnitude of Ay will be directly 
related to the rate of active Na transport and the resistance of the tissue to Cl 
diffusion. [The resemblance of this line of reasoning to that underlying the 
diffusion potential between two solutions of the same salt (Section 3.1) 
should be noted. In that case, Ay was determined by the difference between 
anion and cation mobilities and the magnitude of the concentration ratio. 
In this case, there is no concentration difference but Ay is determined by 
properties of the system related to the maintenance of electroneutrality.] 
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INNER 
SOLUTION 


Figure 7.5. Illustration of the application of the short-circuit tech- 
nique to the study of a membrane that is capable of actively trans- 
porting Na from the outer to the inner solution; Cl transport is 
assumed to be entirely diffusional. 


When sufficient current is passed across the membrane to abolish the open 
circuit Ay, clearly, Ati, = Afic) = 0. This is referred to as the short-circuit 
condition, inasmuch as the electrical potential difference across the mem. 
brane has been “shorted”; it is as if the membrane were bypassed with an 
infinitely high conductor (infinitely low resistor). Because Na is actively 
transported, a net flow of Na from the outer to the inner solution will persist 
in spite of the fact that Afiy, = 0. On the other hand, because the flow of Cl 
is entirely diffusional, Jc, will cease when Ajigy = 0. 

How, then, is electroneutrality maintained under short-circuit conditions? 

Let us consider how an external current is passed across a membrane, sur- 
rounded by electrolyte solutions, from an electromotive force (the battery) 
via electrodes immersed in the two solutions. In solid-state or metallic con- 
ductors, current is carried by electrons, but in solutions it is carried by ions 
(electrolytic conduction). It follows that at the interfaces between the 
current-passing electrodes and the solutions in Figure 7.5, processes take 
place that convert ions into electrons and vice versa; these are referred to as 
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electrode reactions. Now, let us assume that the current-passing electrodes 
consist of silver wires or plates that have been electroplated with Cl, that is, 
silver-silver chloride (Ag/AgCl) electrodes. At the cathode the electrode reac- 
tion is 


AgCl+e —>Ag+Cl” 


so that an electron driven down the wire from the negative pole of the bat- 
tery generates a Cl ion that enters the solution. At the anode, the opposite 
teaction takes place, that is, 


Cl + Ag —> AgCl + e~ 


so that a Cl ion is “plated out” of solution; the electron released is drawn to 
the positive pole of the battery and the circuit is completed! 

Thus, when sufficient current is passed across the tissue to abolish AW and 
tender Aficy = 0, Cl transport through the membrane ceases but there is a net 
displacement of Cl from the outer to the inner solution resulting from the 
electrode reactions. In essence, for every Na actively transported across the 
membrane from compartment o to compartment i, one Cl in the outer solu- 
tion (which contains the anode) is plated out and the electron that is released 
flows through the ammeter to the positive “sink” of the battery. At the same 
time an electron derived from the negative “source” of the battery interacts 
with AgCl at the cathode resulting in the release of Cl ion into the inner solu- 
tion. Thus, the rate of flow of electrons through the extemal circuit is pre- 
cisely equal to the flow of Na from the outer to the inner solution through 
the membrane, so that the /,, is equal to the active Na current and electro- 
neutrality is maintained. 

Obviously, the same reasoning can be applied to a membrane capable of 
actively transporting several ions having different valences in different direc- 
tions. It is a simple matter to show that the /,, (i.e., the rate of electron flow 
in the extemal circuit) is always the algebraic sum of the ionic currents 
through the membrane. 


74. A caveat regarding electrical circuit models 


In their important paper dealing with the use of equivalent electrical 
circuits to model ionic movements across membranes, Finkelstein and Mauro 
(1963) point out a significant shortcoming of this approach that must be kept 
in mind, Briefly, the problem arises from the fact that Kirchhoff’s laws, which 
are used to solve electrical circuits, deal exclusively with the flow of electrons 
but do not recognize different ionic species; that is, these laws do not dis- 
criminate between an electron carried by Cl and one carried by HCO3 nor do 
they discriminate between the current carried by the flow of a monovalent 
anion in one direction and that carried by the equal flow of a monovalent 
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cation in the opposite direction; these are electrically equivalent. But, artifi- 
cial and biological membranes do make these distinctions! 

For example, in Figure 7.1(b), there is no formal way of preventing the cur- 
rent J; from flowing across the resistor R j; according to Kirchhoff’s laws, 
currents are currents and resistors are resistors and the subscriptsi,j,k,...,n 
that designate different ionic species are unrecognized. With reference to 
Figure 7.4, we employed equations (7.25) and (7.26) to force J; through R; 
and J; through R; and then imposed the physical restriction that J; = -J;; 
these impositions are not natural consequences of a formal solution of the 
equivalent electrical circuit model. 

Some of the possible misleading consequences of applying equivalent elec- 
trical circuit models to multi-ionic systems have been discussed by Finkelstein 
and Mauro (1963) for single plasma membranes and by Schultz and others 
(1977) and Schultz (1979) for epithelial tissues. The purposes of these 
caveats are not to disparage or discourage the use of this valuable approach to 
modeling electrophysiologic properties of membranes. Instead, it must be 
emphasized that they must be employed with circumspection and that, often, 
restrictions must be placed on the model in order to transform a collection of 
wires, batteries, and resistors into a realistic portrayal of a biological system. 


Notes 


Chapter 1 


. There are many excellent texts on classical thermodynamics. My 


favorites are Lewis and Randall (1961), Glasstone (1956), Denbigh 
(1966), and Guggenheim (1950), the latter being a very rigorous 
treatment. For an unusual and delightful approach see Angrist and 
Hepler (1967). 


. We will emphasize closed systems (i.e., those that do not exchange 


matter with their surroundings) at constant temperature. 


. The partial differential (0x/0y),),y is the derivative of x with 


respect to y when other variables are maintained constant. Thus, the 
integral of [(0x/0y)y, v, w] dy from an initial state o to a final state i 
is the contribution of the change in y alone to the total change in x. 

In some instances for the sake of clarity we will omit the sub- 
scripts u,v, w,...and simply write (0x/dy). 


. We now refer to a “closed system” that can exchange heat but not 


matter with its surroundings. As will be discussed briefly, matter also 
“contains” entropy so that in an “open system” changes in S, in gen- 
eral, result from exchanges of heat and matter with the surroundings. 


. It is a common misconception that, in general, a difference in temper- 


ature is the driving force for heat flow. This is true only for a “‘closed 
system”? which cannot exchange matter with its surroundings. Since 
matter has a heat capacity, a difference in temperature is a driving 
force for the flow of matter [see equation (1.29)] and the exchange 
of entropy between an “open system” and its surroundings involves 
both the flow of heat and the flow of matter. 

Further, it should be noted that P in equation (1.5) now refers to 
the internal pressure of the system rather than the external pressure 
of equation (1.2), since if the process is reversible, the internal and 
external pressures can only differ infinitesimally. 


. Josiah Willard Gibbs (1839-1903), a chemist, as well as a theoretical 


physicist, was particularly interested in reactions that take place at 
constant temperature and pressure and was not concerned with 
“work” resulting from a change in volume of a solution undergoing a 
chemical reaction; hence the implication that this (PdV) work is 
“useless,” 


. Clausius coined the term entropy which stems from the Greek, 


entrope-the act of turning toward or turning into. He clearly recog- 
nized that this function reflects the “directionality” of spontaneous 
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processes and is a measure of the degradation of a system’s ability 
to perform work. 


. Boltzmann’s statistical mechanical treatment of the second law made 


the increase in entropy which accompanies an irreversible process a 
highly probable occurrence rather than the absolute certainty en- 
visioned by Clausius. Planck initially rejected this probabilistic ap- 
proach and undertook his studies on ‘‘black-body”’ radiation, in part 
to establish the absoluteness of the concept of irreversibility. In 
the end, he was forced to accept the statistical basis of the second 
law and his studies, which were designed for an entirely different 
purpose, led to the quantum theory (M. J. Klein, Physics Today, 
19: 23-32, 1966). 


. It must be emphasized that all potentials (intensive properties) must 


be referred to an arbitrary “reference point,” or “ground state.” 
Thus, when we refer to a pressure, it is generally with respect to that 
at sea level which is assigned the value 0 mm Hg. The electrical 
potential at a point is with respect to “ground,” which is assigned 
the value 0 volts. And, as is well known, temperature on the centi- 
grade scale is referred to the freezing point of water (at atmospheric 
pressure), or on the Kelvin scale to the hypothetical point where the 
volume of an ideal gas vanishes (-273.16 C). A difference between 
two pressures is meaningful only when each pressure has the same 
reference point; the same is true for a difference in electrical poten- 
tial and a difference in temperature. Similarly, there is no such thing 
as an absolute electrochemical potential. All values of i; are with 
reference to a “standard (ground) state,” where T= 0, P= 0, y =0, 
and c; = | in a given solvent; clearly, under these conditions ji; = u;°. 
A difference in electrochemical potential difference is meaningful 
only when reference is made to the same standard (ground) state. 
For general references see DeGroot (1952), DeGroot and Mazur 
(1962), Katchalsky and Curran (1965) and Prigogine (1961); each is 
a classic in its own right. 


Chapter 2 


. Thus, J; is defined as positive when directed from compartment o to 


compartment i; this convention will be strictly adhered to through- 
out. It should be noted that for charged species the current flow 
attributable to J; is z; ¥J;, so that for anions the current J; is nega- 
tive when J; is positive. 


. For a discussion of the effects of unstirred layers see Lakshminara- 


yanaiah (1969) or Helfferich (1962). 


. We will only be concerned with constant or “steady-state”’ velocities. 


Under these conditions the total force acting on a particle is zero, so 
that the driving force is ‘“‘balanced”’ by retarding, frictional forces; the 
mobility is a measure of these frictional forces. For example, for dif- 
fusion of a spherical particle with radius 7; in a fluid whose viscosity 
is 7, vu; = (1/6mn7;)f;. Thus, u; is inversely proportional to the vis: 
cosity of the surrounding fluid and the radius of the particle (the 
Stokes-Einstein relation). (See L.G. Longsworth in Electrochemistry 
in biology and medicine (T. Shedlovsky, ed.), New York: John 
Wiley, 1955, pp. 225-47). 
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. For a solute with a partial molar volume of 20 cm? /mole in a solu- 


tion at atmospheric pressure, 0,P = 20(cm? /mole) X 10° (dynes/ 
cm?) = 2 J/mole. In contrast, RT at 25 C = 2500 J/mole, and fora 
monovalent cation z;,¥ = 96,500 J/mole when w is expressed in 
volts. 


. It is generally assumed that i is distributed at equilibrium across each 


interface. In essence this implies that diffusion through the mem- 
brane is slow compared to partitioning across the interfaces. Under 
these conditions we can write for each interface pw, = B;. From equa- 
tion (1.28) it follows that when z, = 0, 8; = ¢,/c; = exp [(u;° - B;°)/ 
RT]. Thus, 8; is a function of the difference in standard state free 
energy, which (at constant T and P) is a function of the different 
“solvent properties” of the membrane and the adjacent solution(s). 
If the two have the same solvent properties for i, B; = 1. On the other 
hand, if i is a lipophilic molecule, the outer solutions are aqueous, and 
the membrane is lipid, i will be more ‘“‘stable’”’ in the membrane 
phase; accordingly, 1;° <u,° and f; > 1. 


. Equation (2.21) can be derived from more general considerations 


when z; ¥ AW <RT (Parlin and Eyring, 1954; Kimizuka and Koketsu, 
1964) and when the potential energy profile is symmetrical about 
the midpoint of the membrane and the external electrical potential 
difference is less than the height of the energy barrier (Hall et ai., 
1973). An expression having a similar form has been derived by 
Jacquez (1971) for an arbitrary, continuous potential profile. 


. Integrated forms of the Nernst-Planck equation assume that the 


energy barriers within the membrane that must be overcome by the 
diffusing particle form a “continuum.” These integrations are not 
generally valid when there are discrete energy barriers which result in 
discontinuities in the electrochemical potential of the diffusing 
species i along the diffusion pathway; hence the term discontinuous 
diffusion (Hille, 1975). 


Chapter 3 


. The assumption that “bulk” or “macroscopic” electroneutrality is 


strictly preserved in solutions is central to all of our subsequent con- 
siderations. This, so-called Nernst-Planck electroneutrality condition 
can be expressed formally as follows: 


> 240,=- D5 2-c. 
c a 


The (overwhelming) justification for this assumption is discussed by 
Guggenheim (1950) and Sten-Knudsen (1978). In essence it can be 
shown that minute departures from electroneutrality would result 
in enormous changes in electric fields many times greater than those 
encountered in most electrodiffusional systems. For example, as dis- 
cussed by Sten-Knudsen (1978, p. 64), if the anion and cation con- 
centrations of a solution differed by 1 mEq/liter, the gradient of the 
electric field in that solution (d? /dx?) would be approximately 
10% V/m?, The electrical potential differences encountered in most 
electrodiffusional systems are less than 1 V and changes in the 
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electric field are not detectable. Thus, the electroneutrality condi- 
tion is an eminently reasonable assumption. 


. Formally, the assumption is that 6; = 6, = B_ and that 6; is indepen- 


dent of concentration. This assumption infers that the membrane is 
“neutral” and does not discriminate between the anion (-) and the 
cation (+). This assumption is not valid for ‘‘charged membranes,” 
which discriminate between anions and cations; a formal treatment 
for the diffusion of a single salt across charged membranes has been 
derived by Teorell (1953) (see Section 3.5). 


. It is a relatively simple matter to extend this derivation to the case 


where the anions and cations are multivalent and the valences have 
the same absolute value; that is, z, = |z_|. 


. Jacquez and Schultz (1974) have derived a general expression for 


the membrane potential of symmetrical cells under steady-state 
conditions which makes no assumptions regarding electrical potential 
profiles, partition coefficients, and so forth, and is valid in the pres- 
ence of “rheogenic pumps” (see Section 6.6). This expression is a 
generalization of one derived earlier by Mullins and Noda (1963). 


. There is some confusion in the literature over the terms transference 


number and transport number. Macinnes (1961, p. 59) and Robin- 
son and Stokes (1959, p. 44) do not distinguish between the two, 
whereas Helfferich (1962, p. 326) does. Strictly speaking, ¢; is de- 
fined as the fraction of the total current carried by the ith ion under 
conditions of uniform concentration driven by Ay alone (i.e., in the 
absence of interactions between i and the flow of solvent or other 
solutes). However, as discussed by Kirkwood (1954), it is quite 
legitimate to refer to the overall t; of a membrane that separates two 
solutions having different concentrations of i. Under these condi- 
tions ¢; is a complex function of the distributions of all diffusible 
ionic species in the membrane and the compositions of the surround- 
ing solutions [see equation (3.19)]; nonetheless, it can be defined 
empirically for any given set of conditions. 


. For a rigorous derivation, see MacInnes (1961, pp. 461-65). 
. Obviously, all membranes are electrically neutral! The terms 


“neutral” and “charged” simply distinguish between the manner by 
which electroneutrality is assured. In a “‘neutral” membrane, the 
structural components are either uncharged or bear an equivalent 
number of oppositely charged groups; thus, electroneutrality is 
maintained by fixed or unvarying elements of the membrane. On 
the other hand, a ‘‘charged’” membrane possess structural com- 
ponents that are capable of dissociating into a preponderance of 
fixed positive or fixed negative charges, and electroneutrality is 
maintained by the presence of free, unfixed counterions. The nature 
of these counterions will vary depending upon the compositions of 
the surrounding solutions (see Lakshminarayanaiah, 1969). 


Chapter 4 


. The mole fraction of i is the number of moles of i in a system divi- 


ded by the total number of moles of all substances. In a mixture of 
perfect gases at total pressure P, the partial pressure of gas i is 


P,= cjRT = x,P 
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(Glasstone, 1956). Thus, using the perfect gas analogy, once again, 
x, is proportional to c; and either term can be used in writing expres- 
sions for [;. 


. See Section 7.3 for a discussion of electrode reactions and Macinnes 


(1961) or Lakshminarayanaiah (1969) for discussions of electro- 
chemical cells and reversible electrodes. 


. See discussion in Section 7.3. 
. Recall that cj! > ¢;° so that Ac; < 0 and J; is negative. 


Chapter 6 


. For general discussions of carrier-mediated transport and carrier- 


kinetics see Wilbrandt and Rosenberg (1965), LeFevre (1975), 
Jacquez (1964) and Stein (1967). An introduction to carrier- 
kinetics can be found in Neames and Richards (1972). 


. See Blumenthal and Katchalsky (1969) for a treatment of a system 


where the interaction between the carrier and the transported solute 
can take place within the membrane and a discussion of the effect of 
carrier association-dissociation rates on transport. 


. For more general treatments of cotransport, including cases where 


the carrier complex is charged, see Stein (1967), Wilbrandt (1972), 
Heinz et al. (1972), Jacquez (1972), and Geck and Heinz (1976). 


. The approach we will take in this section follows that of Rosenberg 


and Wilbrandt (1955). 


. More general treatments of the kinetics of rheogenic carrier-me- 


diated transport are given by Lauger and Stark (1970), Hladky 
(1972), Hall et al. (1973), and Stein (1977). An analysis of the 
effect of electrical potential differences on ion-linked cotransport is 
given by Geck and Heinz (1976) and Jacquez (1972). 


. It should be emphasized that whereas the carrier model is perhaps 


the simplest explanation for “uphill” countertransport and cotrans- 
port, positive (cotransport) as well as negative (countertransport) 
coupling between tracer flows can, in principle, result from mem- 
brane heterogeneities and need not involve “carriers” (Kedem and 
Essig, 1965; Patlak and Rapoport, 1971; Li and Essig, 1976, 1977). 


Chapter 7 


. For an elegant and scintillating discussion of the application of elec- 


trical circuit models to biological membranes, see Cole (1968). 


. Clearly, when Ay = 0 and c;? = c;' the unidirectional fluxes across the 


membrane are simply J;% = Je = P;c;. Thus, under these (“short- 
circuit”) conditions G; = J,%(z;?¥?/RT). When z;= +1, J;! expressed 
in microequivalents per centimeter squared hour is almost numer- 
ically equal to G; expressed in millimhos per centimeter squared. 


3. From equation (7.4) it follows that the integral conductances are 


1 


G; = om 
J (dx /ujz,° FC) 


136 


Notes to p, 127 
and 


on 23 


1 


i 
| (dx /ujz,?F *¢;) 
oO 


Thus, in general, the overall t; is a function of the concentrations 
and distributions of all i in the membrane which, in turn, are influ- 
enced by the concentrations in the surrounding solutions and the 
partition coefficients (see Section 3.3). Clearly, when both solutions 
have identical compositions 


ae _ 452 ;7 Bye; 
(Cee ; 
> Mi217B¢; 
i 


which is identical with equation (3.19). 


. The direction of an electrical current is defined as positive when it 


flows from the anode to the cathode. Thus, strictly speaking, the 
Ts¢ is in the same direction as the active flow of i and, in fact, is the 
flow of i across the membrane. Accordingly, in essence, when 
Aw = 0, the active transport of i draws an equivalent amount of cur- 
rent from the battery through the external circuit to complete the 
internal circuit. However, 1 have found that it is difficult for “‘begin- 
ners’”’ to appreciate this concept whereas they readily accept the 
notion that the J,, ‘‘bucks out” the current resulting from a4;. This 
notion is contrary to the conventional definition of the direction of 
current flow but is conceptually equally valid. 
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Frequently used symbols 


affinity of a chemical reaction (calories/mole) 
concentration of species i in bulk solution (see Figure 2.1) 
concentration of species i in the membrane (see Figure 2.1) 
mean concentration of species / 

diffusion coefficien? of species i (cm? /sec) 

total electromotive force (volts) 

electromotive force of species # (volts) 

the Faraday (96,500 coulombs/equivalent) 

Gibbs free energy (calories) 

partial ionic conductance of species i (mhos/cm?) 

total membrane conductance (mhos/cm?) 

total current (A/cm?) 

current attributable to the flow of i (A/cm?) 

flux of species i (moles/cm? sec) (see Figure 2.1) 

rate of volume flow (cm/sec) 

hydraulic conductivity (filtration coefficient) (cm? /dyne sec) 
pressure (dynes/cm?; atm) 

permeability coefficient of membrane to species i (cm/sec) 
degree of coupling 

universal gas constant [8.13 J/(K mole), or 1.99 cal/(K mole), or 0.082 
liter-atm/(K mole)] 


, resistance to species i (ohms cm?) 


total membrane resistance (ohms cm?) 
absolute temperature (K) 

time (sec) 

mobility of species i [moles cm? /(cal sec)] 


V_ volume (cm?) 


partial molal volume of species / (cm3 /mole) 


x thickness of membrane (cm) 


valence of species / 
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B; partition coefficient of species i 

wu; chemical potential of species i (cal/mole) 

i; electrochemical potential of species i (cal/mole) 
™ osmotic pressure 

go; reflection coefficient of species i 

y electrical potential (volts) (see Figure 2.1) 


Index 


active transport, 82-3, 87, 124 
primary, 83, 106 
secondary, 83 
activity, 12 
coefficient, 12-13 
affinity, 89, 94 
anomalous osmosis, 66, 82 
antiport, see countertransport 


chemical potential, 3, 48 
chemi-osmotic coupling, 83 
coefficient 
activity, 12-13 
coupling, 16 
cross, 16 
diffusion, 24, 40 
of electro-osmotic pressure, 72 
elemental, 73 
partition, 14, 24, 133 
permeability, 24, 39-41, 63, 75, 78 
reflection, 59, 74, 78 
straight, 16 
competitive inhibition, 39, 99 
conductance 
partial ionic, 119, 123 
slope, 120 
conjugate 
driving force, 3 
property, 3 
constant-field flux equation, 25 
cotransport, 103, 124 
countertransport, 101 
coupling 
chemi-osmotic, 83 
coefficient, 16 
degree of, 91 
negative, 19, 86, 91 
osmo-osmotic, 83 
positive, 86 


diffusion 
coefficient, 24, 40 
discontinuous, 33-9, 133 


exchange, 85, 99 
facilitated, 97 
first law of, 24 
potential, 42, 46-51 
restricted, 40 
single-file, 84-5 
dissipation function, 17 
Donnan (or Gibbs-Donnan) equilibrium, 
54 
Donnan ratio, 54 


efficiency, 20, 93 
maximum, 93 
electrical conductance, 71 
resistance, 28, 117 
electrochemical potential, 8-11 
standard state, 11, 132 
electromotive force, 69, 117, 123 
electroneutrality, 42, 52, 128, 133 _ 
electro-osmosis, 70 
electro-osmotic 
permeability, 71 
pressure, 72 
electrophoresis, 70 
energy 
Gibbs free, 6, 8 
internal, 3 
entropy, 5-6 
equation 
constant-field flux, 25 
Gibbs, 6, 17 
Goldman-Hodgkin-Katz, 25, 44-5 
Henderson (diffusion), 52 
Kedem-Katchalsky, 60 
linear phenomenological, 15-17 
Nernst equilibrium, 44 
Nernst-Planck, 21-3 
Planck, 52 
Renkin, 40, 67 
of state, 1 
Ussing, 31-3, 84-7 
van’t Hoff, 58 
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equilibrium 
conditions of, 13-14 
Donnan, 54 
Nernst, 14 

exchange diffusion, 85, 99 


facilitated diffusion, see facilitated transfer 


facilitated transfer, 85, 97-9 
flux-ratio equation, see Ussing equation 


Gibbs 
equation, 6 
free energy, 6, 8 
Goldman-Hodgkin-K atz equation, 25, 
44-5 


Henderson (diffusion) equation, 52 
hydraulic conductivity, 59, 74, 79 


independence principle, 29-30, 33 
inhibition 

competitive, 39, 99 

trans-, 99 
integral resistance, 117 


level flow, 92 
mobility, definition of, 21 


negative coupling, 19, 86, 91 
negative osmosis, see anomalous osmosis 


osmosis 
anomalous, 66, 82 
electro-, 70 
osmotic pressure, 58 
effective, 60 


partial molal entropy, 10 
free-energy, 9 
volume, 10 

partition coefficient, 14, 24, 133 
of ‘‘charged”” membranes, 54-6 

permeability coefficients, 24, 39-41, 63 
of composite membranes, 75, 78 
electro-osmotic, 71 

Planck equation, 52 

potential 
chemical, 3, 48 
diffusion, 42, 46-51, 56 
electrochemical, 8-11 
interfacial, 25-44 
Nernst equilibrium, 14, 28, 44 
open-circuit, 126 
reversal, 28, 111-13 


standard electrochemical, 11, 132 
streaming, 72 
property 
conjugate, 3 
extensive, 2 
intensive, 2 
state, 1 


reciprocal relations, 18-19 
Tectification, 26, 38, 78, 120 
reflection coefficient, 59 

of composite membranes, 75, 78 
Renkin equation, 40, 67 
resistance, 28, 117 
reversal potential, 28, 111-13 
theogenic transport, 109, 124 


short-circuit current, 126 

single-file diffusion, 84-5 

slope conductance, 120 

solute drag, 86, 124 

solvent drag, 63, 65, 82, 86, 124 

static head, 91 

streaming current, 70, 72 
potential, 72, 124 

symport, see cotransport 


thermodynamics 
first law of, 3-4 
second law of, 4 
transeffects, 38, 98 
transference number, 48-50, 71, 123, 134 
transinhibition, 99 
transport 
active, 82-3, 87, 124 
theogenic, 109, 124 
co-, 103, 124 
counter-, 101 
primary active, 83, 106 
secondary active, 83 
transport number, see transference 
number 
transstimulation, 98 


ultrafiltration, 62 
Ussing equation, 31-3, 84-7 


van’t Hoff equation, 58 


work 
function, 7 
useful, 6-7, 19-20 


zero current condition, 43, 123 


This book provides an introduction to the basic principles 
employed for the analysis and description of solute and 
ments across biological membranes. 

Tt opens with a brief survey of the thermodynamic 
erning the transitions of a system from one state to another 
from the transfer of matter within the system. Subsequent chapters 
deal with specific transport processes such as simple and restricted 
diffusion of ions and nonelectrolytes, the origins of transmembrane 
electrical potential differences, osmosis, solvent drag, electrokinetic 
processes, and carrier-mediated flows. The final chapters deal with 
the coupling between solute transport and exergonic chemical reac- 
tions and the use of equivalent electrical circuit analogs for the de- 
scription of ion transport across cell membranes. 

The book presents fundamental principles in an easily readable 
monograph and provides references to the more detailed literature. 
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